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Eindimensionale Probleme 
bei einem nichtlinearen Elastizititsgesetz 


Von FRIEDRICH JINDRA, Stuttgart 


1. Einleitung 


' Die klassische Elastizitatstheorie geht vom Hookeschen Gesetz aus, das 
einen linearen Zusammenhang zwischen Spannungen und Verzerrungen an- 
‘nimmt. Dabei werden die Verzerrungen als unendlich klein behandelt. Nun gibt 
es aber auch Werkstoffe, die schon bei diesen kleinen Verzerrungen ein nicht- 
lineares Verhalten zeigen und also zur Beschreibung des elastischen Verhaltens 
ein nichtlineares Elastizitaétsgesetz erfordern. 

Hier beniitzen wir ein nichtlineares Elastizitatsgesetz fiir kleine Verzer- 
rungen, das samtliche elastischen Probleme umfasst und die klassische Theorie 
als Grenzfall enthalt. Auf ein rechtwinkliges (x, y, z)-System bezogen, lauten 
die Spannungs-Dehnungs-Beziehungen dieses nichtlinearen Gesetzes in allge- 
meiner Form}): 


1 9 
€y = 3K. R(S9) 09 te OG g(t) (» a So) z 


=~— 
— 
— 
eee 


1 2 
fs "3 R(S9) 69 + ei. g(tG) (6, — o%) » 


1 i a 1 2 
kd eae ol g(t) Tene Lee 8 (06) Tyz » Pex ~ G 8(l6) Tea - 


Dabei bedeuten oz, 6,, 0, die Normalspannungen, T,,, Ty, Tz, die Schubspan- 
mungen, &,€,, €, die Dehnungen und Yay Wyo Yzx die Scherungen, K und G den 
Kompressionsmodul und den Schubmodul, sy und fy die dimensionslosen Span- 
nungsgréssen 

1 


i So= = mit 5 — 3 (Cy + Oy + a) , 
ee 2 fp 2 
ge omit B= 3 7 [Gx — %)? + (oy — 0)? + (0% — %)*] (2) 


: + thy +13, + tel. 


1) H. Kauperer, Ing.-Arch. 17, 450 (1949). 
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Die Funktionen A(s,) und g(¢?), die als sogenannte Kompressionsfunktion und 
Schubfunktion neu eingefiihrt werden, und die im Falle des Hookeschen Ge- 
setzes den Wert Eins annehmen, werden als Potenzreihen von sy bzw. von 2 


t is q 
angesetz B(So) = 1+ Ay Sp + Ba SB + hg 8 +E ---, 1 
(3) 

g(t) =14+ gb +b +8 t-:-- | 


Die Funktionen (3) sind fiir jeden Stoff durch geeignete Versuche zu ermitteln. 5 

Hier handelt es sich um Lésungen von elastomechanischen Problemen, die 
nur von einer Koordinate abhangen. Legt man, kartesische Koordinaten zu- 
grunde, so ergibt sich sofort der einfache Zugversuch. Durch die experimentelle 
Verwirklichung dieses einfachen Spannungszustandes kann man die zwei will- 
ktirlichen Funktionen (3) des Elastizitatsgesetzes (1) bestimmen. Wesentlich) 
interessanter sind die eindimensionalen Probleme mit Zylinder- oder Kugel- 
koordinaten. Bei allen Zahlenbeispielen zeigt es sich, dass die Spannungsver-— 
teilung von den Werten der linearen Theorie erheblich abweicht, und zwar ins— 
besondere an den Randern, so dass der aus der technischen Praxis haufig be- 
kannte erhebliche Abbau der Spannungsspitzen schon durch sehr kleine Ab- 
weichungen des Elastizitatsgesetzes von der Linearitat erklart werden kann_ 


2. Der Zugversuch 


Der Zugversuch eignet sich zur Bestimmung der Funktionen (3), wenn mani 
dabei ausser der Langsdehnung auch noch die Querdehnung in Abhiingigkeit 
von der Spannung misst?). 

Ist bei einem Zugversuch von den Spannungen nur o, =o + 0, so folgt: 
nach (2) 


eo oe 
Damit erhalt man nach (1) fiir die Dehnung in Zugrichtung 


1 il o 1 2 o* | 
fa else lax) te asm (4) 
und fiir die Querdehnung | 
ae Oe i ep o 1 OR fe | 
OU araek ee toe) 2G eras he (3) 
Aus (4) und (5) ergeben sich dann die Beziehungen 
o SG Patan) 2G . 
A soe) = —— (f& +28), else) ee eee (6) 


woraus sich sowohl K und (sy) als auch G und g(¢2) ermitteln lassen. 


') H. Kauperer, Ing.-Arch. 17, 460 (1949). 
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| Bei Versuchen?) an gezogenem Kupfer ergab sich der in Figur 1 wieder- 
gegebene Verlauf der Langs- und Querdehnung abhiangig von der Zugspan- 
nung. Die Hookesche Gerade ist als Tangente der Kurven im Nullpunkt ge- 
-strichelt eingetragen. Man findet daraus 
K = 1,37-10%kg/em? und G=0,46-10% kg/cm?. 
‘Wertet man die Versuchsergebnisse nach (6) aus, so darf innerhalb der Ver- 
- suchsgenauigkeit die Kompressionsfunktion (so) konstant gleich Eins gesetzt 


kg/em?47 
1000 


on 


800 


600 


400 


200 


0 2 4 6 &-10-* 
Figur 1 
Spannungs-Dehnungsdiagramm fiir gezogenes Kupfer. 


werden, und die Schubfunktion g(¢$) kann in guter Naherung durch die ersten 
zwei Glieder der Reihenentwicklung (3) dargestellt werden, also 


R(s) =1, g()=1+4 &4 (7) 
* mit g, = 0,18-10°. 


3. Das dickwandige Rohr 


Wird ein dickwandiges, zylindrisches Rohr, dessen Lange gross gegentiber 
seinem Durchmesser ist, durch gleichmissigen Innendruck beansprucht, so 
liegt in guter Naherung ein ebenes Spannungsproblem vor. Fiir die axiale 


Spannung wird also vorausgesetzt: 
o,=0. (8) 


iner Metalle und ihy Gang im periodisch en 


1) K. Craus, Bestimmung dey Poissonschen Zahl re 
System, Dissertation (Stuttgart 1953), S. 39. 
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} 
Die Belastung ruft einen drehsymmetrischen Spannungs- und Formanderungs- | 
zustand hervor. Das Problem hangt nur von der Koordinate 7 ab. | 


Das Kraftegleichgewicht an irgendeinem aus der Zylinderwand herausge- 
schnittenen Element (Figur 2) liefert die Beziehung zwischen der Radialspan- 


Figur 2 


Dickwandiges Rohr unter Innendruck. 


nung o, und der Tangentialspannung o, 


: 
Op =O, +7 ae (9) 


Mit den Hauptspannungen o, und o, errechnet man nach (2) die invarianten 
Spannungsgréssen zu 


2 


1 P . 
S=o (rte), B= oar (Hos — a, 0). (10) 


Mittels der radialen Verschiebung wu lassen sich die Dehnungen in radialer und 
tangentialer Richtung darstellen durch 


du u 
ep = ay , Ep = y . (11) 


Der Zusammenhang der Dehnungen mit den Spannungen folgt dann aus dem 
nichtlinearen Elastizitatsgesetz (1) zu 


du ih Hl 
= Gr 9K ME Orb on) + 5G AB) 2or—o,). | 


1 ik A 
Eo = | =p ae PS) (Ont o,) + 6G &(4%) (20, —o;) . | 


Durch Entfernen von u aus den beiden Gleichungen (12) und mit (9) ergibt sich | 


5 as. 
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die Differentialgleichung fiir die Radialspannung 


Lod ; ae , 
SE Gy [elso) (20+ 7 o)] + SE Gy Lele) (or + 2 7 oy) | 


; (13) 
mit : 5 
80 = oe (2 or +1701) » = 9G: (02 + 10,0, + 77 0;2) . (14) 
Striche bedeuten Ableitimgen nach 7. 
Hierzu treten die linearen Randbedingungen: 
o,=0 ftir y=7,,, G=—P tir 7=7%; (15) 


wobei 7,, den 4usseren und 7, den inneren Halbmesser bedeutet. 
Beschrankt man sich fiir die Anwendungen auf den Sonderfall (7); SO ver- 
einfacht sich die Differentialgleichung (13) zu 


(207202 + 2A70,0,+A07 +1) 176, | 
(16) 
+ (40 P2012 + 5A70,0,+ 3102 + 3)0,=0 | 
mit der Stoffzahl Baer 
shoei (17) 
Gk +G)G 


Wir wollen hier die Integralkurve von (16) in der Umgebung der Stelle 
y = yr, durch eine Potenzreihe darstellen. Zu diesem Zweck setzen wir in (16) 


s= i= é 
Vo 
und machen fiir - 7 » 
: a(t) =6,(1— 5) = ab) 
den Ansatz = 
G(s) = DG, S* = CS + CS? + O80 +" (18) 
pal 


wobei wir das erste Glied cy = 0 setzen und somit von vornherein die erste 


Randbedingung (15) erfiillen. 
Geht man mit (18) in (16) ein, so lassen sich die Koeffizienten der Reihen- 


entwicklung durch Vergleich entsprechender Potenzen von s der Reihe nach 
berechnen : 


21421 a)e,-Bt4Ad)a=0, 
6 (1+ 218) —8%+ 2484-204) aQt BAG=%, 
12 (1+ 24 c;) Cy — 15 cy + 3A (24 Ce 0 — 28 ¢1) Os 

+ A (16 cg — 108 2 ¢ + 95 c2) cg —12AG=0, 
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20 (1+ 2A c2) cs — 24, + 4A (32 c, c — 36 Ch) ey 
+ A (72 Cy cy + 96 c3 — 412 c, cy + 173 Cf) Cg 
— A (88 c2 — 212 c, ce, + 71 ci) ce =0, 


0 (1+ 2A c?) cg — 35¢5+ 5A (40 cy c, — 44 CR) C5 
+ A (240 cs ¢, + 160 c3 — 668 c, c, + 273 ci) cy 
+ A (180 ¢3 ¢, — 372 Cz ¢, — 476 c3 + 734 Cyc, — 116 c®) cg 
+ 4 (150 ¢ — 135 c, ¢,)¢,= 0; 


42 (1+ 2A ci) cy —48c,+ 6A (48 cy c, — 52 cf) Ce 
+ A (360 cy c, + 240 c2 — 984 cy cy + 395 c2) cs 


+A (192 cy cy + 576 cz cg — 1164 cg c, — 744 c3 + 1120 c, c, — 171 C3) c, 
+ A (108 c3 — 828 cs cg + 615 cy c, + 754 c2 — 430 cy cy) cg — 841 8=0, 


36 (1 + 24 ci) cg — 63 c, + 74 (56 c, c, — 60 €) c, 
+ A (504 cg cy + 336 c2 — 1360 cy c, + 539 C2) Cy 


+ A (448 cy Cy — 888 c, c, + 504 c? — 2524 C5 cy 
+ 1834 cg cy + 1124 c2 — 622 cy cy) Cy 
— A (468 c} — 1235 cs cy + 336 cy c, + 395 c2) cg = 0, 


usw. 
Die Auflésung dieser Gleichungen ergibt: 


Cy 


a S, (3+4a), 


—ta~ (12 + 5la + 8402+ 4403), 


ss a Ne (60 + 396 a + 1312 a? + 2224 a3 + 1920 at + 640 a?) 
2 e We “y7 (360 + 2955 a + 14877 a® + 42870 a®+ 76072 at 
+ 81120 a5 + 48016 a® + 11872 a?) , 


C.= a5G Fr (2520 + 22365 a + 158325 a? + 650724 a3 
+ 1761036 a4 +3191968 a5 + 3877552 a® 
+ 3023680 a? + 1369920 a8 + 271 104 a 


+ A (560 c, ¢, + 840 cs cg — 1672 cy cy — 1068 c3 + 1582 cy c, — 236 c?) cs 


| (19) 


to 
ig 


4 


¥ 
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usw. mit den Abkiirzungen 

R= Ac, av = 14°20... ~ (19a) 


Die Lésung ist nun der zweiten Randbedingung (15) anzupassen. Man erhalt 


Delta g) = he. (20) 


womit eine Bestimmungsgleichung fiir die Konstante c, gewonnen ist. 
Die Tangentialspannung o, findet man endlich nach (9) und (18) zu 


3, (s) = H(9) — (1-9) Be Se, % 0-9) SY byes (21) 


u=0 
Fiir die beiden Ringspannungen an den Stellen y= 7, und 7 = 7% folgt 
as — = Vo \# 
GAO) ee Cy 5 a, — 79) =—py—— YF w+ )) Cu+1 (ae) (22) 


Ya Yq p=0 


Fiir 2 > 0 geht unsere Lésung (18) in die bekannte Lésung der linearen 
Theorie tiber. Nach (19) findet man namlich fiir die Koeffizienten der linearen 
Theorie 


Damit wird 


ot”) = —> of (1-74). (23) 


2 
eee a soe (24) 
mit der Abkiirzung 
o= <s ‘ (25) 


Die Spannungskomponenten der linearen Theorie nehmen dann in der Tat die 
bekannten Werte an: 


= Sees fe) ne ~ (1+ | (26) 
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samt den Randwerten 


Zz, e7+1 
ihe) Meee = Ae ? (rye = Po FE er 


Als Zahlenbeispiel wahlen wir Kupfer von Ziffer 2, fiir welches man 
A = 0,255-10-§ cm*/kg? errechnet. Figur 3 zeigt den Verlauf der Radialspan- | 


kg/cm? 4 Or, % 


Figur 3 


Spannungsverteilung in dickwandigem Rohr. 


nungen o, und o* sowie der Tangentialspannungen o, und o* langs eines Halb- 
messers (und zwar fiir 9 = 2,0 und fiir Po = 610 kg/cm?2). Die Spannungsvertei- 
lung wird infolge der Nichtlinearitit gleichmiassiger, die Spannungsspitze am 
Innenrand wird herabgesetzt. Figur 4 stellt die Abhangigkeit der Tangential- 
spannung o, am Innenrand 7 = 7, vom Innendruck Po fiir verschiedene Werte 
des Halbmesserverhiltnisses @ dar. Mit steigendem 0 weichen die Kurven 
(0,),-r, immer mehr von den Geraden der linearen Theorie ab, die gestrichelt 
eingetragen sind. Der Einfluss der Nichtlinearitat wird also bei grosseren Halb- 
messerverhaltnissen starker. 
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kg/om*4 (Ap) = ry 


1200 


Po 
0 100 200 300 400 500 600 700 kg/cm? 
Figur 4 
Anderung der Tangentialspannung am Innenrand ftir dickwandiges Rohr. 


4. Die Hohlkugel 


Legt man Kugelkoordinaten zugrunde, so ergibt sich hier die Berechnung 
der Spannungsverteilung in einer dickwandigen Hohlkugel, die durch gleich- 
massigen Innendruck oder Aussendruck beansprucht wird. Die Losung dieses 
Problems ist bereits an einer anderen Stelle!) mitgeteilt worden. Unter Be- 
schrankung auf den Sonderfall (7) gelingt es, die Differentialgleichung fiir 
die radiale Spannung geschlossen zu integrieren. Auf die Wiedergabe der For- 
meln mége hier verzichtet werden. 

In Figur 5 ist die Verteilung der Radialspannung o, und der Tangential- 
spannung o, fiir Kupfer von Ziffer 2 tiber dem Halbmesser aufgetragen, und 
zwar wieder fiir das Halbmesserverhaltnis 0 = 7,/7) = 2,0 und fiir den Innen- 


- druck #) = 1000 kg/cm*. Die Spannungswerte der linearen Theorie of und of 


sind gestrichelt eingetragen. Man sieht auch hier, dass eine kleine Abweichung 
des Elastizitatsgesetzes von der Linearitat eine recht erhebliche Verminderung 
des Wertes der Tangentialspannung am Innenrand verursacht. Figur 6 stellt 
die Abhangigkeit der Tangentialspannung o, am Innenrand v = 7, vom Innen- 
druck , graphisch dar. Die Abnahme der Tangentialspannung wird um so 
grésser, je grésser das Halbmesserverhiltnis 9 = 74/7 ist, das heisst, je grosser 
die Spannungsspitze der linearen Theorie am Innenrand ware. 


1) F, Jinpra, Ing.-Arch. 22, 411 (1954). 
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kg/cm? Op, Op 


Figur 5 
Spannungsverteilung in dickwandiger Hohlkugel. 


800 


600 


400 


200 


Po 
0 200 400 600 800 7000 kg/cm? 


Figur 6 


Anderung der Tangentialspannung am Innenrand fiir dickwandige Hohlkugel. 


@ 
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be Summary 


; Using a general nonlinear elasticity law for small strain, solutions of some 
_elastomechanic problems are stated depending on only one coordinate. If Car- 
_tesian coordinates are applied, the equations of the simple tensile test are imme- 
diately given, by which the two arbitrary functions of the nonlinear elasticity 
_law can easily be determined. In cylindrical and polar coordinates, the stresses 
are calculated in thick tubes and thick hollow spheres subjected to a uniform 
internal pressure. Numerical examples demonstrate a considerable deviation of 
stress distribution from linear theory, especially at the boundary, so that the 
decrease of stress maximum known in technical practice can be explained by 
; small deviations from linearity. 


(Eingegangen: 5. August 1954.) 


Einfache Resonanzversuche am symmetrischen Kreisel 


: Von Fritz KIRCHNER, K6ln?) 


Die wachsende Bedeutung gyromagnetischer Messmethoden in der moder- 
nen Atomphysik, von denen die eindrucksvollen Kernresonanzversuche von 
F. Biocu und seinen Mitarbeitern am bekanntesten geworden sind, liess es 
wiinschenswert erscheinen, das diesen Methoden zugrunde liegende physika- 
lische Geschehen durch einen méglichst einfachen mechanischen Modellversuch 
~ anschaulich zu demonstrieren. Der Blochsche Kernresonanzversuch besteht 
 bekanntlich in folgendem: Atomkerne, zum Beispiel Protonen, die ein bestimm- 
tes mechanisches Impulsmoment und ein magnetisches Moment besitzen, fiih- 
ren in einem verhidltnismassig starken magnetischen Feld eine Prazessions- 
- bewegung aus, deren Frequenz durch ihr «gyromagnetisches Verhaltnis» und 
die Feldstarke des magnetischen Feldes bestimmt ist. Lasst man auf die praze- 
- dierenden Atomkerne zusitzlich noch ein sehr schwaches magnetisches Wech- 
selfeld quer zu dem konstanten Magnetfeld einwirken, dann erfolgt nunmehr 
die Priizession aller Atomkerne in gleicher Phase und lasst sich infolgedessen 
makroskopisch mit Hilfe einer Induktionsspule als rotierende Magnetisierung 
nachweisen, falls die Frequenz des zusitzlichen schwachen Wechselfeldes mit 
der Priizessionsfrequenz der Atomkerne iibereinstimmt. 

Nach verschiedenen Vorversuchen hat sich die in Figur 1 abgebildete Ver- 
suchsanordnung ergeben, mit der sich diese Erscheinungen in sehr einfacher 
Weise mechanisch veranschaulichen lassen. Ein symmetrischer Kreisel, der 
durch einen eingebauten kleinen Motor mit veranderlicher Umdrehungszahl 
angetrieben wird, ist in einer in der Institutswerkstatt hergestellten karda- 
nischen Aufhangung befestigt, und zwar so, dass seine Figurenachse durch 


1) I, Physikalisches Institut der Universitat. 
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Figur 1 


Versuchsanordnung zum Studium der erzwungenen Schwingungen des symmetrischen schweren 
Kreisels. 


den Drehpunkt der Aufhangung hindurchgeht (Figur 1). Der Schwerpunkt des 
Kreisels fallt nicht mit dem Aufhingepunkt zusammen; seine Lage auf der 
Figurenachse kann aber mit Hilfe eines langs der Figurenachse verschiebbaren 
kleinen Zusatzgewichts in gewissen Grenzen verdéndert werden. Das von der 
Schwerkraft herriithrende Drehmoment entspricht dem Drehmoment, das die 
Atomkerne im konstanten magnetischen Feld des Blochschen Versuchs infolge 
ihres magnetischen Moments erfahren. Die Schwerkraft verursacht eine Pra- 
zession der Kreiselachse, deren Umlaufsgeschwindigkeit sich in bekannter 
Weise aus dem Quotienten der Betrage von Drehmoment und Impulsmoment 
ermitteln lasst (vgl. unten), die aber nach einiger Zeit infolge der Lagerreibung 
abklingt, wobei der Schwerpunkt seine tiefste Lage erreicht. 

Mit einer in horizontaler Richtung gespannten und an der Figurenachse des 
Kreisels befestigten Spiralfeder (vgl. Figur 1) lasst man nun auf den Kreisel zu- 
sdtzlich eine schwache, aber periodisch sich andernde Kraft in horizontaler 
Richtung, also senkrecht zur Schwerkraft, einwirken; diese periodisch sich 
andernde Federkraft entspricht dem magnetischen Wechselfeld des Blochschen 
Versuchs. Die periodische Spannung der Feder erfolgt mit Hilfe einer kleinen 
Exzenterscheibe, die iiber ein Getriebe von einem Motor in Umlauf versetzt 
wird; die Umlaufsfrequenz, und damit auch die Frequenz der «erregenden 


. 
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Kraft», lasst sich in weiten Grenzen andern (zwischen etwa 0,1/s und 20/s). 
Die sich unter dem Einfluss der periodisch verdnderlichen Kraft ausbildende 
_erzwungene Prazessionsbewegung der Kreiselachse lasst sich zwar am Kreisel 


selbst ganz gut verfolgen. Um sie aber einem grésseren Auditorium sichtbar zu 
-machen und um gleichzeitig auch die Phasenbeziehungen zwischen dieser 


Prazessionsbewegung und der erregenden Kraft demonstrieren zu kénnen, ist 
auf der oberen Spitze des Kreisels und auf dem Schieber, an dem die Spann- 


feder befestigt ist, je ein kleines Spiegelchen angebracht. Die beiden Spiegel- 


a, 


chen sind schwach gekriimmt (Galvanometerspiegel) und reflektieren aus dem 
Lichtkegel einer hellen Lampe je ein schmales Biindel heraus, so dass man auf 
dem Projektionsschirm in vergréssertem Mafstab den zeitlichen Ablauf der 
-Prazessionsbewegung und der erregenden Kraft verfolgen kann. 
Schon bei den Vorversuchen zeigte sich, dass ausser der Resonanzerschei- 
nung, die dem Blochschen Versuch entspricht, mit der beschriebenen Anord- 


nung auch interessante Beobachtungen tiber die Phasenverhiltnisse bei der 


erzwungenen Prazessionsbewegung gemacht werden konnen. Es erschien des- 


halb wiinschenswert, parallel mit der weiteren Ausgestaltung der Versuche 


auch die Theorie der erzwungenen Schwingungen des symmetrischen Kreisels 
nicht nur fiir den Resonanzfall, sondern auf méglichst allgemeiner Grundlage 
auszubauen. Dieser Aufgabe hat sich R. WIEBELITZ unterzogen; er berichtet 


- dariiber im einzelnen im nichsten Aufsatz. An dieser Stelle soll deshalb nur 


eine kurze Ubersicht tiber die Versuchsergebnisse gegeben werden, die man er- 
halt, wenn man, etwa mit dem kleinstméglichen Frequenzwert der erregenden 
Kraft anfangend, die Frequenz der Kraft langsam steigert. Wenn man wartet, 
bis eine etwa vorhandene Eigenprazession abgeklungen ist, so besteht die 
schliesslich allein iibrigbleibende «erzwungene Schwingung» stets darin, dass 
die Figurenachse des Kreisels auf einem elliptischen Kegel — im Resonanzfall 


auf einem Kreiskegel — prazediert. Solange die Frequenz der Erregung kleiner 
"ist als die Frequenz der «Eigenprazession», hegt die gréssere Halbachse der 


Ellipse in der Richtung der erregenden Kraft); nach dem Passieren der Reso- 
nanzstelle kehrt sich das Verhaltnis der beiden Halbachsen um, so dass nun die 
gréssere Halbachse senkrecht zur erregenden Kraft liegt. Steigert man die 
Frequenz der Erregung weiter, so wird die Ellipse immer flacher und geht 
schliesslich in eine Gerade senkrecht zur erregenden Kraft tiber. Die kritische 
Frequenz, bei der also die Figurenachse des Kreisels nur noch senkrecht zur 
Kraftrichtung hin und her pendelt, ist unabhangig von der Umlaufsfrequenz 
des Kreisels und stimmt mit der gewohnlichen Pendelfrequenz des stillstehen- 
den Kreisels iiberein?). Steigert man die Frequenz der erregenden Kraft tiber 
diese kritische Frequenz hinaus, dann beschreibt das Ende der Figurenachse 
wieder eine Ellipse, die nun aber im entgegengesetzten Sinn durchlaufen wird 


1) Bei den Versuchen rotiert der Kreisel stets in gleichem Sinne um seine Figurenachse. 
2) Vel. R. WIEBELITZ, im folgenden Aufsatz S. 362. 


| 
=> ~<a ‘| 
| 
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als vorher. Bei weiterer Steigerung der Frequenz der erregenden Kraft ergibt | 
sich ein zweiter Resonanzfall, der der «Nutation» des Kreisels entspricht und | 
bei dem die Figurenachse ihren Resonanzkreiskegel im umgekehrten Sinne | 
durchlauft als im ersten Resonanzfall. Bei noch héheren Frequenzen der er- | 
regenden Kraft geht der Prazessionskegel wieder in einen elliptischen Kegel | 


a 
e) 
e) 
ea = 
Ta G2 


Figur 2 
Maximaler Ausschlag « der Figurenachse in Richtung der erregenden Kraft als Funktion der Erre- 
gungsfrequenz (O: Messwerte; ausgezogene Kurve: nach Gleichung (6) der folgenden Arbeit von 
WIeEBELITz berechnet; Werte der Konstanten: Q, 314s", O, 854 g cm?2, ©, 6600 g cm?2, 
Gs 334000 dyn cm, ky 4,32 p, s,; 5cm). 


uber, wobei die gréssere Halbachse der Ellipse wieder, wie im Gebiet der 
niedrigsten Erregerfrequenz, in der Richtung der erregenden Kraft liegt. 

Figur 2.gibt das Ergebnis einer Versuchsreihe wieder, die WIEBELITz in der 
oben angegebenen Weise durchgefiihrt hat. Als Abszisse ist dabei die Frequenz 
der erregenden Kraft aufgetragen, als Ordinate die Amplitude des Winkelaus- 
schlags in Richtung der erregenden Kraft, den die Figurenachse des Kreisels 
bei der erzwungenen Prazessionsbewegung der betreffenden Frequenz zeigt. 
Aus der Abbildung ist ersichtlich, dass die gemessenen Werte recht gut mit den 
theoretisch berechneten (vgl. den folgenden Aufsatz von WIEBELITZ) tiberein- 
stimmen. 

An dieser Stelle soll noch auf einen hiibschen Versuch hingewiesen werden, 
der durch die theoretische Behandlung des Problems angeregt worden ist. 
Stellt man den Kreisel so auf, dass die Drehachsen der Kardanringe in der Ruhe- 
lage parallel und senkrecht zur Richtung der erregenden Kraft sind, und wahlt 


Op eA ee rN a 


a 
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man die Erregerfrequenz so, dass bei freier Beweglichkeit aller Kardanringe 
nur eine Pendelung der Figurenachse in der Ebene senkrecht zur erregenden 
Kraft, also kein Ausschlag in Richtung der Kraft stattfindet, dann erhalt man 
eine sehr ausgepragte Resonanzpendelung in Richtung der Kraft, wenn man 
mit der Hand einen Kardanring so festhalt, dass die Figurenachse in der Rich- 
tung senkrecht zur Kraft nicht mehr pendeln kann. 

Die oben beschriebenen Versuche wurden bei der Physikertagung am 
26. April 1954 in Goslar vorgefiihrt ; gleichzeitig wurden die theoretischen Er- 
gebnisse des folgenden Aufsatzes von WIEBELITz mitgeteilt. 

Bei den bisherigen Versuchen blieb die Richtung der erregenden Kraft 
konstant. Um die erzwungenen Schwingungen des Kreisels auch bei «zirkularer 
Erregung» zu studieren, wurde an der Kreiselachse eine zweite Spiralfeder be- 
festigt und tiber die gleiche Exzenterscheibe wie die erste Feder periodisch 
gespannt, aber so, dass ihre Kraft zwar auch in der Horizontalebene, aber 
senkrecht zu derjenigen der ersten Feder und mit einer Phasenverschiebung 
von + 90° auf den Kreisel einwirkt. Bei dieser Art von Erregung erhalt man 
je nach dem Umlaufsinn der Richtung der erregenden Kraft entweder nur den 
ersten oder nur den zweiten Resonanzfall, wie es ja nach dem tatsachlichen 
Bewegungsablauf bei linearer Erregung — entgegengesetzter Umlaufsinn in den 
beiden Resonanzfallen! — auch zu erwarten war. 

Die genaueren Einzelheiten uber die Abhangigkeit der Amplitude der er- 
zwungenen Kreiselschwingungen von der Frequenz der Erregung und tiber die 
Phasenbeziehungen zwischen der erzwungenen Schwingung und der erregenden 
Kraft sind aus der ausfiihrlichen Diskussion der Lésungen der Differential- 
gleichungen, durch die der Bewegungsablauf beschrieben wird, ersichtlich (vgl. 
die folgende Arbeit); dabei wird auch der Einfluss der Reibung — nur die 
Reibung in den Lagern der kardanischen Aufhingung, nicht die Reibung im 
Lager des Kreisels selbst, tritt bei den Versuchen in Erscheinung, da ja der 
Kreisel durch einen Motor angetrieben wird — naher untersucht. Im Hinblick 
auf den doch nicht ganz unbetrachtlichen mathematischen Aufwand, den die 
genauere theoretische Behandlung des Problems erfordert, mag hier aber noch 
darauf hingewiesen werden, dass man die beiden Resonanzfalle auch mit ganz 
elementaren Mitteln anschaulich interpretieren kann, wenn man sie, wie man 
es ja auch bei den gewohnlichen erzwungenen Schwingungen tun kann, als 
«entdimpfte Eigenschwingungen» auffasst, deren Dampfung eben durch die 
«erregende Kraft» gerade aufgehoben wird. Der erste Resonanzfall entspricht 
dann der gewohnlichen Prazession des schweren symmetrischen Kreisels unter 
dem alleinigen Einfluss eines konstanten Drehmoments, das hier von der 
Schwerkraft herriihrt; die Prizessionsfrequenz, die der «Larmorfrequenz» in 
der Atomphysik entspricht, ergibt sich unter der Annahme}?), dass der Dreh- 


1) Diese Annahme ist mit grosser Naherung erfiillt, wenn der Kreisel hinreichend schnell um 


seine Figurenachse rotiert. 
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impulsvektor in Richtung der Figurenachse zeigt, in bekannter Weise (vgl. 


Figur 2) aus aM | 
ay dp _ 
ES sai de Sree Saas 
Zu ; 
_ 49 _Mglsnd Mg 
2, rim NSE ee” Noes (1) 


Der zweite Resonanzfall, der der Nutation des symmetrischen Kreisels ent- 
spricht, lasst sich ebenfalls als Eigenschwingung anschaulich interpretieren. Um 
das Wesentliche dieser Art von Eigenschwingung des Kreisels mit einfachsten 
Mitteln zu erfassen, ist es zweckmissig, als schematische Darstellung ein Han- 
telmodell zu wahlen (Figur 3). Die Hantelachse stelle die Figurenachse des 
Kreisels dar; die Hantel rotiere um irgendeine Achse, die durch den Punkt D 
der Figurenachse (Drehpunkt der kardanischen Aufhaéngung) hindurchgeht. 
Dann liefern die Zentrifugalkrafte (vgl. Figur 3) ein resultierendes Drehmoment, 
das stets senkrecht zur Kreiselachse, aber bei den in den beiden Figuren 2 und 
3 angenommenen Verhdltnissen entgegengesetzt zu dem im ersten Resonanzfall 
von der Schwerkraft herrithrenden Drehmoment liegt. Das von den Zentrifugal- 
krdften herrithrende Drehmoment fiihrt demnach ceteris paribus zu einem 
Umlauf der Kreiselachse auf einem Kreiskegel im umgekehrten Drehsinn als 


D D 


Figur 3 
Zur anschaulichen Ableitung der « Eigenfrequenzen» des schweren Symmetrischen Kreisels; 
(Mt Drall um Figur-Achse); a: D Drehpunkt der kardanischen Aufhangung. b: «Hantelmodell» 
des symmetrischen Kreisels zur Veranschaulichung der Nutation: D Drehpunkt der kardanischen 
Aufhangung; mD Figurenachse des Kreisels; DO Nutationsachse; 3 Zentrifugalkraft auf m als 
Ursache fiir die «Nutation», (Die zweite Teilmasse m der Hantel ist nur der Einfachheit halber 
in die Drehachse gelegt!) 


ee 


4 
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im Falle der «Prazession» (erster Resonanzfall); dies ist die als Nutation be- 
‘kannte Bewegung des kraftefreien symmetrischen Kreisels. Die Umlaufsfre- 
-quenz lasst sich in analoger Weise hinschreiben wie im ersten Resonanzfall ; 
‘man hat dazu lediglich das von der Schwerkraft verursachte Drehmoment 
durch das Drehmoment zu ersetzen, das von den Zentrifugalkraften herriihrt. 
Fir unser Hantelmodell ergibt sich also: 


|M| = m Qy o-2lcosd 


pnd Q dp | Mt | m Q%, (2 1)? sind cosé 
N~ dt ~ |N]sind — N sin® 
oder 
' N 
Qy m (2 1)? cosé (2) 


und Q, Qy = (mgl)|(m[21]? cos), das heisst, das Produkt aus den beiden 
Resonanzfrequenzen ist unabhangig von der Kreiselfrequenz. Der im Nenner 
stehende Ausdruck m (2 /)? ist aber das Tragheitsmoment @, des Kreisels um 
eine senkrecht zur Figurenachse durch den Drehpunkt D hindurchgehende 
Achse; anderseits ist m g 1 nichts anderes als das von der Schwerkraft m g her- 
riihrende Direktionsmoment D, des physischen Pendels, dessen Schwerpunkt 
den Abstand J vom Drehpunkt D hat. Fiir kleine Ausschlage gilt bekanntlich 
fiir die Frequenz 2, des physischen Pendels: 

Dy mgl 


Me By 7m aN* 


Wir erhalten demnach auf elementarem Wege die gleiche Beziechung 
Q, Qy = Oe, 
die sich unten fiir kleine Ausschlage allgemein aus der Lésung der Differential- 
gleichung ergibt. 
Summary 


Investigations are carried out on forced oscillations of a heavy symmetrical 
gyroscope — rotating with constant angular velocity about its axis ofsymmetry —as 
a function of an external linear or elliptical sinusoidal periodic force acting in 
horizontal direction on this gyroscope. The purpose of this work was to find a 
mechanic model for the phenomenon of nuclear induction. Experiments with the 
arrangement shown in figure 1 prove that the forced vibrations left over, after 
the natural vibrations have ceased, consist in a rotation of the body‘s axis on an 
elliptic cone. It is demonstrated how sense of rotation, orientation of the principal 
axis of the ellipse and phase difference between forced vibrations and acting 
external torque vary with frequency. Two frequencies of resonance are found, 
one at the natural frequency of precession and one at that of nutation, the first 
case corresponding to nuclear induction. An elementary calculation of the two 


frequencies of resonance is given. 
(Eingegangen: 21. Mai 1954.) 
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Zut Theorie der erzwungenen Schwingungen 
des symmetrischen Kreisels 


Von RUDOLF WIEBELITZ, K6ln!) 


In der alteren Literatur hat man sich durchweg nur mit dem kraftefreien 
und dem schweren Kreisel beschaftigt, nicht aber mit Kreiselbewegungen 
unter der Einwirkung periodischer Krafte. Erst in neuester Zeit ist dieses 
Problem ausfithrlich behandelt worden, und zwar von BRAUNBEK?). Den Anlass 
dazu gaben die ganz analogen Vorginge bei der Kerninduktion?). Vorher hat 
lediglich SCHULER*) wiederholt darauf hingewiesen, dass sich selbst dieses 
Problem mit Hilfe einer von A. FOppL stammenden Naherungsmethode lésen 
liesse und dass es dabei 2 Resonanzstellen gebe. 

Die vorliegende Arbeit ist von Herrn F. KIRCHNER angeregt, um die Theorie 
zu den in der vorstehenden Arbeit besprochenen Experimenten zu liefern. Sie 
beschrankt sich daher auf die Bewegungen des schweren symmetrischen Krei- 
sels, der mit konstanter Geschwindigkeit um seine Figurenachse rotiert und 
auf den senkrecht zur Schwerkraft sinusf6rmige periodische Krafte einwirken. 
Da die Untersuchungen von BRAUNBEK, soweit sie sich auf das vorliegende 
Problem beziehen (Fall IT bei BRAUNBER), sich auf Nadherungen fiir den spe- 
ziellen Fall des »schnellen» Kreisels ohne Beriicksichtigung der Reibung und 
ohne Berticksichtigung der Nutation beschranken, diirfte es von Interesse sein, 
das Problem noch einmal mit anderen Methoden und in allgemeinerer Form zu 
behandeln. 


I. Beschrankung auf kleine Winkel 
a) Das Differentialgleichungssystem 


In dieser Arbeit soll ein sehr anschauliches Gleichungssystem benutzt wer- 
den, das von A. Fépr stammt®: Mit Hilfe der Haupttraigheitsmomente @, und 


1) I. Physikalisches Institut der Universitat. 

*) W. BraunBek, Der symmetrische Kreisel mit zeitlich periodischem Richtmoment, ZAMM 33, 
174-188 (1953). 

3) Siehe vor allem F. Biocn, W. W. Hansen und M. PACKARD, Nuclear Induction und The 
Nuclear Induction Experiment, Phys. Rey. 70, 460-485 (1946). Ferner F. Brocu and A. SIEGERT, 
Phys. Rev. 67), 522-527 (1940) und R. K. Wanesness and F. Biocu, Phys. Rev. 89, 728-739 (1958). 

4) M. SCHULER, Kreisellehre (6. Kap. in Miiller-Pouillets Lehrbuch der Physik, 11. Aufl. 1929, 
LI, Teil, 1. Band, S. 763-767), und M. Scuuter, Einfiihrung in die Mechantk, Teil II (1951), S. 
157-159. 

5) A. Féppr, Vorlesungen tiber technische Mechanik, Bd. 4, 10. Aufl. (R. Oldenburg, Miinchen 
und Berlin 1944), S. 226-233; Bd. 6,6. Aufl. (R. Oldenburg, Miinchen und Berlin 1944), S. 181-182. 
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2) ', (|| parallel zur Figurenachse und | senkrecht dazu) lasst sich der Dreh- 
impulsvektor Jt folgendermassen darstellen 


4 N= 2 Oin+ 2, On, = 2,O,n+, |n, ay) (1) 


(2,, 2) = Winkelgeschwindigkeitskomponenten, n = Einheitsvektor in Rich- 
tung der Figurenachse, n, derjenige senkrecht dazu). 
_ Die dusseren Drehmomente sollen in einem rechtwinkligen Koordinaten- 
system (x, y, z) mit senkrechter z-Achse formuliert werden, dessen Mittelpunkt 
‘im Drehpunkt liegt. (Einheitsvektoren i, j, f). Zur Schwerkraft —m gf (m 
— Kreiselmasse) tritt dabei eine sinusférmige Kraft, die in der Form h, j sinw¢ 
+ k,icoswt angesetzt wird. Das Vorzeichen (--) gibt dabei den Umlaufsinn 
des Kraftvektors um die Figurenachse an. — Dabei sind folgende Sonderfalle 
zu unterscheiden: 1. lineare aussere Kraft (zum Beispiel k, = 0), 2. zirkulare 
“Kraft (k, = ks) und 3. elliptische Kraft (h,,2 beliebig). Die zu diesen Kraften 
gehorigen Radiusvektoren seien: —s n(|s| = Schwerpunktsabstand)?), s,1 
‘und s,n. Dann gilt unter Vernachlassigung der Reibung nach dem Dreh- 
impulssatz 


d*n 


. d 
Q) 9 = ae [n, aE 


= +k, s,[n, i]cosmt+ k, s,[n, j]sinwt (2 
2 


+mes{n, f]. 


Dabei ist das Glied (dQ,/dt) 0, n weggelassen worden, da im Experiment 2, 
durch den eingebauten Kreiselmotor konstant gehalten wurde’). 

Division durch 9, und Komponentenzerlegung ergibt das nichtlineare Dif- 
ferentialgleichungssystem fiir die Komponenten 1,, ”, und ”, von 1 


N, - Aysy 


s Fé Wes of 
Be — Me My + My a — 9 Thy = oe n,sinwt, 

N, - fe te 2 Ro Ss 
GM ta He — Notte + Oy Me = Oe N, COSM@t, (3) 
Ni i 55 — kysy : Ry Sg 

fas My Mo era. Ny, Sin w t -- ron Ny COS @ t 


(N, = 2,9, = Drehimpulskomponente um die Figurenachse, 27 = [mg s]/O,). 
Q, ist fiir s > 0?) gerade die Kreisfrequenz des Kreisels fiir den Fall, dass er, 
ohne zu rotieren, einfach als physikalisches Pendel schwingt. 


1) Denn dn/dt = (2, n), also [n, dn/at] = [n, (2, n] = 2 — Qyn. i. 

2) Wenn der Schwerpunkt unterhalb des Drehpunktes liegt, ist s also positiv, falls n «nach 
oben» zeigt (n, > 0). — ai. 

3) Ohne Beriicksichtigung der Reibung ware es unabhangig davon = 0, wie man durch skalare 
Multiplikation von (2) mit 1 erkennt. 


a | 
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Dieses System vereinfacht sich wesentlich, wenn man sich auf «kleine Win- 
kel», das heisst auf kleine Abweichungen der Figurenachse von der Vertikalen: 
beschrankt. Dann weicht namlich die Kugeloberflache, auf der sich die Spitze 
von n bewegt, nur sehr wenig von der Ebene z = 1 ab, und man kann nahe- 
rungsweise 1, ~ 1 setzen. Nach dem Vorbild von SCHULER sollen hier ferner 
n, und 1, durch die Drehwinkel « und # der Figurenachse um die (+ x)-Achse: 
bzw. um die (+ y)-Achse ersetzt werden. Legt man dabei fest, dass in der: 


(x z)-Ebene « = 0 ist und in der (y z)-Ebene B = 0, so folgt 
Mm, ~sinpBy B und n, y —sina Y —a.- 


| Durch Einsetzen in (3) ergibt sich dann das System mit konstanten Koeffi- 
zienten: 


Ny op Masry & as gp! Nes 2 yee 
On p= On sinwot; B+ 228 6. c= = : coswmt. | 
(4) 


Nach SCHULER?) ist es ferner zweckmAssig, statt @,; zwei verschiedene Haupt- 
tragheitsmomente ©, und @, einzufiihren®). 0, zum Beispiel ist dabei das Trag- 
heitsmoment des Kreisels um die #-Achse fiir den Fall, dass die Figurenachse auf 
die (y z)-Ebene beschrankt ist. 

Genau genommen muss ferner noch die Reibung beriicksichtigt werden, die 
vor allem in den Lagern des kardanischen Gehanges zu suchen ist, wahrend die 
Luftreibung demgegeniiber klein ist. Setzt man die Lagerreibung proportional zur 
Winkelgeschwindigkeit an, so kommen auf der rechten Seite von (2) die Glieder 


a+ Q2a%+4 


2 
— LR, ([n, 1), ty) 1, hinzu (n,, = Einheitsvektoren in Richtung der Kardan- 
m=1 


achsen). Bei Beschrankung auf kleine Winkel ergeben sich daraus auf der linken 
Seite von (4) die Zusatzglieder 


Tie se : Tice ; 
—— % = 0,0 bzw. - P= 036. 
0, 1 0, B 2 
Unter Beniitzung der Abkiirzungen 
GS BSL pia giEe Saas pun ae 
aE Ox,y nig On, oad Oxy 


wobei die jeweils ersten Indizes zusammengehéren (und ebenso die jeweils zweiten) 
erhalt man aus den allgemeinen Gleichungen SCHULERS ty 


a+ 22 a+ 0, %+ a, B = —b, sinwt; B+228+0,8—ai= + bycoswt. (5) 


1) M. Scuuter, Kreisellehre (6. Kap. in Miller-Pouillets Lehrbuch der Physik, 11. Aufl. 1929, 
1. Teil, 1. Band, S. 763-767), und M. SCHULER, Einfithrung in die Mechantk, Teil 11 (1951), S. 
157-159. 
: 2) 0, und oF unterscheiden sich in den Tragheitsmomenten der jeweils mitbewegten Kardan- 
tinge und eventuell durch verschiedene Lagen der betreffenden Drehachsen (Abbildungen bei 
SCHULER). 


3) S, und s,, = Abstande des Schwerpunkts von den Drehachsen, meist ist Sy = Sy = S. 
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b) Lineare periodische Kraft ohne Beriicksichtigung dey Reibung 


e Die Gleichungen (4) stellen einen besonderen Fall der Differentialgleichun- 
gen fiir gekoppelte Schwingungen dar. Daher haben auch die Losungen grosse 
Ahnlichkeit mit denjenigen bei Schwingungen. Fiir k, = 0 lauten sie 


a =Csinwt+ Asin(Q,¢+ gy) + Bsin(Q,¢+ 9) | 

6 
B = Dcoswt + A sign (Q) 22) cos(Q,¢ + y) — B sign (Q) cos(Q,,¢+ y) | 
‘mit 

oe 25 — w? As. p NO, © me 
| ae Nee Oy (Q2— w)?— (NP/O}) wo Oy ” 


(sign = Vorzeichen, also signa = +1 falls'a20;"2, = Kreisfrequenz der 
_ Prazession, Q,, = diejenige der Nutation; A, B, y und y sind die Integrations- 
konstanten). 
| Da die Eigenschwingungen — wie noch zu zeigen ist — durch die Reibung 
-gedampft und daher nach einer gewissen Zeit praktisch abgeklungen sind 
- (wenigstens fiir s > 0), ist es experimentell leicht méglich, die eigenschwin- 
-gungsfreie erzwungene Schwingung [= stationire Lésung; erstes Glied der 
rechten Seite von (5)] zu untersuchen. Die Figurenachse des Kreisels bewegt 
sich dabei im allgemeinen auf einem elliptischen Kegel, ihre Spitze also auf 
einer Ellipse. Dieser elliptischen Bewegung werden durch die Eigenschwin- 
- gungen, namlich durch die Prazession (2. Glied) und die Nutation (3. Glied), 
zwei kreisformige Bewegungen iiberlagert. Diese Schwingungen sind bei der 
Fépplschen Naherung alle unabhangig voneinander. —- Wenn der Nenner der 
stationairen Loésung verschwindet, liegt Resonanz vor. Das ist gerade bei den 
Eigenfrequenzen 


N, NG site aes 
ai nog ae Nal 4 VRE + OF} 


der Fall. Aus diesen Formeln ersieht man, dass 2, Q,, = Q2 ist. Fiir schnelle 
Umdrehungen des Kreisels um seine Figurenachse, das heisst fiir grosse Werte 
von Q, und N, , wird naherungsweise 
220 mgs N 2, | O 
20, mes. gy IN _ 11 


0, Sess = - 
z |Ni | | 2, | Oy 


Aus Q) > ce folgt also auch Q,, > oo. Das wird der Grund sein, weshalb bei 
BRAUNBEK keine Resonanz mit der Nutation erwahnt wird. 

Fiir den Fall, dass der Schwerpunkt oberhalb des Drehpunktes liegt (Q?< 0), 
kénnen die beiden Resonanzstellen zusammenfallen [fiir (N,/O,)? = 4 |23]] 


: 


| 
| 
366 RubDOLF WIEBELITZ ZAMP | 


oder ganz verschwinden [fiir (N\/O,)?<4 |Q2|]. Doch ist in diesem Falle | 
die experimentelle Untersuchung schwierig, da Prazessions-Eigenschwingungen 
kaum zu vermeiden sind, und diese — wie noch zu zeigen ist — unter dem Einfluss 
der Reibung exponentiell anwachsen. { 
Fiir die stationare Lésung ergibt sich ferner aus den Amplituden |C| und 
|D| von « und f, dass die Hauptachse der Ellipse fiir |w| <2, und |@| >, 
in y-Richtung, also in Richtung der dusseren periodischen Kraft liegt, dagegen’ 


GD 


Sis eee 


Figur 1 


Die Amplitudenfunktionen C (Goer ) und D( ) der stationadren Schwingung in Abhangigkeit 
von der Kreisfrequenz w der dusseren periodischen Kraft nach Formel (6). (Werte der Konstanten: 
0), = 10sec, N, /O,, = 10sec, kys,/Q2) = 1sec—®.) 


fir 2,<|o|<, quer dazu in «-Richtung. Der Umlaufsinn der Figuren- 
achse um die (+ z)-Achse ergibt sich aus den Vorzeichen. Er ist positiv oder 
negativ, je nachdem sign { Q) (w? — Q2) }=+1 ist. Fir |o| <Q, stimmt 
er mit dem Umlaufsinn der freien Prazession und fiir |w| > 2, mit dem- 
jenigen der freien Nutation liberein, wie aus den Vorzeichen der Eigenschwin- 
gungen ersichtlich ist. Dazwischen liegt der Fall der Querschwingung (« = Oy; 
der erstaunlicherweise gerade fiir |o| = Q, auftritt. Die Phasendifferenz zwi- 
schen dusserer periodischer Kraft und a-Schwingung, die experimentell eben- 
falls leicht zu untersuchen ist, ergibt sich auch aus den Vorzeichen. Sie ist 
0 oder +, je nachdem sign(—C) = +1 ist. An den Resonanzstellen und bei 
der Querschwingung gibt es daher Phasenspriinge. 

Den quantitativen Verlauf ersieht man aus Figur 1. Dabei geben die ab- 
soluten Betrage von C und D die Amplituden der betreffenden Winkel an, 
wahrend die Vorzeichen die Phasenverhiltnisse ausdriicken. Vorzeichenwech- 
sel bedeuten also Phasenspriinge. 


es 
(a 
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es Uber die in der Arbeit von BRAUNBEK behandelten subharmonischen Re- 
sonanzen lasst sich auf Grund der hier besprochenen Naherung noch nichts 
-aussagen, da diese nur die Grundschwingung erfasst, nicht die Oberschwin- 


-gungen. 
’ c) Elliptische periodische Kraft ohne Beriichsichtigung der Reibung 


Fiir 0; = @. = 0 lauten die stationdren Lésungen*) von (5) 


, b(Q- wo) Fao. 
oat): ; p Sl t 

4 1 (Q2— o) (Q2— o%}) — a, 4, w no 

‘ b, a, w F by (QZ — w?) 

= (22 — w) (Q2 — w) — a, a, © cosmt. | 


Die stationare Schwingung lasst sich also durch Uberlagerung der beiden Schwin- 
gungen erklaren, die durch die beiden periodischen Krafte erzeugt werden, wenn 
sie einzeln wirken. Diese Einzelschwingungen entsprechen genau dem schon be- 
handelten Fall (Abschnitt Ib), nur dass fiir Q) + Q5 und a, + a, Unsymmetrien 
hineinkommen. 
Beispiele: 1. 0< },<0,; Q)>0; Q2 = Q} > 0; a, = a,> 0; positiver Um- 
laufsinn (oberes Vorzeichen). Die Resonanzstellen, das heisst die Nullstellen des 
Nenners, sind an den alten Stellen geblieben, doch die Querschwingungen haben 
sich verschoben. Es gibt jetzt zwei Arten von Querschwingungen: Entweder 
wird «= 0 (an den Stellen w = Qa) oder B = 0 (fiir w = Qg). Wenn 8, von 0 
bis b, wachst (fiir festes b,), so fallt Qzy — wie man leicht nachrechnet — monoton 
'von 2, bis 2,, wahrend 2g monoton von 0 bis 2, wachst. 
: 2. Zirkulare Aussere Kraft (b,; = b2), 23 = Q7; 4, = 42. Die «Querschwingun- 
gen» fallen mit gewissen Resonanzstellen zusammen, so dass man kiirzen kann: 


by } b, 
., 2 ee ee . 2 sails, oe See es ee : 7 
Oy Q: . ee Simo er? py = ss 2 s : cos w t (7) 


Es gibt also keine Querschwingungen mehr; die Kreiselspitze bewegt sich viel- 
_ mehr durchweg auf einem Kreis. Wie man leicht nachrechnen kann, tritt dann 

und nur dann Resonanz fiir wo = +2,,, ein, wenn die entsprechende Eigen- 

schwingung denselben Umlaufsinn hat wie die aussere periodische Kraft. 


d) Elliptische periodische Kraft wnter Beriicksichtigung der Reibung 


Die allgemeine Lésung des Gleichungssystems (5) lautet: 


Bs sin(w t+ 6) + Ape~“? sin(Q,¢+ 9) + Ae e~4ntsin(Q,¢+ y), 


3 


p= ie sin(wé+ 6) + Bye“? sin(Qy¢+ e+ m) + B,e 4”! sin(Q, t+ y+na) - 
3 

W, ist eine nach unten beschrankte Funktion, so dass die Amplituden an den 

Resonanzstellen nicht unendlich gross werden kénnen. W, und W, hangen von 

b, und b, ab und konnen fiir sehr kleine Werte von b,,,s0 klein werden, dass selbst 

die Resonanzstellen in den Bereich «kleiner Winkel » fallen. 


1) Die dazugehoérigen Eigenschwingungen sind von Scnuter (Kreisellehre) diskutiert. 
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Angaben zur stationaren Lésung: 


W, = VU(2j — 0) b, F a, @ by)? + @3 wb, 
es Vi(22 — w) by F a, w bi\? + Qf w? BF, 
W, =P? + Q2 + 4a, ay 0; 0, oF 


Qe = (22 = 0) 5 (22 — 0) 1 0, 


mit 


P= (22 — @?) (23 — w?) + 01 0, w? = a, a, w?. . 


Die Phasenwinkel sind bestimmt durch 


: 1 Seat oe 
sind = W, W, {i(25 = ©)? 0, + 0, 03 w? + ay ay Q, 07) why F QO. a, o ci 


(8) 
sine= 4 { Pa, w b, F ([( 22 — w?)2 + 9? w*] (22 — w?) | 
po — a, a, (Q2 — w) w*) by}. 

- cosé erhalt man aus dem Ausdruck fiir sine und ebenso + cose aus sin 0, 

wenn man alle Indizes 1 durch 2 sowie ¥ durch y ersetzt, und umgekehrt. (Also 

W, <—> W,, 2,<—> Q,, b, <> by usw.) 
Fiir den Spezialfall 2, = 2), 4, = 4, 0, = Q, und b, = 0 sind 6 und e in 

Figur 2 graphisch dargestellt. Die geraden Streckenziige entsprechen den Phasen 


Figur 2 


Die Phasenwinkel 6 (~~-) und ¢ ( 


gezogenen geraden Linien deuten di 
Konstanten: Q, = Q, 


) sowie ihre Differenz ¢e— 6 (sss ) nach Formel (8). Die aus- 

e Phasen ohne Beriicksichtigung der Reibung an. (Werte der 

= 10sec™, a, = a = 10sec, 01 = @2 = 1sec™, db, = 0, also b, beliebig 
+ 0.) 


ohne Berticksichtigung der Reibung!). Wie bei der harmonischen Schwingung 
bewirkt die Reibung also auch hier, dass keine Phasenspriinge auftreten. Zwischen 
dem gréssten und dem kleinsten Wert von « liegt eine Differenz von 2 z. Deshalb 
war die Angabe zweier Winkelfunktionen (sine und cose) erforderlich, wahrend 
man bei den harmonischen Schwingungen mit einer Angabe (meist tge) aus- 


1) Zu 6 gehort zum Beis 


Piel der Streckenzug (0,2), (Q,,,7), (Q,,, 0), (Q »9), (2,20), (Qy, 70), 
(Q,,, 0), (00, 0). » D a Fa ia 


ae 
ee 
Vol. VI,1955 Zur Theorie der erzwungenen Schwingungen des symmetrischen Kreisels 3609 


Ve 


kommt. Da durchweg ¢ — 6 + 41/2 z ist, liegt die Ellipse, auf der sich die Spitze 
4 der Figurenachse bewegt, schief?). 
: Angaben zu den Eigenschwingungen: 

Die Eigenfrequenzen 2, und 2, sind (wie bei den harmonischen Schwin- 
_ gungen) gegeniiber dem reibungsfreien Fall verschoben. Ausserdem treten jetzt 
_ Dampfungsexponenten A, und A, auf. Die Werte %,.= —A, +19, und 
)%,,——A, +7 2, sind die Lésungen der charakteristischen Gleichung 
ct (01 + Qs) #8 + (a, ay + 01 Qo + QF + QF) x? + (0, QF + Os Q3) H+ QF QY= 0. 


_ Natiirlich lassen sich die Lésungen dieser Gleichung ohne weiteres explizit hin- 
_ schreiben, aber im allgemeinen Fall nur in recht untibersichtlicher Form. 

; Die Gréssen B,, und 1 ,, sind keine Integrationskonstanten, sondern sind 
- bestimmt durch B,, = Aj, |Up, | bzw. m,2 = arg (Up, n), wobei 


U. as (—Apy.n 4% ery) 
age (Aan aU; a fg ah =i Q2 (Ap, a a 25.) a 25 


ist. (Die jeweils ersten Indizes gehéren zusammen.) Im allgemeinen bewegt sich 
_ also die Spitze der Figurenachse bei den Eigenschwingungen nicht mehr auf einem 
_ Kreis, sondern auf einer schiefen Ellipse. Nur fiir Q? = Q?2, ay = ay, 01 = OQ. wird 
&|U,.,|=1, also By.= A;,2, und m= 1/2" sign (Q, 22), n, = —1/2 a sign (Q)), 
so dass ahnliche Verhaltnisse vorliegen wie ohne Beriicksichtigung der Reibung. 

Dies ldsst sich leicht sehen, wenn man die ¥,; (i=1,..., 4) ausrechnet, was 
allgemein fiir Q2 = 2? durchgefiihrt werden soll. Durch die Substitutionen 
x= Q22y,2= y+ (1/y) wird die charakteristische Gleichung auf eine quadratische 
Gleichung zuriickgefiihrt, aus welcher sich folgende Lésung ergibt 


Qi+ @ &y 
i oy ei ar Ver 


é 


+ 2V2 (of + of — 2 41 ay — 8 23) — 26 (ai + 2) VT 


mit 

T = (01 — 2)? — 441%; S2= qe 1. 
Fiir T => 0 gibt es je nach Grosse des Radikanden der letzten Wurzel analog zur 
harmonischen Schwingung einen Schwing-, Grenz- oder Kriechfall. Fiir T< 0 
formt man besser etwas um und erhalt 


a paw Bt, Bea aS +i VCE I BEV Rs OVS 


R= 8 03 + 2 4,04, — of — O83 
und : : ‘ 
S = (422 + aya, + 01 G2)? — 4 23 (01 + @2)*- 


Fiir den Dampfungsexponenten A, erhalt man daraus 


Qi + Q2 Lowe Se > Ort Oo 
4,5 zV R+2ysz2%4 
1 Y—R + 2V(4Q2 + a, 4, + yz= 0 
eee ae ( ave ee) HD OT.03)— ss 
1) Abgesehen wor der unmittelbaren Umgebung von | o| = Q, ist das jedoch normalerweise 


im Experiment kaum zu merken, da der Reibungseinfluss meist viel geringer ist, als in Figur 2 
(der Ubersichtlichkeit halber) angenommen wurde. 


ZAMP V1/24 


’ 
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je nachdem 22 2 0. Je nach der Lage des Schwerpunktes vergréssert oder ver- | 
kleinert die Prazessionsreibung also die Amplitude der tiberlagerten Prazession — 


exponentiell (vgl. Abschnitt Ib). 


e) Festhalten eines Kardanringes 


Aus dem System (5) lasst sich ferner die Kreiselbewegung leicht fiir den 
Fall berechnen, dass man den adusseren Kardanring festhalt, so dass sich der 
Kreisel nur noch in einer Richtung bewegen kann. (Er hat dann also nur noch 
einen Freiheitsgrad.) Diese Richtung wird einfachheitshalber als Richtung 
einer Achse, zum Beispiel der x-Achse, angenommen. Dieser Fall entspricht 
dem Grenziibergang 23 > oo. Denn bei diesem Grenziibergang wird ja das 
konstante Drehmoment um die y-Achse so gross, dass sich der Kreisel nicht 
mehr um diese bewegen kann. Dividiert man die zweite Gleichung von (5) durch 
2; und fiihrt den Grenziibergang durch, so folgt 6 = 0. Aus der ersten Gleichung 
aber ergibt sich die Gleichung der erzwungenen harmonischen Schwingung 
a+ 220+ 0,%=—b, sinw t. Wegen der Beschrankung auf kleine Winkel ist 
dies gleichzeitig die Gleichung der erzwungenen Pendelschwingung des Kreisels 
fur den Fall, dass dieser gar nicht rotiert. Fiir lineare periodische Kraft (b, = 0) 
an der Stelle der Querschwingung erhalt man daraus unter Vernachlassigung 
der Reibung folgendes interessante Ergebnis: Verhindert man die Querschwin- 
gung durch Festhalten des ausseren Kardanrings, so geht die Kreiselbewegung 
in eine Langsschwingung iiber, die sich sogar zur Resonanz aufschaukelt. 


II. Exakte Behandlung der zirkularen periodischen Kraft ohne 
Bericksichtigung der Reibung 


Die exakte Lésung des Systems (3) bereitet im allgemeinen grosse Schwie- 
rigkeiten. Fiir zirkulare dussere periodische Kraft (ky = ky, = k; s) = sy =) 
erhalt man jedoch eine einfache stationire Lésung: 


i : Co ks ; 

ate DiC at Ole oo" | : 
Cy pied (9) 

n —— SINO 43) Gy 


oe: 6 a ae (Nj/O,) oO, 


Die Kreiselspitze beschreibt hierbei also durchweg einen Kreis. - Obwohl Ce 
eine Integrationskonstante des homogenen Systems (3) ist, kann es nicht belie- 
big gewahlt werden, sondern muss aus der Zusatzbedingung 


2S a2 2 ae 
||? = m2 + m2 + m2 = 1 


berechnet werden, die besagt, dass n ein Einheitsvektor ist. Daraus folgt die 


: 


J 
= 
“a 
% 
~ 
c 


I 
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Gleichung vierten Grades fiir Cy 


(Co a = (1- C3) (@- Cott =e 0). 


iese lasst sich wesentlich einfacher nach w auflésen als nach C: 


¥] 
' rin, N 
O=-- 76,0, ° *\) 6c of 


QO» : in S 
Co YO yrs C 


+ (10a) 


j (e:,. = +1). Eine mit Hilfe dieser Gleichung erhaltene Kurve zeigt Figur 3. 
Zu jedem Werte von w gibt es demnach 2 bzw. 4 Werte von Co, das heisst, es 
_ gibt ebensoviele stationare Lésungen. [Denn wenn man C, kennt, berechnet 


man nach (9) auch leicht ”, und m,.] — Die Zweige der Kurve in der Nahe von 


‘eee 


Figur 3 
) der stationaren Losung in Abhangigkeit von der Kreisfrequenz 
) und zum Vergleich dazu nach (10b) der konstante 
) in Abhangigkeit von der Kreisfrequenz 
Q, bzw. Q,, der betreffenden Eigenschwingung. (Werte der Konstanten: N,/O-= 10 sec, 
2 = 10 sec-2, ks/O, = 10sec~*.) 


Der konstante Wert n, = Co ( 
@ der zirkularen dusseren Kraft nach Formel (10a 
Wert ”, = C der freien Prazession oder Nutation (--------~ 


C, = 1 entsprechen der Naherungslésung (7). Dazu kommen jetzt noch die 
Resonanzkurven. Sie sind, wie es auch bei andern nichtlinearen erzwungenen 
Schwingungen der Fall ist, nach grosseren oder kleineren Frequenzen hin ge- 
neigt. Zwischen den beiden benachbarten Kurvendsten herrscht dabei eine 
Phasendifferenz von 180°, die im reibungsfreien Fall nicht tiberbriickt werden 
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kann, wenn man die Eigenschwingungen ausschaltet (was durch geeignete 
Wahl der Anfangsbedingungen méglich ist). Unter dieser Voraussetzung muss 
sich daher der Kreisel beim Durchgang durch die Prazessionsresonanz umkeh- 
ren. Wenn die Spitze vorher nach oben zeigte, so zeigt sie nachher nach unten. — 

Zur Erlauterung der stationdren Lésung sollen noch die freien Eigenfre- 
quenzen’), also die Lésungen des Systems (3), fiir k; = kg = 0 angegeben wer- 
den). Fiir die Prazession lauten sie: 


Ny = + V1 — C2cos (Q,¢+ 6); 
nm, = — V1 — C3 sign (Q) 22) sin(Q,t+ 8); z= Cy. 


Die Gleichungen fiir die Nutation erhalt man daraus, wenn man Q, durch Q,, 
und die Konstanten C, und 6 durch andere Konstanten C,, und ¢ ersetzt. Die | 
Eigenfrequenzen Q,,,, sind dabei keine absoluten Konstanten, sondern hangen 
VOU Ny, = C4 abe 


xe N, i / Ni or 
25,m = lete|+ 402,01. ° C 


(Das obere Vorzeichen und C, gehdren zu Q,.) — Fir grosse Werte von 2, 
(gross im Vergleich zu den andern Gréssen) erhalt man daraus unter Beriick- 
sichtigung von C,, = cos die Formeln (1) und (2) des vorigen Aufsatzes von 
F. Krrcuner. — Tragt man C,,,, als Funktion von Q,,, auf (siehe Figur 3), 
so liegen die erhaltenen Kurven gerade zwischen den Resonanzkurven der sta- 
tiondren Lésung. Der seltsame Verlauf der Resonanzkurven ist also durch die 
Eigenschwingungen bedingt. 


(10) 


Dn | 


III. Zur Braunbekschen Naherungslésung 


Fiir den Fall der linearen oder elliptischen periodischen Kraft bietet das 
System von F6ppr grosse Schwierigkeiten, so dass man hier auf die von 
BRAUNBEK angegebene Naherungsmethode angewiesen ist?). 

Nach BRavunBEX ist es fraglich, ob die von ihm angegebenen Reihen kon- 
vergieren. Durch eine Anderung der Rekursionsmethode ist es jedoch méglich, 
die Konvergenz unter geeigneten Bedingungen zu beweisen und gleichzeitig 
die Prazessionseigenschwingungen und die elliptische periodische Kraft mit zu 
erfassen. In Analogie zu BRAUNBEK soll in Gleichung (1) fiir grosse Werte von 


1) Wegen der Nichtlinearitat kénnen diese Lésungen nicht einfach zu der partikularen Lésung 
addiert werden, um die allgemeine Lésung zu erhalten. 

2) Vgl. A. Forpt, Vorlesungen tibey technische Mechanik, Bd. 6 (R. Oldenburg, Miinchen und 
Berlin 1944), S. 230ff. 


3) Ww. BRAUNBEK, Der symmetrische Kreisel mit zeitlich periodischem Richtmoment, ZAMM 33, 
174-188 (1953). 


mi 


ay 
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_ Q) das 2. Glied vernachlassigt werden, so dass n ~ Jt/(Q, @,) wird. Setzt man 
dies in (3) ein, dividiert durch N\/@, und setzt 


Rk 2 
wa ba ee ee Bas 0. = 6, 


mgs Ny Rg Sy 


so erhalt man das Differentialgleichungssystem 


N.— Q,N,=—af Q,N,sinwt, No Q, Nz, = +4 Q,N, coset, | 


(11) 
| N,=«f8 2,N,sinwt— aN, 2,coswt | 


| (BRAUNBEK schreibt P,,... statt N,,... und mz, statt Q,). Macht man den 
Ansatz (46) von BRAUNBEK}), so erhalt man die Gleichungssysteme 


2 
in = B Qn bn1S8inwt F Qyqn1coswt (n=1,2,...). 


Fiir 2 = 0 verschwindet die rechte Seite, und man erhalt 
p= A,cos(2,¢+ 9), %= —A,sin(Q,¢+ 6), %=Co- 


Entsprechend sieht allgemein die Lésung des homogenen Systems von (12) 
aus. Variation der Konstanten ergibt dann die Rekursionsformeln 


t 
Ds @, | rea{-B sinw x cosQ, (f — x) 4 cosw xsinQ, (¢ — x) }dx 
° + A, cos(Q,¢+ d,); 


t 
Ga 2, | ma{B sinw x sinQ, (t — x) + cosw x cos M2, (¢ — x) }dx (13) 
‘ — A, sin(2, ¢ + 0,)3 


t 
t= 2, | {B bn SiN@ %F qn cosw x}dx+ Cy. 
0 


Die Integrationskonstanten 4,,, 6, und C,, sollen so gewahlt werden”), dass auf 

den rechten Seiten keine Grdssen (Q, ¢ + 6) bzw. in der letzten Gleichung 

keine Konstanten mehr vorkommen, falls n > 0 ist. - Die ersten Glieder der 
1) Entwicklung nach Potenzen von a: Ng = Z py a"; Ny = 2 qo"; Nz = 2 ry 


2) Diese Festsetzung bedeutet keine Einschrankung, denn die Gréssen, die jetzt wegfallen, 
wiirden lediglich einer Entwicklung von Ag bzw. Cy in Potenzreihen nach « entsprechen. 
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Entwicklung erhalt man hieraus leicht durch Einsetzen, zum Beispiel 


bi = C, 2, Sar gr Sree 


Ay 22 CG oe 2, 
pos ; es : . = = ae = cos[(2 w + 2,)¢+ 6] 


—Q 2 
a8 fey : 3 = = = ae cos| (2 m — Q,) t — oh. 


Hierbei treten ausser den schon von W. BRAUNBEK behandelten Oberschwin- 
gungen also noch Summen- und Differenzschwingungen zwischen erzwungener 
Schwingung und Prazessionsschwingung auf. (Ob entsprechende Glieder auch 
bei der Nutation eine Rolle spielen, kann hier nicht entschieden werden, da die 
Braunbeksche Naherung ja grundsatzlich auf die Diskussion der Nutation ver- 
zichtet.) Diese Erscheinung ist zu vergleichen mit den Summations- und Dif- 
ferenzténen, die beim Auftreten zweier Schallwellen auf einen nicht-harmoni- 
schen Oszillator (zum Beispiel aufs Trommelfell) entstehen. Auch dort hat man 
es ibrigens nach der Helmholtzschen Theorie mit nichtlinearen Gleichungs- 
systemen zu tun. 
Fir das allgemeine Glied der Entwicklung kann man nun ansetzen: 


n=1 
Pena = (OP By Agn—1,2y41 cos (2 w+1) ot 
w=0 


Jan—1 = Co a, Bont, 2441 8in(2 w+ 1) at 


pao 
one gies Anca cos{[(2u4+1)@+ 2,]t+ 6} 
+ Con-t2n41 CoS{ [(2 4 + 1) Q,]t— 6}, 
bon = Ay 3, Agno, Cos[|(2 uw + Q,) t+ 6] uo 
= 
+ Aono, cos[ (2.40 — Q,)¢— é], 
Jan = Ay 3, Bono, in| (2 uw + Q,) t+ d] 


de gh oy Sin[(2 wm —Q,)¢— 6], 


ue 


oll 
b, 
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Setzt man dies in die Rekursionsformeln (13) ein, so erhalt man 


Bann = S22 3 Cana, (GOREN E  cosaut iat — | 
: ee pg 005(2¢— 1) 4), 
fans = “exe By Con ay ce fei eres ie sin(2 +1) wt 
er) 
Bai ns a (-B Asniay + Bonso,) Eat to + Ot 9} 
+ (B Aan ay + Bon, 2,) set [eters = 
+ (A’, B’)-Glieder , 
Ce ara renee ae ie 
x cos{[(2 4 + 2) a+ Qz]¢+ aa tnets oe a a 
x cos[(2 uw + 2%) t+ 6] + (C’)-Glieder , 
Jan+2 = ae oy Con41, 2u41 = eee, = ae v 
x sin{[(2 w+ 2) @ + Q5]i+ os ae Seri : a 
x sin| (2 wo + Q,)1+ 6] + (C’)-Glieder , 
Aa — — ee Eten, ea vee FE cos2(ut+1) ot 


B Agnst,2 ur + Bo nst,0u+1 


gee cos (2 uw t) 3) 


a 


+ 


(Q% = Q,). Die (A’, B’)-Glieder haben denselben Aufbau wie die unmittelbar 
davorstehenden Ausdriicke (A, 2,/2) X...) und entstehen aus diesen, wenn 


1) Die Glieder mit noakeat + 6, und die Konstante in 72, ,9 sind noch in diesen Formeln ent- 


halten, werden aber durch die Integrationskonstanten zum Verschwinden gebracht. 
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man A;,, durch Aj, 1, B,,, durch B,,,, 25 durch —2Q,") und 6 durch —6 ersetzt. 
Ebenso erhalt man die (C’)-Glieder aus den jeweils davorstehenden Summen 
durch dieselben Einsetzungen und Cj, statt C;, 4. 

Zunachst ergibt sich mit Hilfe der Rekursionsformeln (15) durch vollstandige 
Induktion, dass sich die £,,, Ym und 7, (m = 1) alle in der Form (14) darstellen’ 
lassen. Dariiber hinaus lasst sich mit Hilfe von (15) die absolute Konvergenz 
der Reihen nachweisen. Definiert man zu diesem Zweck 


de = pe Lane ? Re, 3 Kero , 


und schliesst man schmale Intervalle | (2+ 1) o + 2,| < A um die Reso-. 
nanzstellen herum aus, so folgt aus (15) 


, Os Papas 
“H=0 


Rea’ 


p| Max(| |, 1) | 2p| Max(| |, 1) 
A A 


@) 
P, nL = | Ren ? Ove = 


BANCO INE na : 
Renan See tL (Pas + Cans) SME Ry, 


mit 


a, — 20] Max(|B) 1) 

VA 
Daraus ergibt sich nach dem Quotientenkriterium, dass Y Ry, «?” fiir «< My} 
konvergiert, und aus dem Majorantenkriterium folgt dann nach den obigen 
Abschatzungen auch die Konvergenz fiir Fane &7"+1 und 2'Q5n4, 02", und 
die Reihen fiir die stationaren Lésungen (A, = 0): 2’ ony, «?"+1 usw. konver- 
gieren in diesem Gebiet absolut und in jedem abgeschlossenen Teilgebiet gleich- 
massig. 

Fiir grosse Werte von |«|, etwa fiir |o| =|2,|/(1 — a4) mit 0< 4 <1, 
kann man offensichtlich M, durch M 1=M,VA(A w) ersetzen. Die Reihen 
konvergieren mit wachsendem || also immer besser. 

Beriicksichtigt man die Prazessionseigenschwingung, so muss man zusiitz- 
lich diejenigen Intervalle |m w+ Q,| = A ausschliessen, fiir die m gerade ist, 
da hier die Nenner der Summen- und Differenzschwingungen nicht nach unten 
beschrankt sind. Unter dieser Zusatzvoraussetzung folgt dann 


| 2,| Max (1, | B)) 
Rina = Mtg Mita Mee ae 
iui espa V4 Min (A, |o]) 


Daraus ergibt sich die Konvergenz der A,-Reihen und wegen M, <= M die 
absolute Konvergenz der Gesamtreihen fiir « << M-1. 


1) Die Glieder 22, bleiben stehen. 


a 
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_ Die Reihen konvergieren also, wenn « klein genug ist, mit Ausnahme 
kleiner Umgebungen der Resonanzstellen und des Nullpunktes. Durch geeignete 
‘kleine Wahl von « kann man in jedem vorgegebenen Punkt mit Ausnahme 
“von w = 0 und der Resonanzstellen selbst Konvergenz erzielen. Resonanzstellen 

sind dabei: 

1. |w| = |Q,| fir alle Schwingungen; 

2. |@| =|Q,|/(2% +1) fiir die Ober-, Summen- und Differenzschwin- 

-gungen von einer gewissen Ordnung ab; 

3. |a| = |2,|/(2 n) fiir die Differenz- und Summenschwingungen ab einer 
-gewissen Ordnung. 

Auch in der Nahe der Resonanzstellen kann natiirlich die Amplitude nicht 
beliebig gross werden. Deshalb miissen immer A, und C, so gewahlt werden, 
dass |9| = Nj, = const ist. Von den drei Integrationskonstanten Ao, Cy und 6 

‘sind also nur noch zwei aus den Anfangsbedingungen bestimmbar. 

__ In dem oben angegebenen Bereich der absoluten Konvergenz kann man die 

-Reihen auch umordnen und zum Beispiel die Glieder mit gleichen Argumenten 

der Winkelfunktionen zusammenfassen. Die stationdre Lésung wird dann durch 
normale Fourier-Reihen dargestellt. 


Summary 


The object of this paper is to give a theory of the experiments described in 
the foregoing paper. The applied system of non-linear differential equations by 
A. Féprx is first solved approximately for ‘small angles’, that means small 
deviations of the axis of the gyroscope from the vertical direction, the cases 
without and with consideration of damping force being discussed separately. All 
characteristic features of the motion observed in the experiments are explained. — 
After that, an exact particular solution is given for the case of a circular periodic 

force. An approximate calculation is carried out for the case of a periodic elliptical 
force acting on a fast rotating gyroscope. The solution contains an infinite number 
of resonance frequencies (subharmonic resonances and resonances with sums and 
differences of the forced and natural frequencies). 


(Eingegangen: 21. Mai 1954.) 
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Sur la thermodynamique des processus irréversibles 


Par KyriLtLE Poporr, Sofia) 
Quatriéme mémoire 


Dans une série de publications [1], [2], [3]?) nous avons considéré la thermo- 
dynamique des processus irréversibles dans un systéme adiabatiquement isolé et 
nous avons établi que, dans ce cas, les relations phénoménologiques d’ONSAGER 
sont des intégrales premiéres d’un systéme d’équations différentielles du type 
newtonien, ot l’expression de l’entropie, autour de létat d’équilibre thermo- 
dynamique, joue le réle d’une fonction potentielle des forces. 

Dans ce qui suit nous allons considérer le cas important ou le systéme 
échange de la chaleur avec le milieu ambiant, la température et la pression 
restant constantes et ot la fonction @ joue le réle d’une fonction potentielle des 
forces. Nous appliquons les résultats ainsi obtenus a la théorie des phases. Les 
développements mathématiques présentant une grande analogie avec les déve- 
loppements mathématiques de notre troisiéme mémoire cité, nous allons pour- 
suivre ces recherches et tirer quelques résultats nouveaux, concernant les pro- 
priétés des intégrales des équations différentielles définissant les processus dans 
les deux cas. Enfin nous traitons la question des barriéres énergétiques. 


Thermodynamique des processus irréversibles dans le cas ou la 
température et la pression restent constantes 


Considérons le cas ou le systéme échange de la chaleur avec le milieu 
ambiant, la température et la pression restant constantes. Ces processus sont 


caractérisés par la fonction a ‘ 
@=S-—TP* (1) 


ou S est l’entropie, U l’énergie totale interne, V le volume total, f la pression 
et 7’ la température absolue du systéme. Dans ce cas ® croit d’une maniére 
monotone et atteint son maximum a l'état de l’équilibre thermodynamique. 

Soient &,, &,...,&, les variables, autres que T et #, définissant l’état du 
systéme, et 9 £9... £9 leurs valeurs a l’état d’équilibre. Posons 


i 
A® = DE) — GE) = ee eH Sigs Mi Xe» AVEC p= k;, 4 =F, — & (2) 
vk 


1) Institut mathématique de l'Université de Sofia. 


2) Les chiffres entre crochets renvoient a la bibliographie, page 386. 


J 
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“ou 2 g;;, x; x, est une forme quadratique définie positive. En suivant la mé- 
_thode de notre troisiéme mémoire, sans apporter aucun changement, on 
-montre que, dans ce cas encore, les relations phénoménologiques 


: y : (0) —A® 
e = Ly Xy fee g tte te XS vee X,= “SS i Lay £59 3) 
sont des intégrales premiéres du systéme 
f a*%, 
: AP a ea, a, (4) 
_intégré das les conditions x,(+ oo) = 0. 
De (4) on tire immédiatement 
iL 12, / 12 
Se ag ee) AO = 0. (5) 
; Application 4 la théorie des phases, T et p restant constants 


Avant de poursuivre notre exposition nous voulons rappeler une propriété 
_ des fonctions homogénes. Soit f(x, v) une fonction homogéne de x et de y de 
- degré m. On a, comme on le sait, 


Peikee y) 
ax 


0 ) 
+y He) — mf, 9) 


On montre facilement que 0/(x, y)/0x, Of(x, y)/0y sont des fonctions 
homogénes de degré (m — 1) et que par conséquent les secondes dérivées par- 
- tielles de f(x,y) sont des fonctions homogénes de degré (m — 2) et ainsi de suite. 

Considérons maintenant un systéme de deux composants en deux phases, 
étudié par PLANCK [4] et soit 


M,=Mi+M!], M,=M,+ MM, (6) 


les masses M, et M, des composants, partagées parmi les phases (’) et ("). 
La fonction caractérisant |’état du systéme est ici 


Cope Phy (1 bis) 
T 
elon a S—=—o +S", U=U'+U",V=V'+V". La température T et la 
pression p étant constantes et les mémes pour les deux phases, l’on aura 
PD = Dp’ i @". 
S,U, V étant des fonctions homogénes de la masse du dégré un, on a le méme 
pour ®. De plus la composition relative de @’ étant la méme pour M,, Mj 


. 
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} 


et pour » Mj, n M3, ©’ est une fonction homogéne de degré wn par rapport 
a Mj et M3. Par conséquent ; 
Op’ OD’ 
om, “OM, 
sont de degré nul et 
0?’ 0?D’ G2Ds 
Mi?’ 0M dM,’ 0M2 


de degré moins un (— 1) par rapport a Mj, My. On a ainsi 


OD’ ‘ MM. ; OD’ ? ; OD” r Op" i 
om hel =y(c’) , OM, pr(c’) , ome ~ ¥i (c"), om ~ Ye ’) 
010 inl 02D” 1 (7) 


aT ie We" 
amo, ~ ua ? °) omrome — uz ? ©) 


ou l’on a mis 
hh maa M; Uke. om My 


CPM ee 


D’autre part, 0@'/0M; étant une fonction homogeéne de degré nul, on obtient 


020’ © NM Oe Mi 026’ 1 (ec) (8 
amie — yee 8 (8) Qe = oer ot omy = ee ee 
et de méme 
0°D” M; Wet 07D” an 1 We Al 
oMs? ~— MaaS yp (c"), oman = ar Ds (9) 


Formons aprés cela A4®. En posant 
@(M,, M,) = ©'(M,, M;) + ©" (M1, M2) 
et en tenant compte de ce que 
6M, =—dM;=m,, 6M,=—6M,; =m, 


on obtient d’abord, en choisissant pour M{, M3, Mj, Mg les valeurs de ces 
variables correspondant a l'état de l’équilibre thermodynamique, 
00’. 08" 44 88 aot 
OME GME a> NT mes ooh a 


=0. (10) 


Ces équations, qu’on peut écrire aussi, ayant égard a (7), sous la forme 


vile’) — py(c") = 0, yale’) — yh (c") = 0 (11) 


\ 


déterminent les concentrations c’, c” a l’état de l’équilibre thermodynamique. 
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a Ainsi, autour de away thermodynamique, l’on obtient 
Bio +1 ( ews “o,) 0M? + 2 ( 020’ ao" ] 
a Ze 0M? nice 0M; 0M; a oM}{ el 
f ' 020’ 020" 
d x 6M! 6M) + vee ae sr) 5M;:| x 
! 
4 
1 se gee Mz 
=F ilar Pe) + we e"(c")| mi 


Beer ag rn Ms \* Tap ( (LS Ms \* 
= 4 [- ioe) (Fe - Fe) — Ma oC (are — aie) |- 
(—A®@) étant une forme quadratique, positivement définie, on voit d'ici que 


— g'(c’) et p"(c") sont positifs. 
; Avec cette expression de A@ les équations (4) deviennent 


— dm, 0(—A®) = M, ee M3; pM y (c’) gp’ (c”) 
mae om, ec 9 (C!) + pre PC ) a woes Ley | see 
j d?m, 7! 0(—A®) gy’ (c’) yp” (c”) y’(c’) py” (c”) 

Sat ~~ om, ar + ae Jot Pare + uy J me. (13) 


_ (—A®@) étant une forme quadratique positivement définie, les racines 7?, 73 
de l’équation algébrique en 7? 


Ma mien, Ms he hat OE) 

mie PC) + gre JA eat sl re ) 0 (14) 
gc’) , (c") p(c’) , (c’) ie 

nen (= M; + ae ) = Mi + Mi r 


sont positives et par conséquent deux des racines de cette équation en 7 sont 
négatives et les deux autres positives. Soient 


foot, eh fe rae 
Les intégrales de (13) vérifiant les conditions 
limm,() =0, limm,(t)=0 pour ¢> + 0° 
sont de la forme 
Mm, = 0,6 + Cie, m,=C, fi gre Cape e*. (15) 
ou C, et C, sont définies par les conditions initiales 


f= 0 wh =C, + CEs mi = C, B, + Cz Be- 


382 KyYRILLE POPOFF ZAMP 
Pendant la durée du processus l’on aura, suivant (5), 


(S2)'+ (Ge) =-249. (16) 


L’élimination de C, e”’, C, e’ des équations (15) et de l’expression de 
dm'/dt d’un part et de dm,/dt d’autre part conduit aux rélations phénoménolo- 
giques avec’ L,, = L,,. 

Aprés tout ce que nous avons dit dans nos mémoires précédants il est 
superflu de donner les expressions explicites de m,, de m, et des relations 
phénoménologiques correspondantes, l’intégration du systéme (13) ne présen- 
tant aucune difficulté. 

Considérons le cas spécial ot la seconde phase ne contient que le composant 
My, c’est-a-dire ot l’on a Mj = 0, Mj = M,. Dans ce cas l’ona 


Om” 0?@” 02D” 


& he He 
pt Noe TUE ally ome —%> aM? omy = 9» omze —0- 
De 

OD” Wt OD” Ud Y 
omy “4 + oir Mz = ® 
on tire 


Op” 
YW Ui as = 
oC Mee ome = C = const. 


Les équations 


CO nee GOs See 0 
OM, 0Mz ee eg 


a 


déterminant 1’état de l’équilibre thermodynamique deviennent maintenant 


OD’ Om’ 


OC 2 “ome 


et l’équation (12), en négligeant les termes de degré > 2 en m, et m,, devient 


i> ane 
ea A Le 


s 


Ainsi le systéme (13) se réduit dans ce cas a l’équation unique 


d*m M any 
a Tee aa p(c’) my . (17) 


Intégrée dans les conditions de limite m,(+ co) = 0, cette équation donne 


-V(MiMi)e (cyt my oe 
mM 1 } ee a f 
1=Ce aa My, ue % fe". 
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Propriétés des intégrales des systemes du type (4) 


Nous avons vu que le systéme (4), caractérisant des processus irréversibles 

dans le cas ou T et p restent constantes est du méme type que le systéme 

_ d’équations différentielles correspondant au cas d’un systéme adiabatiquement 
isolé. Dans les deux cas les forces X; proviennent des formes quadratiques 
positivement définies, soit de 


1 
P= Z dy 8 Xi %y%, avec 8in= Exi- (18) 


Ainsi tous ce qui va suivre se rapporte aux deux cas et aux cas analogues se 
présentant dans l’étude des processus thermodynamiques irréversibles. 

Souvent on sera obligé de donner des poids différents aux forces X,. Ainsi, 
par exemple, dans |’étude de l’osmose le passage des molécules a travers une 
paroi poreuse dépend, entre autre, du parcours moyen des molécules, ce qui 
cause une sélection des masses. Ici encore nous nous bornons au cas de deux 
variables indépendantes pour mettre au claire Vinfluence du poids sur la 
forme des relations phénoménologiques. 

Considérons le systéme d’équations différentielles 


a*x% Ole 
ne ae Ed + S12 % = XL» 
2 (19) 
2 
ae =1 5 =A Bo X +A Bon Xp = AX, avec hee 
et les intégrales correspondant aux conditions de limites 
limz,(t)=0, limx,(?)=0 pour ?>-+oc. 
En posant 
gee xp = Bev" 
on est conduit, pour la détermination de 7 et de f, aux équations 
fr — 77+ Bb =9, A go + B (Abo — 77) = 9, (20) 
d’ou l’on tire pour rv? l’équation 
fat" S12 ie 
- | = rt — 7? (A Boo + B11) +A (811 822 — Sie) = 0, (21) 
A gor Ws I | 


ce qui donne 


| Ot ee 
ee rea Le > 


,4 
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Pour A = 0 la valeur de l’expression sous le signe du radical est positive. Cette 
expression ne s’annule pour aucune valeur réelle de 2 tant que gy. + 0. 
En effet de ; F 
(£2 5 Ss.) EU gie Oe (23) 


qu’on peut écrire sous la forme 
A? Bon — 2A (611 Bon — 2 G2e) + 4, = 0, 
-on tire 
= 2 eS Tae 
hp Ae 811 822 — 2 Ss Vee 811 §22 ; (24) 
Or comme gy, 47 + 2 B15 %1 % + Yoo % est une forme quadratique positive- 
ment définie, l’on a 


Fu Soa ia Of iq = 0, Son > 0. 


1° Dans le cas ott gy, = 0 les équations (19) se réduisent A 
ax ad*x 
ea bie a a 4 ban He (25) 
et les relations phénoménologiques correspondantes deviennent 
4% = 1, X,, wy = Ly, X, 
douse E,, = 0. 
2° Dans le cas ot g,. + O la valeur de expression sous le signe du radical 
de (24) sera négative; 1, et A, seront complexes et par conséquent les racines 


ri, v3 déterminées par (22) seront réelles et positives pour 2 > 0, comme on 
s’en rend compte en écrivant (22) sous la forme 


2 Siu + Ag 811 + A Boe \? 2 
h.2= 2 = +/( = 2 ad A (811 S02 — geo) - 


Ainsi les quatre racines de (21) sont réelles: 


HA 0; 0,0 emp, 44 =7q. 


Comme dans le processus physique on a lim x, = lim x, = 0 pour ¢ + + oo, 
il faut mettre C,; = C,=0 et prendre les intégrales du systéme (19) sous la 
forme 


% = Cre"? 4+ C, et! i= C, ye oe Car. ee | 
y= C, Bye" + Cy B, e%', te = Cir, Bye Cory bse". | 28) 

Avec ces valeurs de *, et %, on obtient pour les expressions de Ay etx, 
Ay = bin %1 + fia Hy = Ce (G12 + 12 81) + Cy e™* (ey + Bip Bs) , 7) 


Xe = Bis Hy + Soy %y = Cie (gia + Bon By) + Cye” (B15 + Boo Bs) « 
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L’élimination de C, e”’, C, e’*' des expressions de xj, X,, X, d’une part et des 
expressions de %9, X,, X, d’autre part conduit aux relations phénoménolo- 
- giques 


figes Lig Ky to Lig Ay, = 1, 2 


~ 
- 


_En suivant le chemin de notre premier et troisitme mémoire on trouve L,.=L,,. 
Dans le cas de racines doubles l’expression sous le signe de radical de (22) 
doit étre nulle et par conséquent l’on aura g,, = 0. Les équations (19) se réduisent 
aux équations (25) déja considérées. 


x 
4 


Les barriéres énergétiques envisagées au point de vue de la théorie des 
processus irréversibles 


La théorie des processus irréversibles exposée dans ce travail éclaire d’une 
maniére nouvelle la question des barriéres énergétiques. Pour ne pas compliquer 
_lexposition considérons le cas n = 2. Les équations (4) s’écrivent dans ce cas 


a*x, d* x5 bj 
ee Bie eae gy ST Baa Xen, AVEC Fi2 = Bin « (4 bis) 
L’intégration de ces équations par des expressions de la forme x, = « e"*, 
%_ = B e’* donne pour 
= 2 
e=82, B=—-u-—7). 


En désignant par 7, et 7, les deux racines négatives l’on a pour les intégrales 
s’annulant pour ¢ = + co 
Yet 


X= Cy Sie emt 4 Co S12 gr; Ke = —Cy (Bu — 77) ent — Cy (8 —7 5) ee”. 


En désignant par x?, x$ les valeurs de x,, ¥, pour ¢=0, on aura pour déter- 
miner les constantes d’intégration Cy, C, les équations 


w= Ci Sig Cokie, xd = —C, (81, — 77) — Cy (811 — 4) - 


Les valeurs de C, et C, ainsi déterminées donnent, en général, pour les dérivées 
dx,|dt, dx,/dt au moment t = 0 des valeurs différentes de zéro. Pour que la 
valeur de dx,/dt, par exemple, soit nulle au moment ¢ = 0 il faut qu’on ait de 


plus 
0 = Cy 1% S12 + C272 S12 - 


Par conséquent, pour que les trois derniéres équations soient compatibles il 


faut qu’on ait 
ay a §12 
x3 fil, a Ree 


ZAMP V1/25 
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On trouve de méme pour que le flux de x, au commencement soit nul la condi- 


tion 

ik ane a ~ 
a ae 12 doit étre satisfaite. 
5 812 


Dans le cas de » > 2 on pourrait avoir des barriéres correspondant aux racines 
multiples. 


BIBLIOGRAPHIE 


[1] Kyritte Poporr, Sur la thermodynamique des processus ivvéversibles, ZAMP 3, 

42-51, 440-448 (1952); 5, 67-83 (1954). ’ 

[2] KyRriLLe Poporr, L’osmose a la lumiére de la théorie des processus theymodyna- 
miques ivvéversibles, Ann. Phys. [12], 9, 261-268 (1954). 

[3] KyRILLeE Poporr, C.r. Acad. Sci. Paris 233, 648-649 (1952); 236, 785-786 
(1953); 236, 1640-1641 (1953); 237, 698-700 (1953); 238, 331-333 (1954). 

[4] Max PLancK, Vorlesungen tiber Thermodynamik, 9e éd. (Walter de Gruyter, 
Berlin et Leipzig 1930). 


Zusammenfassung 


In einer Reihe von Publikationen in der ZAMP und in denC.r. Pariser Akademie 
der Wissenschaften haben wir gezeigt, dass die phdnomenologischen Relationen 


Be EE ee mit Lj, = Ly; , (1) 


im Falle eines adiabatisch isolierten Systems erste Integrale eines Systems Diffe- 
rentialgleichungen von dem Typus 


ad*x,; 
“ap (2) | 
sind, wo a(— AS) 
xX; = Ox (3) 


und wo 4S den Ausdruck der Entropie im thermodynamischen Gleichgewichts- 
zustand bedeuten. In der vorliegenden Arbeit wird der Fall betrachtet, in dem das 
System, bei konstantem Druck und konstanter Temperatur gehalten, Warme 
mit der Umgebung austauscht. Es wird gezeigt, dass auch in diesem Falle die pha- 
nomenologischen Relationen (1) Integrale eines Systems der Differentialgleichun- 
gen (2) sind, in denen 


(4) 
und 
@= S$ 


darstellen. Die so entwickelte Theorie wird auf die Phasentheorie angewendet. 
Im weitern folgt eine Untersuchung des Systems 
aH, = ax, ; 
5 Ge 4s, Dems Ou, 
in welcher nachgewiesen wird, dass die ersten Integrale dieses Systems zu pha- 
nomenologischen Relationen (1) mit L,, = L., fiihren. 
Zuletzt wird die Frage nach den energetischen Barrieren erértert. 
(Regu: le 26 juillet 1954.) 
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: 

Bestimmung der Eigenwerte und Eigenvektoren einer Matrix 
mit Hilfe des Quotienten-Differenzen-Algorithmus 


Von HEINZ RUTISHAUSER, Ziirich!) 


Der Verfasser hat in zwei friiheren Arbeiten?) Theorie und einige Anwendun- 
gen des Quotienten-Differenzen-Algorithmus (QD-Algorithmus) behandelt. Die 
vorliegende Arbeit befasst sich mit einer Erweiterung des QD-Algorithmus, 
welche das im Titel genannte Problem zu lésen gestattet. Zuvor behandeln wir 
nochmals kurz die grundlegenden Eigenschaften des QD-Algorithmus: 


Wenn eine Potenzreihe /(z) = / s,,/z**1 durch ihre Koeffizienten gegeben ist, 


. 0 
so kann man, von den Quotienten gq!” =s,,,/s, ausgehend, mit Hilfe der 
Formeln 


erry) 
(v+1 (v+1 eee +1) 
Bt ae Ge (*) 


o 


ein zweidimensionales Schema von Zahlen q®, e” aufbauen: 


So 
qi 
S, e 
qi? Ge) 
Bike Al ead & 
o  @  @® 
er ar 


Die in der obersten Schragreihe dieses Schemas stehenden Zahlen q"’, e{”) 


treten auch in der Kettenbruchentwicklung der Funktion f(z) auf, es gilt nam- 
lich die formale Identitat: 


So | g | ef | | ey | 
lea y 5 e | « ; 2 Zz 


Eine der wichtigsten Anwendungen des QD-Algorithmus beruht darauf, 
dass man das QD-Schema der Funktion /(z) auch von oben nach unten auf- 


1) Institut fiir angewandte Mathematik der ETH. 

2) Der Quotienten-Differenzen-Algorithmus, ZAMP 5, 233-251 (1954), und A nwendungen des 
Quotienten-Differenzen-Algorithmus, ZAMP 5, 496-508 (1945). Diese beiden Arbeiten sind im fol- 
genden immer als I und II zitiert; ein Hinweis wie etwa (II, 7) weist auf die Formel (7) der Arbeit II. 
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bauen kann, wenn die Funktion f(z) bereits als S-Kettenbruch gegeben ist, das 
heisst, wenn man die oberste Schragzeile des QD-Schemas kennt. Die Rekur- 


sionsformeln hiefiir entstehen durch Umstellung der Formeln (*) (vgl. II, § 5): 
(v) ot) 
Io+1 °c 


i 
fr ado + ef — etn, etn = Bae 
co 


Auf diese Weise lassen sich leicht die Pole von f(z) berechnen, da diese nach 
I, § 7, die Grenzwerte der g-Kolonnen des QD-Schemas sind: 
Ap gS (or) 2, 2. 2¥. 


v—>co 


Die vorliegende Arbeit zeigt nun, dass dieses Prinzip auch auf Vektoren 
tibertragen werden kann, das heisst, man kann Rekursionsformeln angeben, 
durch die eine dem QD-Schema analoge Anordnung von Vektoren x) erzeugt 
wird, die fiir » + oo gegen die Eigenvektoren v, der Matrix A konvergieren. 
Diese Rekursionsformeln sind ein diskontinuierliches Analogon zu den Formeln 
(24) und (25) in, [3)4). 


1. Die Bestimmung der Eigenwerte 


Zur Bestimmung der Eigenwerte einer Matrix A (oder eines Integralopera- 
tors) kann man auf die in I, § 2, beschriebene Weise einige Schwarzsche Kon- 
stanten berechnen und dann das QD-Schema aufstellen. Die so erhaltenen 
Werte q’”, gf) geben dann einigen Aufschluss tiber die Lage der beiden ersten 
Eigenwerte, doch werden genauere Aussagen aus den in II, § 3, genannten 
Griinden nicht zu erhalten sein. Es empfiehlt sich vielmehr die Anwendung des 
QD-Algorithmus in der progressiven Form (II, § 5), insbesondere wenn man 
noch weitere Eigenwerte sucht. 

Die fiir diesen Zweck notwendigen Anfangswerte g), e) beschafft man sich 
im Falle eines Eigenwertproblems am besten mit Hilfe des Biorthogonalisie- 
rungsalgorithmus von C. Lanczos [1], [2], [3]; die Rechenvorschrift fiir dieses 
Verfahren lautet wie folgt: 

Rechenvorschrift A: Beginnend mit einem beliebigen Vektorpaar +,, y,?) 
berechne man rekursiv fiir o = 1, 2, nee ee 


i ; ; Ae 
Xo 44 A Xo Ay Xe Bean Ng-1) Vor. =A Vo — Ky ples (Pee Vesa 


ci 

) 

gnitte fig 0, Bleep NE ae @) 
(¥%¥o-1) Vo-1) (¥a, Vo) 

1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 401. 

*) Bei ungliicklicher Wahl von x, und y, kénnen allerdings die in [8], § 4, beschriebenen Falle 
eintreten. Wir machen aber hier und auch weiterhin die stillschweigende Voraussetzung, dass die 
Anfangsvektoren ftir den BO-Algorithmus immer in allgemeiner Lage seien. 

3) Fir o = n (Reihenzahl der Matrix A) miissen diese Formeln %n41 = Vn4i = 0 ergeben, was 
man als Kontrolle benititzen kann. Demgegeniiber sind die ebenfalls im letzten Schritt berechneten 
Grossen o%, und [Pa wesentlich fiir die weitere Rechnung. 


: 
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Damit hat man zwar erst die Jacobische Form der Matrix A, aber durch 


“Umkehrung der Formeln (I, 19) ergeben sich aus den «,, B, sogleich die ge- 
‘suchten Anfangswerte g, e®): 


0 Py 
= a,, q (0) ’ | 
, qy 
ft) Pe gl) o _ Po 
q Ap — Ey; i) (0)? 
2 
NS oer kta gore OY 2) 
ee 0 0 By-a 
Gi yi Oo ety = a 
Wn-1 
0 0 
q) ae Cee ae , 


‘Die weitere Rechnung verlauft nach II, § 5. 


2. Das Problem der Eigenvektorberechnung 


Wir gehen aus von der Aussage, dass die durch den BO-Algorithmus aus 

%1, ¥, erzeugten Vektorsysteme %1, %, +--+, %n und 1, Vo) «++, Yn «annahernd die 

Eigenvektoren der Matrizen A und A* sind, falls die Anfangsvektoren %1, V1 

geeignet gewahlt wurden » (vgl. [3], Satz 7, S. 54). Zwar sind die dort verwen- 

— deten Anfangsvektoren x(t) = e4"? x) und y(t) = e4*t2 y. fir die gegenwar- 

tige Untersuchung nicht zweckmissig, jedoch kann man in vollig analoger 
Weise den folgenden Satz beweisen: 


Satz 1: Falls die Eigenwerte von A verschiedene Absolutbetrage besitzen, gilt 
fiir das mit dem BO-Algorithmus aus den Anfangsvektoren 


20) = AM of, ye — (A*)M 9 (3) 


Gu) ae x(2e) vy (2H) 5 yt (2H): 


2 val 


erzeugte Biorthogonalsystem x 
Fiir > 0 streben die Vektoren x2, y2" (¢=1,2,...,m) gegen die Eigen- 


vektoren der Matrix A bzw. A*?). 


Die praktische Bedeutung dieses Satzes ist zundchst gering, da der BO- 
Algorithmus mit den Anfangsvektoren (3) bei grossem numerisch instabil ist. 
Die Aufrundungsfehler haben dann naimlich eine mit wachsendem o zuneh- 
mende Verfalschung der Vektoren x?" und y°") zur Folge, welche so stark sein 
kann, dass die Aussage von Satz 1 vollig illusorisch wird. 


1) Diese Konvergenzaussage bezieht sich nur auf die Richtung der Vektoren. 


f 
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“se 


Angesichts dieser Tatsache scheint es zweckmassig, den BO-Algontinag 
nicht mit einem nach (3) vorbehandelten Vektorpaar zu beginnen, sondern mit: 
einem Paar von Vektoren x", y), an welchen alle Eigenvektoren gleichmassig | 
beteiligt sind1), Allerdings sind dann die x, y (6 =1,2,...,n) zwar ein 
Biorthogonalsystem, jedoch auch nicht angendhert die Eigenvektoren der 
Matrix A bzw. A*. Wir werden aber in § 3 zeigen, dass man aus diesen x, of | 
die den Anfangsvektoren (3) entsprechenden Biorthogonalsysteme x“), you 
rekursiv berechnen kann, ganz ahnlich wie man die g"), e mit Hilfe der F or- 


meln (II, 4) rekursiv aus den g"), e) erhalt. Dabei ist diese rekursive Berech- 
nung der «?” und y?” durchaus numerisch stabil.’ 


3. Rekursive Berechnung der Vektoren x24), y2u) 


Wir ziehen nun auch die mit dem BO-Algorithmus verbundenen Polynome | 
o(z)) in den Kreis unserer Betrachtungen. Diese kénnen unter Verwendung 
der in (1) auftretenden Gréssen «, und f, wie folgt rekursiv definiert werden: 


Pal) =0, Peel, 


(4) 
Po+r(@) = — %41) Pal?) — Bo bas(2) (= 0,1,...,n—1). | 


Diese Polynome gestatten eine explizite Darstellung der Vektoren x,, y, 
(o=1,2,...,); es gilt namlich®) 


Apes P(A) Me Mg N= ~,(A*) V1: (5) : 


Im weiteren wollen wir auch die Polynome p?”) in Betracht ziehen, die in 
analoger Weise aus den Anfangsvektoren x2") = 4« x und ye") = (A*)# yO 
hervorgehen. Hieriiber gilt: 


Satz 2: Man bilde mit den Anfangsvektoren x), y) einerseits die Schwarz- 
schen Konstanten s,, = (A* x, yy) und mit diesen das QD-Schema nach L, Ss. 
sowie das P-Schema nach 1, § 4. Fiihrt man anderseits mit den Anfangsvektoren 
xP"), y 2") aus (3) den BO-Algorithmus durch und bildet analog zu (4) die Polynome 
p?”), so sind diese identisch mit den gleichnamigen Polynomen in der (2 Lu)-ten 
Schrégzeile des obenerwihnten P-Schemas. 


Bewets: Bildet man Schwarzsche Konstanten sy mit den Anfangsvektoren 
af"), y?#) (an Stelle von x, y{), so erhalt man wegen s* = Sona ein QOD- 


1) Bin Vektorpaar, dessen Komponenten Zufallswerte sind, wird dieser Bedingung in der Regel 
gentigen. Noch besser ware freilich eine Vorbehandlung durch Multiplikation mit A—! bzw. (A*)-1, 
die die Anteile der absolut kleinsten Eigenwerte verstarkt. Bei positiv definiten symmetrischen 
Matrizen kann statt dessen Vorbehandlung mit einem Verfahren von Karusu[4] empfohlen werden. 
Dieses beniitzt A nicht. 

2) Diese Polynome sind in [3], § 5, mit Poss (A) bezeichnet. 


3) Vel. etwa [5], Formel (9), S. 983. Der Vektor §; entspricht dabei unserem Lira 


ovens 
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"Schema und ein zugeordnetes P-Schema, welche offenbar auch durch Ab- 
“schneiden des QD- (bzw. P-) Schemas der s,, oberhalb der mit 2 «4 numerierten 
' Schragzeile entstehen: 


~ 


; So 
i 0 
ie qh 
“J 
1 ey 
1 
| qy” is 
| Sp : et) el) 
, (1 
4 q2) 
2 2 
of, 9) sts, 
2m) 
ap” 
2 
ay ra So 41 ay z ‘ 
2 2 
qe oe 
elu) 


Damit sind aber die p?”) des P-Schemas nach I, § 8, auch die Polynome von 
Lanczos fiir die Anfangsvektoren x{’" ), y2"), was zu beweisen war. 

Nach (5) lassen sich nun auch die in Satz 1 auftretenden Vektoren x"), yO“ 
durch die Polynome des P-Schemas ausdriicken: 
(24) _ p?")(A) ah) — Alt p?")(A) x0) ’ 


Xo+1 


2 g 2 (6) 
yee) — pen(a*) ye") = (Ar) py (A*) yy. 


rauf, dass es fiir die hier auftretenden 


Der Vorteil dieser Darstellung beruht da 
die wir zum Teil schon in I, § 4, her- 


Polynome ~°” Rekursionsformeln gibt, 
geleitet haben. Aus (I, 5) folgt namlich 


2 pee) = PMG) +90" boa) » 


o-1 


ferner aus (I, 6) 


ce ee 2 peat N(2)=—2 gabe pene Aia) 
Daraus durch Addition 
Ppee a) = pS ez) + ge) pez) — 2 eft W peural(z) , (7) 


_, nund fiir alle wz, jedoch unter Beachtung von e\” = 0 


und zwar fiiro = 1, 2,.. 
die Matrix A an Stelle der Variablen z ein und 


und ®), ,(z) = 0. Setzt man in (7) 


i 
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multipliziert dann das ganze mit dem Vektor x(?”), so folgt mit (6) 


2u+2 2 
Apo, A) = Age Ga a = Aa hg ae 


und weiter wegen A x@“) = x(?u+2); 


2 2u+1 2 +2) 
a ted ee Colne ae | 
firo= 1,2, ..., # und fir alle be (8) 
Dabei ist immer e?“*”) = 0 und x?“ = 0 zu setzen. | 


Diese Formel erméglicht offenbar die rekursive Berechnung der in Satz 1 ge- 
nannten Vektoren «?”), x2”, ..., x2” — und damit die Berechnung der Eigen- 
vektoren — aus den x), wenn man sie fiir jedes feste ys (u = 0, 1,.. .) nachein- 
ander fiir o = 1, 2, ..., » anwendet. Analoge Formeln gelten fiir die Vektoren y, 
doch werden diese im weiteren Verlauf der Rechnung nicht mehr gebraucht. 
Die in (8) auftretenden Skalare g™, e®) lassen sich natiirlich ebenfalls rekursiv 
aus den g), e!°) berechnen (vgl. II, § 5). Wir wollen jedoch auch die in den 
Formeln (II, 4) nicht auftretenden Gréssen «2, B°*) berechnen, welche zu den 
Vektoren x?”, y@“) gehéren. Diese hingen mit den q°"), e2 nach II, § 7, wie 
folgt zusammen: 


ones yi een 2u+1 2u+2 (2@+1) ,(2u+1 
gfe +2) gH +1) 4 o(2n+1), pet) _ g@u- ) @(2u+1) (9) 


wahrend sich umgekehrt g?”, e?“) aus «2, BE” nach (2) ergeben. Insgesamt 
erhalt man so zur Berechnung der Eigenwerte und Eigenvektoren einer Matrix 
A die Rechenvorschrift B (siehe Seite 393). Weil diese die Kenntnis der durch 
den BO-Algorithmus berechneten Gréssen %, Bg, X, Voraussetzt, muss vor- 
gadngig die Rechenvorschrift A durchgefiihrt werden. 

Auf Seite 394 zeigen wir als Beispiel die Bestimmung der Eigenwerte einer 
dreireihigen Matrix mit Hilfe der Rechenvorschrift B; das Beispiel erlaubt 
allerdings keinen auch fiir grosses 1 schliissigen Vergleich mit anderen Methoden. 
Die Matrix sei 


gp eas al | 
A= Wa "2. 3" 
L326 


Wir verfahren erst nach Rechenvorschrift A mit x7 = y, =(1, 0, 0) und dann 
nach Rechenvorschrift B. Die Tabelle auf Seite 394 zeigt fiir 4=0 die als 
Ergebnis von Rechenvorschrift A erhaltenen Ausgangswerte 1°, «°, ° und 
weiter die mit der Rechenvorschrift B erhaltenen Gréssen. Man erkennt fiir 
4 = 4 bereits die ungefahren Eigenwerte und Eigenvektoren, 

Das Beispiel zeigt ferner, dass die Vektoren x?") nur der Richtung nach 
konvergieren, wihrend ihre Langen zum Teil gegen 0, zum Teil gegen co streben. 
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Rechenvorschrift B+) © 


Nip ot 38 
b+ = be 

BS 24+), o> e2H) 
| ge) + eH) ee e2u+1) > gear) (cere oe a) 

(24) ,(2 4) 
| Pe By geet) a ed cg t= 1) i 

ele 

Vv } ; 
ES oe 

geutd) =e eu +1) > our?) Petal Soe.) 

2u+2 

ae geet) > p +2) 
Y 
VI 
| 2u+2 
| pS 40H); 0> 424 ) ag 
a 


2 2u+1) ,(244+2) (2u+2 
babe + gb aG0 — bt” ans? > af 


1) Das sogenannte «Ergibt-Zeichen» ass bedeutet, dass die rechts davon stehende Grdsse in der 
angegebenen Weise zu berechnen ist. 0 > y bedeutet, dass der Index yz den Wert 0 annehmen soll. 
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_ Man kénnte jedoch an Stelle der Vektoren x?” die Vektoren 


ae i) 


&? Pie 
wv g? )g®)... ea 2) 


-betrachten, welche im eigentlichen Sinne konvergieren, und die Arbeitsvor- 


schrift B auf die Berechnung der &?”) umstellen. Dies ist jedoch nur eine Frage 


_ der Rechentechnik und soll hier nicht behandelt werden. 


Wenn die Matrix A komplexe Eigenwerte besitzt, so haben diese paarweise 


i gleiche Absolutbetrage, so dass Satz 1 nicht gilt. Wenn aber /,,.; = Am 2 ist, 


das Polynom 


jedoch keine weiteren Eigenwerte cea Absolutbetrag besitzen, so kon- 
vergiert wenigstens die von x@”),, x”), aufgespannte Ebene fiir ju oo gegen 


die von den Eigenvektoren v,, ,; und v,,,2 = V,,,, aufgespannte Ebene, wahrend 


— [ot + oe A+ Eee? ae das 5 Boe, 


. gegen ein Polynom mit den Wurzeln 4,,.., und /,, ,. konvergiert. Die Eigenvek- 


toren v,,,, und v,,,. erhalt man dann als Grenzwert von 
2m) (2y) 
toes + [A sa oa Kin , 
je nachdem man /,,.,, oder /,,:» fiir A einsetzt. 


Betrachtung des Abbrechfehlers 


Wenn man den durch die Rechenvorschrift B gegebenen Prozess bei einem 
bestimmten yu abbricht, so sind die Vektoren «°” nur angenahert die Eigenvek- 


_ toren von A. Den Abbrechfehler kann man am einfachsten durch die Betrach- 


tung der Beziehungen 


Ag PP S728) a Pel OM) 4p BRM Be (10) 


“ abschatzen, wobei natiirlich f°” = 0 und x\?"| = 0 zu setzen ist. (10) beruht 


darauf, dass die x°"), Yen fiir jedes ein Biorthogonalsystem nach Lanczos 
_bilden und die «?”), B?# ee Elemente der zugehérigen Jacobischen Form von 
A sind. Wenn man die y?“) in der Rechnung mitfiihrt, kann man iibrigens die 
Biorthogonalitat der beiden eorioiocshtinis x?” und y@") als Rechenkontrolle 


bentitzen. 


4. Ein quadratisch konvergentes Verfahren zur Eigenvektorbestimmung 


Der QD- -Algorithmus in der bisher angegebenen Form (ausgenommen in 
pi. § 7) konvergiert linear, das heisst, die in I, Satz 3, behauptete Konvergenz 
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der q” >A, und der e”) > 0 (fiir » > 00) ist ungefahr so gut wie die Kon- | 
vergenz einer geometrischen Reihe; dasselbe gilt fiir die Vektoren x?) der 4 


Rechenvorschrift B. Wenn aber die Matrix A symmetrisch ist, so kann sie den 


in §3 beschriebenen Prozess zur Berechnung der Eigenwerte und Eigenvektoren — 


wesentlich beschleunigen. Das Prinzip dieser Konvergenzbeschleunigung ist 
dasselbe wie in II, § 7; es beruht auf einer Verschiebung des Koordinatennull- 
punkts der 4-Ebene. Fiir das Folgende nehmen wir an, die Eigenwerte seien 
einfach und der Grésse nach geordnet: A, > A, > 13> +--+: > An- 

Wir gehen aus von den Grossen x), «!°, B®, die man mit Hilfe des BO- 
Algorithmus von LANCZOs aus einem Anfangsvektor1) ermittelt hat. Aus die- 
sen wiirde man nach Rechenvorschrift B schrittweise «”, B®, x°), dann «'*, 
B“, x) usw. berechnen. 

Nun ist aber der BO-Algorithmus im wesentlichen invariant gegentiber 
einer Spektralverschiebung, das heisst, wenn man ihn mit demselben Anfangs- 
vektor x, aber mit A = A—tE an Stelle von A durchfiihrt, andern sich nur 
die «, ZU &, = «, — t; die x und # bleiben unverandert. Wenn man also nach der 
Berechnung von «,, 8,, %> (gemass Rechenvorschrift A) die Gréssen «, = a, — t 
(o=1,2,..., m) bildet und mit diesen in die Rechenvorschrift B (Feld II) ein- 
geht, verlauft der ganze Prozess wie wenn man mit der Matrix A —?£ an 
Stelle von A gerechnet hatte: An Stelle der B?”, die ungefaéhr wie (A, .1/A,)"“ 


gegen 0 konvergieren, erhalten wir jetzt andere Werte Bev, die wie 
[(Asu1 — #)/(A. — #)]*“ gegen 0 streben. 


Fiir ¢ setzt man vorteilhaft einen Naherungswert fiir den kleinsten Eigen- 


wert /,, ein, dann wird (A,, — #)/(A,1 — #) klein, so dass pe ) rasch gegen 0 kon- 
vergiert (und ebenso «?”) gegen A, — t, sowie *?”) gegen v,). Damit wird die 
Konvergenz des Verfahrens zwar besser, bleibt aber linear. Es lasst sich jedoch 
bereits nach einem Schritt auf Grund der neuen «- und f-Werte ein besserer 
Naherungswert fiir J,, angeben, was man fiir eine weitere Spektralverschiebung 
verwenden kann. Auf diese Weise wird die Konvergenz immer schneller, bis 


schliesslich B?"} so klein ist, dass *?”) praktisch ein Eigenvektor ist und #2”), 
x2, ..., #2" einen beziiglich A invarianten Teilraum aufspannen. Man kann 
in diesem Fall «,, 8,1, %» einfach weglassen und dasselbe Rechenverfahren auf 
die verbleibenden Gréssen anwenden. Dies bedeutet, dass man nach Abspaltung 
des Eigenvektors v, im Teilraum der verbleibenden Eigenvektoren weiter- 
rechnet und so nacheinander alle Eigenvektoren erhalt. 

Die Zusammenfassung dieser Uberlegungen ergibt eine Rechenvorschrift C, 
die aus den Feldern I, II, IV, V, VI der Rechenvorschrift B mit den auf Seite 398 


angegebenen Erginzungen besteht. Diese Rechenvorschrift liefert nacheinander | 


1) Wenn man bei symmetrischer Matrix A ein Paar von gleichen Anfangsvektoren x, = yy 
wahlt, so liefert der BO-Algorithmus %¢ = yo fiir alle o, man braucht also nur mit einem Vektor- 
system zu arbeiten, und dieses ist orthogonal. 
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die Eigenwerte A; und die zugeh6rigen Eigenvektoren v; (j = n,m —1,..., 2, 1) 
(siehe Feld IX der Arbeitsvorschrift C). 
Praktische Uberlegungen haben gezeigt, dass man als Koordinatennullpunkt 


zweckmassig eine untere Schranke fiir die noch verbleibenden Eigenwerte von 


A wahlt: Ist der Koordinatennullpunkt nach « Schritten insgesamt um ¢ ver- 
schoben worden, so verschiebt man ihn im nachsten Schritt um 6,, wo 6, eine 
méglichst gute untere Schranke fiir die Eigenwerte der Jacobischen Matrix 


(24) 2 
x Bi = 
(2 2 
AN ple ee oh BS AS 
pee = 1. baa 2s. (11) 
aa Pra ae 
IS ny (24) (24) 
Se n-1 n—-1 
RS 2 
i! ol) 


ist. Feld IIIa der Rechenvorschrift C beniitzt eine Formel von GERSCHGORIN [6] 
fiir die untere Schranke der Eigenwerte von (11), doch gibt es noch andere und 
zum Teil wohl bessere Méglichkeiten. 


Nachtrag bei der Korrektur. 

Numerische Experimente lassen es als vorteilhaft erscheinen, zuerst nach 
den in § 1 erwahnten Methoden die Eigenwerte allein zu bestimmen und dann 
erst die Rechenvorschrift C anzuwenden, wobei die 6, so zu wahlen sind, dass 
t, nacheinander die n Eigenwerte durchlauft. Auf diese Weise erhalt man die 7 
Eigenvektoren in 7 Schritten. 


5. Eigenwerte und Eigenvektoren unendlicher symmetrischer Matrizen 


Fiir unendliche Matrizen bzw. Integraloperatoren kann man natitirlich den 
BO-Algorithmus von Lanczos?) nicht zu Ende fihren, wie dies fiir die Durch- 
fiihrung der Rechenvorschrift B oder C notwendig ware. Praktisch gilt dasselbe 
auch schon fiir endliche Matrizen grosser Reihenzahl. 

Bekanntlich kann man aber bei symmetrischen und positiv definiten Ma- 
trizen (bzw. Operatoren) auch dann eine vorlaufige Information tiber die gréss- 
ten Eigenwerte und die zugeh6rigen Eigenvektoren erhalten, wenn man nur ein 
unvollstindiges Lanczossches Orthogonalsystem und die zugehérigen «, B be- 
rechnet. Man kommt jedoch in die bereits in § 2 genannten Schwierigkeiten, 
sobald man versucht, genauere Resultate zu erzielen, indem man den Anfangs- 
vektor x, durch Multiplikation mit A vorbehandelt (sogenannte « Voriteration »). 


1) Vgl. Fussnote 1, Seite 396. 
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Rechenvorschrift C 


ae 
n—-1>n |- 
von Feld II tr, 
ok gat Ga B) a Dilys | 
|}-<— (Oe | 
[o>] 
F 
Illa Fe 
mtn a Ve Ve po ty 


(mit pe #) os #) 0) 


m | Meek 
Dork On > buss | e+ide | 
a2) — 5 Se!) | 
pew) on HN 
ge) (eo == 1. n — 1) 
ot(2 4) Pa ee) =6, = g’ ) mal Von Feld VI 


(Rechenvorschrift B) 


nach Feld IV 
(Rechenvorschrift B) 


Es diirfte daher von Interesse sein, dass die in § 3 entwickelten Methoden auch 
auf diesen Fall mit Erfolg angewendet werden kénnen. 

Wir setzen im folgenden eine symmetrische und positiv definite Matrix A 
voraus und bestimmen fiir diese ein unvollstandiges Lanczossches Orthogonal- 
system %,...,%», und die zugehérigen Skalare aj, a, ..., Om; Bi ++», Bm—1) 
indem wir mit einem geeigneten x,1) die Rechenvorschrift A bis inklusive o = m 
durchrechnen. Hier ist m eine nach Bedarf festgesetzte ganze Zahl <n; das 
zu entwickelnde Verfahren wird dann die m gréssten Eigenwerte und die zu- 
gehérigen Eigenvektoren liefern. 

Die Gréssen x, Gy tys3 By; ins By AG Xa. s;, SB», SiNd Zusaramen-mit dem 
durch Formel (1) fiir o = m ebenfalls berechneten Vektor x,,,, die Ausgangs- 


1) Fiir die Wahl dieses Anfangsvektors siehe Fussnote 1, Seite 390. 
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werte fiir das weitere Vorgehen und werden daher kiinftig mit dem oberen 


Index 0 versehen. Man erhalt aus den «°) und f° ohne weiteres ¢, ..., q0’; 


—e, ..., e) , (Feld III der Rechenvorschrift B), doch kann man daraus nach 


Vtoren x?),..., x! 


‘den Formeln in Feld IV nur noch gq, ..., g)_,; ef’, ..., ef, berechnen, weil das 


zur Berechnung von q'!) notwendige e fehlt. Ebenso erhalt man aus den q 


{1) i 2). 
e nach Feld V der Rechenvorschrift B nur noch «?),..., 0) 4; BP, ..., Bd 


Dagegen erlauben die in Feld VI stehenden Formeln die Berechnung der Vek- 
®) weil ja x),, am Anfang vorausgesetzt war. 

Um nun einen dem Anfangszustand analogen Zustand herzustellen, miissen 
wir also noch «®), B®), und x«@),, berechnen. Weil die x) ebenfalls noch ein 
Lanczossches Orthogonalsystem!) fiir den Anfangsvektor x bilden, kann 
man die fehlenden Gréssen mit Hilfe der Formeln (1) fiir o = m berechnen. Ins- 
gesamt ergibt sich die Arbeitsvorschrift D (siehe Seite 400). 

Die so berechneten Gréssen «?”, ..., a") konvergieren fiir 4 > oo gegen 


die m gréssten Eigenwerte der Matrix A und die Vektoren x"), ..., x},/ gegen 


die zugehérigen Eigenvektoren. Im Gegensatz zur Rechenvorschrift B verlangt 
Rechenvorschrift D bei jedem Iterationsschritt die Multiplikation einer Matrix 
mit einem Vektor. Dies ist offenbar der Preis dafiir, dass mit einem unvollstan- 
digen Lanczosschen Orthogonalsystem gearbeitet werden kann. Dennoch sticht 


das Verfahren vorteilhaft ab, wenn man bedenkt, dass fiir einen Iterations- 


schritt, der m Vektoren x?” simultan verbessert, nur einmal eine Matrix mit 
einem Vektor multipliziert werden muss. 


Anmerkungen zur Rechenvorschrift D 


a) Es mag rechentechnisch von Vorteil sein, den Vektor x“), (siche Feld 
VIII) gar nicht zu berechnen, sondern den Ausdruck fiir x2”) in die Formel (8) 


- (« = m) einzusetzen; es resultiert dann nach einigen Umformungen die Bezie- 


hung 

Buse) _ ov) oHs9) g2Hs% » (12) 
Diese Formel kann dann fiir o = m an Stelle von (8) zur Berechnung von gure 
verwendet werden (anschliessend an Feld VI der Arbeitsvorschrift D). 

b) Wahrend die Formeln in Feld VI der Arbeitsvorschrift D nicht sehr 
empfindlich auf Rundungsfehler sind, ist die Berechnung von x‘), (Feld VIII) 
in dieser Hinsicht weit mehr gefahrdet, vor allem fiir grosses uw. Es ist daher 
ratsam, den nach den Formeln in Feld VIII berechneten Vektor x24) einer 
« Nachorthogonalisierung» zu unterwerfen, indem man beniitzt, dass die x4) 
fiir jedes ww ein Orthogonalsystem bilden miissen. Diese Orthogonalitatseigen- 
schaft kann natiirlich auch fiir Rechenkontrollen bentitzt werden. 


1) Siehe Fussnote 1, Seite 396, 
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Rechenvorschrift D 


de => oc?) Gee: 102) ae, 900) Vil 


A #2u) berechnen 


| 
| 
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| 
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| 


y (A x 2H), wer) al 
lege OS ULW RE aT Wires ot 
II (4) es) S sft 
Oe ie (79. aoe) |} 
ies) “5, 20) | A #2H) = oe(2H) #2u) oe pew) xe > Oe 
qe") eC re ee 8 
ar ee > ae | 
IV | 
0 . > een | 
qe) + e@ a) « gone) Est gare) 
Cees easy dpa Bots tte) 
gy 7 | 
Y 
y. | 
gbnrd) gener eee) ged 8 ae 1) os 
Ge ee tL OME giae yg, tc te he 2) Pine As 
A 
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a2) be gq) #2u) a een +1) a es g 2+?) (o ed 
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c) Auch bei dem durch Rechenvorschrift D definierten Verfahren ist noch 
eine Konvergenzverbesserung méglich, indem man bei jedem Iterationsschritt 
eine geeignete Spektralverschiebung um 6, vornimmt. Sind die Eigenwerte 2, 


' der Matrix A der Grésse nach numeriert und ist /,, der kleinste noch interes- 


sierende Eigenwert, so sind die 6, so zu wahlen, dass die sich daraus nach 


_ Arbeitsvorschrift C, Feld IIIb, ergebenden ¢, das Intervall 0 < 4 <A,, még- 


lichst dicht ausfiillen. Diese Methode der Konvergenzbeschleunigung ist in der 
Relaxationsrechnung als Uberrelaxation bekannt. 
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Summary 


Two previous papers (in Vol. V) describe theory and some applications of the 
quotient-difference (= QD-)algorithm. Here we give an extension which allows 
the determination of the eigenvectors of a matrix. Let #0), ney. #0) be a coor- 
dinate system in which A has Jacobi form (such a system may be constructed 
with methods of C. Lanczos or W. GIvENs). Then the QD-algorithm allows the 


construction of a sequence of coordinate systems gen, ae ge (iO pa aire ctn) 
which converge for «4 > co to the system of the eigenvectors of A. 


(Ringegangen: 19. Juli 1954). 
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The Indentation of an Orthotropic Half Plane 
By Harry DonaLp Conway, Ithaca, New York, U.S.A.?) 


Introduction 


The first theoretical investigation of an indented half plane appears to have 
been made by SapDowsky[1]?), who calculated the pressure distribution under a 
flat-ended rigid block in contact with an elastic, isotropic half plane, the block 
being subjected to an axial force. Sapowsky ignored the shearing forces between 
the block and the half plane, and thus assumed that the half plane was deformed 
only by normally applied forces. 

Further progress was made by OKuBo [2], who treated the same problem but 
took into account the frictional forces between the block and half plane. OKUBO 
observed that, in so far as the calculation of the contact pressure is concerned, 
the effect of friction could be neglected. Due to the slipping which may occur in 
practise, it is probable that the frictional forces are smaller than those computed 
by OxuBo. For the above reasons and because the inclusion of the frictional 
forces considerably increases the complexity of the problem, further investi- 
gators have chosen to ignore them. 

A very considerable advancement was next made by SCHUBERT [3], who gave 
a simple general solution of the integral equation for the contact pressure. Using 
this equation, he calculated the distribution of pressure under a flat-ended block 
subjected to non-axial force, and also the pressures under a variety of shapes of 
axially loaded rigid blocks. A number of these problems have been discussed 
independently by SNEDDON [4] using a Fourier transform method, and by SEN 
[5], [6] using a complex function method. However, SCHUBERT’s method of solution 
is the simplest. 

Recently, OkuBo [7] has re-examined his earlier work, and has extended its 
results by assuming that the block is flexible but ignoring the frictional forces. 
He has given graphs of the contact pressures and deformation shapes for several 
ratios of the moduli of the block and half plane. 

The present work is concerned with the contact pressure between a rigid block 
and an orthotropic half plane, some physically interesting results being obtained. 


Analysis 


Consider an orthotropic half plane whose plane boundary corresponds with 
the *-coordinate axis, the plane being subjected to a normal concentrated load P 
acting at the origin. Assuming that the x-, y-coordinate axes correspond with the 


1) Cornell University, College of Engineering, Dept. of Mechanics. 
2) Numbers in brackets refer to References, page 404. 
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‘principal axes of orthotropy, it can be shown that the stresses [8] are given by 


se lg P (a4, + %) #7 y ’ 
oS (oh a + V5) (ah a + 9)’ 
a P (a, + %) y? ’ 
ue (Eat VA) (bE VF 
= P (a + %) ¥ y? 
oe a (ey A eV ) a ee *) 


Tt. 


’ 


where a2 a3 = S4,/Sg. and af + 03 = (Sop + 2 Sy2)/So2, the elastic constants being 
_ defined by the stress-strain equations 


ie = S31 % + Sig %y, by = Soo Oy + S125, Yay = Se6 Tay « (2) 


Substituting the first two of equations (1) in the second of equations (2), integra- 
' ting, and assuming that the displacement v = 0at¥= +4, ¥= 0, it is found that 


. f af v2 + y* 2 Ce Saal 
vy 2 m (a, — ay) [ (oi Soo — Sip) ae — (03 See — Sis) log | : 


(3) 
It follows that 


Oe ie bed P (a + M2) Soe 
Ox\y=-0 | Ve, o(*) 7X (4) 
Thus an incremental load p(u) du at wu from the origin gives 
(a, + %) S p(u) du 
De (5) 


and the slope ot the half plane under a block of length 24 is 


A. 
1 : a! p(u) du 


a, os 
2 (a, + a) Soo a, 0(*) ase) *—U 


(6) 


For isotropy, #, = % = 1 and writing Sj, = 1/E in the usual notation, we have 


A 
E aS pi Pu) du (7) 


20 x—U 
=A 


Using a solution due to. HAMEL [9], ScHUBERT [3] has shown that equation (7) 
gives the distribution of contact pressure in the form 


A 
1 Ef vb V4A2— 2 
plu) ep at 4 vol) V dx\ (8) 
nVA?— u* ‘2 at 


It follows that the results for various shapes of blocks found from equation (8) 
will also apply to an orthotropic half plane if E is replaced by the effective mo- 


dulus 
po 


(a + Ot) a 


= 
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i 
The block of greatest interest is that having a flat end, and this will be considered 
in detail. : : 
For axial loading, vj, 9(¥) in the integral in equation (8) is zero and the well 
known result follows = : 
u) = rn (9) 
10) = aes 

If the load is applied with an eccentricity b and contact is made over the entire 
length of the block, uv}, 9(¥) = B (a constant), the slope of the contact surface, 

the orthotropic version of equation (8) becomes 


1 EA 4 
: 10 
1) = [P+ 3 2.8] (10) 
For equilibrium 2 
A 
ay % 5 =) 4Pb 
Po= | plu) udu= 4B, A B fey mE AA 
=A 
Thus 
124 [ qt 
| ae | ieee (11) 
Plu) a VA — u2 a 


If b> A/2, it follows that contact will cease over part of the block. Assuming 
that the contact length is 2a, where a< A, it follows from equilibrium that 


P[b— (4 — ay) = | plu) udu = 7 Esa? p 


and 
o Po Sp aes aresy | 
Plu) aH ea LE me i (12) 
Now p(— a) = 0,.a = 2(A — b) and, from equation (12) 
Sees P 2(A—db) tu 
naa wy Vege (13) 


All the above results, repeated from the work of SCHUBERT (3), are independent 
of the elastic constants of the half plane. Thus it follows that the pressure 
distribution under a flat-ended block subjected to any loading is independent of 
the elastic constants of the half plane, including the special case where Oy = Op 
and the half plane is isotropic. This statement is only true for flat-ended blocks, 
and is not true for blocks having other contact surfaces, although the form of 
the pressure distribution will be the same. 
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Zusammenfassung 


Der Verfasser findet eine allgemeine Lésung fiir die Verteilung des Druckes 
zwischen einem Block und einer orthotropen Halbebene. Es wird gezeigt, dass 
die Verteilung des Druckes fiir den flachen Block mit allgemeiner Belastung un- 
abhingig ist von den elastischen Eigenschaften: der Halbebene und auch genau 

_ dieselbe, als ob die Halbebene isotrop ware. 


(Received: May 7, 1955.) 


On the Oscillating Cylinder Viscometer 


By Att A. K. IBRAHIM and ABDEL Monem I. KABIEL, Alexandria, Egypt?) 


In our paper the theory of an oscillating cylinder viscometer part IED), sla: 
relationship [Ref. 2, equ. (16)} between the magnification factor 09/®, the viscosity 
_» of a Newtonian liquid and the constants of the apparatus was given. This rela- 
tionship, although it covers liquids of any density and viscosity, is not a direct 
relation between 7 and the other constants of the apparatus. A disadvantage 
which has been already found in previous derived equations. 

In the same paper’), the previous relationship was approximated and a simple 
formula for the viscosity 7 [Ref. 2, equ. (19)] was deduced. Unfortunately, this 
formula is applicable only to Newtonian liquids with 7 > 1. 

The object of the present note is to make the first relation [Ref. 2, equ. (16) 
direct and applicable to Newtonian liquids with all values of 7. 

Introduce a dimensionless parameter defined by 


V2BR= = ((2*— 1j2—1}, (1) 


where ¥ may take any value, say 1, 2, 3, 4,.-- 


1) Physics Department, Faculty of Science, University of Alexandria. 
2) The theory of an oscillating cylinder viscometer, Part II, Aut A. K. Ipranim and ABDEL 


Monem I. KaBIEL, Z. angew. Math. Phys. 5, 398 (1954). 


; a, 
* : a r 
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Thus, it is found that 
: J ; 
(te) = pty = eon Ch DI (ap 


®, Ci + Di + 43+ BZ+ Cj + Dz’ 
where 
' 9 
ores 
9 
ie me 
9 Ps Sih 
Ab= {72 ot — 5s (x — 1) | sin 34 yo, 
= eee —) 3..6 
i 3 \ ie Vx 1 V4 o Vx | sin 54 be de 
: 3d 
(OF — yi? w2 + st x } sinh*—,- Vx, 
; 5} 3 — v <a WA 6d 
Di= oe ag Vt He #— 1 ¢ sinh fx, 
and 


9 (Zw? — 1) V(2 ¥— 1)? -1 
8aL R2 9 a? w3 : 


ee 


Equation (2) satisfies the experimental curve only in one point of a certain 
value of (m, and y,). 
Hence, from (1) we have, 


2 
eile ae ie (3) 
Ole # st) hn 1 


from which 7 is calculated. 


O$§ 


1 


20 


7] 
n(Y%) 
Figure 1 
The variation of[@9/@ (n) for paraffin oil. 
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Experimental verification : 

Using the dimensions of OLDROYD’S, STRAWBRIDGE’S, and Toms’ apparatus’), 
and taking * = 10, say, equation (2) becomes 

. ( Ln y ee 10-93 n# 

®,)  -4~ 41-42 nt — 1834 n® + 27620 © 


(4) 


Hence, on plotting the curve represented by this equation on the experimental 
curve obtained by OLDROYD, STRAWBRIDGE, and Toms (see fig. 1), the point 
_ of intersection corresponds to a value , = 3st. 

Hence, substituting in (3) we have 


n= 1-74 gem "st. 


This is in good agreement with what was found in tables. 


Zusammenfassung 


Bei Einfiihrung eines dimensionslosen Parameters ist es méglich, die Bezie- 
hung [1, Gleichung (16)] direkt herzuleiten und auf Newtonsche Fliissigkeiten mit 
allen Werten von 7 anzuwenden. 


(Received: February 19, 1955). 


Mechanismus der Vereisung unterkiihlter Wassertropfen 
durch disperse Keimsubstanzen 


Von Hans RupoLF PRUPPACHER und RAYMUND SANGER, Ziirich®) 


I. Experimentelle Untersuchungen tiber die Eiskeimfahigkeit der Stoffe 


A. Einleitung 


Die Niederschlagsbildung stellt ein sehr wichtiges Problem der physikalischen 
Meteorologie dar. Die Zahl der bis jetzt angestellten experimentellen und theore- 
tischen Untersuchungen ist sehr mannigfaltig, zeigt aber, dass man noch weit 
davon entfernt ist, ein einheitliches Bild tiber den Mechanismus der Phasentiber- 
gange von Wasser geben zu konnen. Vor allem sind die Vorgange, welche im Ver- 
eisungsprozess eine Rolle spielen, noch sehr unabgeklart. Im Vordergrund steht 
dabei die Frage der Fiskeimfahigkeit der Stoffe, das heisst ihre Eigenschaft, die 
Vereisung unterkiihlter Wassertropfen einleiten zu kénnen*). Um vermehrten 
Finblick in die Vorgange zu gewinnen, sind die nachstehend beschriebenen Unter- 
suchungen angestellt worden, welche im Sommer 1954 begonnen und vor kurzem 


abgeschlossen wurden. 


1) OLDROYD, STRAWBRIDGE, and Toms, Proc. Phys. Soc. (London) [B] 64, 44 (1951). 

2) Physikalisches Institut der ETH. 

3) Mit dem Ausdruck « Eiskeimfahigkeit » soll angedeutet werden, dass die Frage, ob es sich um 
Gefrier- oder Sublimationskerne handelt, in dieser Arbeit im wesentlichen noch offengelassen wird. 
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B. Durchftihrung der Versuche 


Die Versuche erfolgten in einer Kiihltruhe der Marke «Hill» nach dem Ver- — 


fahren von SCHAEFER-VONNEGUT?). Der Innenraum der Truhe bestand aus vier 
Kiihlkammern, die nur durch eine schmale Offnung am Boden der Truhe mit- 


einander in Verbindung standen. In der Kammer, in welcher die Experimente © 


durchgefiihrt wurden (Volumen etwa 80 1), war ein Holzgesteli eingepasst, das 
zum Teil mit schwarzem Chiffon und zum Teil mit schwarzer Gaze tiberzogen war. 
In der vorderen Seite der Kiihltruhe befand sich ein dreifach verglastes Beobach- 
tungsfenster. Der obere Teil der Truhe war zur Halfte durch das Kiihlsystem und 
zur Halfte von einer eingepassten Plexiglasplatte abgedeckt, in welcher sich iiber 
der Untersuchungskammer ein Loch zur Einfiihrung der Feuchte und der Impf- 
substanz befand. Pro Tag wurden nur Untersuchungen bei einer bestimmten 


Thermostateinstellung durchgefiihrt. Durch fiinfstiindige Vorkiihlung bei dieser ' 


Einstellung gelang es, bis zu Versuchsbeginn den Temperaturunterschied zwischen 
Boden und Plexiglasdecke der Kiihltruhe auf 10° C zu reduzieren. 

Zur Beobachtung wurde eine Projektionslampe (150 W) verwendet, deren 
Licht mit Hilfe eines beweglichen Spiegels nur im Moment der Beobachtung 
in die Kammer gelenkt wurde. Damit konnte die 6rtliche Erwarmung durch die 
Strahlung auf ein Minimum gesenkt werden. 

Die feine Zerteilung der Impfsubstanzen erfolgte nicht durch Verdampfen, 
sondern durch Luftzerstaubung der fein pulverisierten Substanzen. So war es 
mdoglich, die Zersetzungserscheinungen infolge Warme und Luftsauerstoff zu um- 
gehen und die Kristallstruktur der Substanz vor der Einstaubung eindeutig fest- 
zulegen. Die Substanzen wurden in einem Achatmérser pulverisiert (Teilchen- 
durchmesser um 0,5—-2-10-4cm) und durch einen 10-cm3-Puderzerstauber aus 
Plastik zerstaubt. Substanzen mit blattrigem, weichem oder biegsamem Charak- 
ter wurden zusammen mit Glas, das bis — 18°C nicht keimfahig ist, pulverisiert 
und so in sehr feine Form gebracht. Hygroskopische Substanzen wurden durch 
fortgesetztes Pulverisierungs- und Trocknungsverfahren in feinzerteilte. Form 
gebracht. Nur im Falle von KI, Nal, AgNO, und I, musste das Rauchverfahren 
angewendet werden. Zur Reinigung der Luft in der Kiihltruhe nach den Impf- 
versuchen verwendete man ein Elektrofilter der Marke «Trion». 

Fir die Temperaturmessung sind 5 verschiedene Thermometerarten (Queck- 
silberthermometer, Widerstandsplattenthermometer, Bimetaltorsionsthermo- 
meter, Thermistoren und Thermoelemente) nach 3 Gesichtspunkten gepriift 
worden: Ein gutes und rasches Anpassungsvermégen des temperaturempfind- 
lichen Teiles an eine gegebene Aussentemperatur, eine grosse Empfindlichkeit 
des Anzeigegerates und eine moéglichst geringe Platzbeanspruchung des tempe- 
raturempfindlichen Teiles. Das erste Kriterium ist an Hand von Temperatur- 
angleichungsdiagrammen untersucht worden, die fiir die einzelnen Thermometer 
aufgenommen wurden. Als weitaus am giinstigsten erwiesen sich die Thermo- 
elemente (Kupfer/Konstantan, 0,1 mm Durchmesser) und die Thermistoren 
Philips 83900-1). Zur Gewahrleistung einer méglichst guten Kontrolle der Tem- 
peratur wurden beide Thermometerarten mit Hilfe verschiedener Instrumente 
geeicht und je sechs links und rechts in vertikaler Richtung in der Kammer 
montiert. Um Turbulenzen in der Kammer auf ein Minimum zu reduzieren, sind 
auch bei den Thermistoren sehr diinne Zuleitungsdrahte verwendet worden. 


1) V. J. ScHarrer, Chem. Rev. 44, 291 (1949); B. Vonnecut, Chem. Rev. 44, 277 (1949). 
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Mit dieser Temperaturmessanlage wurde der Temperaturverlauf in der unge- 
-stérten Kiihltruhe nach Einschalten der Beobachtungslampe und nach der Er- 
zeugung des Nebels durch Hineinhauchen kontrolliert. Wahrenddem der Einfluss 
der Lampe nur wenige Zehntelsgrad betrug, zeigte sich nach einem fiinf- bis 
siebenmaligen Einhauchen ein Temperaturanstieg an allen Messpunkten von etwa 
7 bis 12°C. Das Absinken der Temperatur erfolgte anfangs sehr rasch, war aber 
Beeresslich erst nach 7 min auf den Wert der ungestérten Kiihltruhe gesunken. 
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|): ace Einfiihren der Feuchte; J = Impfen; 
B = Beobachten; 1. = erste Eiskristalle; 
v. = vollstandige Vereisung. 


Figur 1 
Temperatur-Diagramm zu Impfexperiment Cus. 


- Da die Impfung in den ersten 2 min nach der Nebelerzeugung erfolgen musste, 

-konnte daraus geschlossen werden, dass trotz der guten Anpassungsfahigkeit 
der verwendeten Thermometer die Hauchmethode sich fiir Keimfahigkeitsunter- 
suchungen nicht eignete. 

Mit Hilfe eines Mikrofliissigkeitszerstaubers war es moglich, den Nebel durch 
direkte Zerstiubung des Wassers zu erzeugen. Infolge der geringen Luftmenge, 
die zur Zerstaubung notig war, und des Wegfalls der Kondensationswarme gelang 
es, die Temperaturschwankung auf etwa 3°C zureduzieren. Damit war eine gute 
Erfassung der Temperatur im Moment der Vereisung gewahrleistet. Um tiber die 
Temperaturverhaltnisse wahrend jedes Impfversuchs genau orientiert zu sein, 
wurden jeweils 4 Temperaturgradientmessungen rasch hintereinander vorgenom- 
men: in ungestérten Verhaltnissen ; nach der Nebelerzeugung; nach der Impfung 
und nach der Beobachtung. Figur 1 illustriert den gemessenen zeitlichen Tempe- 
raturgang in der Kiihlkammer wahrend eines Impfversuches mit CuS, aufge- 
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nommen in verschiedenen Héhen der Kammer!). Der vor der Feuchteeinfiihrung 
und nach der Impfung im Bilde deutlich wahrnehmbare allgemeine Temperatur-_ 
anstieg zeigt das Temperaturverhalten wahrend den Arbeitspausen des Kiuhl- } 
aggregates. ; 

Die Auswahl der Substanzen erfolgte zunachst rein kristallographisch-chemisch, — 
spater hauptsdchlich in Hinsicht der Polarisierbarkeit der Gitterbausteine. : 

Um méglichst genaue Kenntnisse von der chemischen und kristallographi- 
schen Reinheit der vorliegenden Substanzen zu erhalten, wurde ein Drittel der 
Impfsubstanzen réntgenographisch untersucht. Unter diesen befanden sich alle | 
die Substanzen, die im Temperaturintervall 0° bis —10°C Eiskeimfahigkeit © 
zeigten. Die Untersuchung ergab, dass keine chemischen Verunreinigungen vor- 
lagen. Dagegen zeigten sich bei einigen Substanzen kristallographische Verun- — 
reinigungen. Im Falle des AgI- und ZnS-Praparates, in denen zwei kristallogra- 
phische Modifikationen zugleich anwesend waren, gelang es aber leicht, die beiden 
Modifikationen zu erkennen und wenigstens approximativ ihre Mengenverhdalt- 
nisse anzugeben. 

Zur Gewdhrleistung der Sauberkeit der Oberflachen waren nur frische und | 
gut getrocknete Substanzen angewandt worden. Es wurde ferner darauf geachtet, 
Praparate zu erhalten, die unmittelbar vor der Einstéubung in den Nebel im 
Morser zerkleinert werden konnten, so dass frische Spaltflachen vorlagen. 

Teilweise wurden die Impfstoffe selber laboratoriumsmassig dargestellt; 
diese umfassen die Stoffe: CuS, MnS, Ag.S, PbS, CoS, CdSe, ZuSe, Ag,Se, CuSe, 
HgTe, CdTe, InSb, 5iMg,, Cu,O, Ag,O, Ag,Cr,O,, Cul, CuBr, CuCl, TU, PbL,; 
PbCl,, HgI,, AulI, AgCl, AgI. Die iibrigen Stoffe wurden vom B.D.H. Labora- 
tory, London (British Drug Houses Ltd.), von Bender & Hobein, Ziirich (Ver- 
tretung der Merck & Co.), und der Chemischen Fabrik Fluka in Buchs (SG) be- 
zogen. 


C. Messergebnisse 


Die Wirkungstemperaturen, welche fiir die untersuchten Impfsubstanzen ge- 
funden worden waren, sind in Tabelle 1 zusammengestellt. In der ersten Zahlen- 
kolonne sind die Temperaturen jener Kiihlkammerhoéhe angegeben, bei welcher die 
ersten Eiskristalle in Erscheinung getreten sind, in der zweiten Kolonne die Tem- 
peraturen jener Hohe, auf welcher vollstandige Vereisung des Nebels erfolgt war. 

Die beobachteten Wirkungstemperaturen lassen sich in fiinf Gruppen ein- 
orduens, +4" bis =5,52;.—5,5° bis =7,5°:; —10° bis —14°% —14° pis — tye 


1) Fir den zeitlichen Temperaturverlauf in der durch Hinfiihrung der Feuchte, das heisst der 
Trépfchenluft, gestérten Kiihlkammer kann auf theoretischem Wege die Formel 


eek Zs kr eK] 


hergeleitet werden. Hierin bedeuten T die wahre Temperatur der Luft in der Kiihltruhe, Tdi 
Temperatur des Thermometers, Ty jene der Kihltruhenwand (hier als konstant vorausgesetzt), 
I, jene der eingefiihrten Trépfchenluft und t die Zeit. k, bzw. kx sind Konstante, die durch die 
Ausdrticke 


lg2 lg2 
ky = S bzw. kk = =. 
UR 


T* 


gegeben sind, worin a bz w.t% die Halbwertzeiten fiir die Angleichung der Temperatur des 
Thermometers an die Lufttemperatur bzw. der Lufttemperatur an die Temperatur der Kiihl- 
truhenwand darstellen. 
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Wirkungstemperaturen dey untersuchten Impfsubstanzen 


Silberjodid (f-, y-Agl) 
Silbersulfid (Ag,S + Glas) . 
Cuprojodid (Cul) 
Cuprisulfid (CuS) . 
Cupriselenid (CuSe) . . 
Kadmiumtellurid (CdTe) 
Bleijodid (PbI,) 
Vanadiumpentoxyd (Ve O oe 
Mercuritellurid (HgTe) 
Silberoxydul (Ag,O) 
Kupferoxydul (Cu,O) . 
Kadmiumselenid (CdSe). 


Silbernitrat (AgNO, ce : 


Quarz (SiO,). 

Porzellan . 

Grammatit : 

Olivin (Peridot) 

Ibgleiaal a ae 

Goldjodid (Aul) 
Kupferoxyd (CuO) . 
Zinkselenid (ZnSe)... . 
Silberchlorid (AgCl -+ Glas) 
Mercurochlorid (Hg,Cl,) . 
Mangandioxyd (MnO,) 
Mellur (fe): 
Kaliumjodid (KI- -Rauch) 
Natriumjodid (NaI-Rauch) 
Silberbichromat (AB:Cr101) 
Jodoform (CHI;) 

B-Zinn (B-Sn) 

Kieselgur . 

PNICEIVKOULG Sy es 
Bleichlorid (PbCl,) . 
Mercurijodid (HgI,) 


Silberselenid (Ag ,Se + Glas). 


Indiumantimon (InSb) 
Kadmiumjodid (CdI,) 
Stannooxyd (SnO) . 
Berylloxyd (BeO) . . 
Thallojodid (TH + Glas) 
Kadmiumsulfid (CdS). 


Kali-Alaun (KAI(SO,),.°12 H,O) 
Ammonium-Alaun (NH,A1l(SO,)2 


Nicht wirksam bis —18° sind: 
Mg(OH),, Ca(OH)s, 


SrF,, BaF, MgF;, CaF,, 


Erste Eis- Vollstandige 
kristalle°C Vereisung °C 


= 4:0 are 
ey A? = 5,3 
Fe = 15,5 
46 = 14 
242 aE a 
att a 5.1 
a5 2 56,5 
og YY | = '6:6 
255-0 767 
= G2 Sey 2 
eS. 1. =, 631 
=" 5,4) =65 
26:3 770 
A773 =473,0 
2116 Pa ea) 
= 9,0 = 1:250 
123 — 13,5 
ltrs 35 
a = 10,0 
Sg = 1354 
S1e3 ee 
—t,0 11350 
Mire et: pa 
att 235 = 4,0 
=i2,3 — 13,6 
5120 +133 
Sse —ais.5 
= 12:0 14,0 
S020 15,0 
12-0 ee a 
— 10,5 or act 
iors SAO 
S120 =13;9 
BE ey S140 
S125 243.6 
= Y5;5 e170 
=< 15% Lely Ss 
1 1355 = 05,5 
— 14,0 13,0 
S143 1750 
4 46,0 
scare ae dak — 14,6 Zi70 
-12H,0) . . — 15,0 oF 7170 


CdF,, ZnF,, NaF, LiF, ZnO, Sn0,, 


CaO, BaO, SrO, MgO, Al,O,, CuCl, CuBr, HgBr,, NH,Cl, NH,Br, NH,I, kL, Zi, 
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PbS, CoS, MnSgriin, NH,F, NaCl, KCl, Na,SO,-10 H,O, Samorph> Siaist. SL, SiC] 
Biotit, Orthoklas, SiMg,, Al, Cu, Sb, SrSO,, PbSO,, BaSO,, CaSO,, MnSOQ,°1 HO, 
MnCO,, FeCO,, BaCO , ZnCO;, SrCO;, CaCOz, CdCO,, PbCO;, MgCO;, Graphit, 


Koks, Koksasche, Zigarettenasche, Glas. 
Total 104 Substanzen, davon 61 inaktiv und 43 aktiv. 


Wirkungstemperatur unterhalb — 18°. Wird die Anzahl der aufgetretenen Ver- 
eisungen (je eine Vereisung pro Substanz gerechnet) in Abhangigkeit von der 
Temperatur graphisch aufgezeichnet, so ergeben sich fur das Auftreten von ersten 
Eiskristallen Maxima bei —4,5° und —12°, fiir das Auftreten von vollstandigen 
Vereisungen Maxima bei —5,7° und —13,5°. 

Wahrend der Keimfahigkeitsexperimente wurde folgendes Verhalten der 
Impfsubstanzen beobachtet: iy 

1. Wasserlésliche oder durch Wasser zersetzbare Substanzen wirkten bis —18° 
tiberhaupt nicht (CaO, SrO, NH,F, KCl, NaCl, NH,Cl, NH,Br), oder aber sie 
wirkten erst bei Temperaturen tiefer als —12° (CdlI,, KI, Nal). 

2. Wahrend einige Substanzen eine grosse Latenzzeit aufwiesen (PbI, rund | 
30s), zeigte sich bei den meisten sofortiges Einsetzen der Vereisung (5 bis 10s). 

3. Es war nicht festzustellen, dass keimfahige Substanzen besonders trocken 
sein mussten. 


D. Ubersichtsmassige Erérterung der Befunde 


a) Einfluss des kristallinen Aufbaus auf die Eiskeimfahighkeit 


Ordnet man die untersuchten Impfsubstanzen rein nach kristallographischen 
Gesichtspunkten, so kénnen folgende Aussagen gemacht werden: 

1. Die Substanzen lassen sich in 44 Gittertypen einordnen. Obwohl gewisse 
Gittertypen, wie zum Beispiel der Fluorittyp oder Casium-Chlorid-Typ unter den 
aktiven Substanzen nicht vertreten sind, andere, wie zum Beispiel der Rutil- oder 
Steinsalztyp, nur sehr schwach, der Zinkblendetyp dagegen eher starker vertreten 
ist, zeigt sich keine besondere Bevorzugung irgendeines bestimmten Kristallbau- 
schemas. 

2. Auch die Symmetrieverhaltnisse scheinen keine entscheidende Rolle als 
Ursache fiir die Eiskeimfahigkeit zu spielen. Von den 44 keimfahigen Substanzen 
(B-AgI und y-AgI als je eine Substanz gerechnet) weisen nur 10, also etwa 23%, 
hexagonale oder trigonale Symmetrie auf. Die restlichen 77 % zeigen andere Sym- 
metrieverhaltnisse: 14 Substanzen sind kubisch, 5 tetragonal, 13 orthorhombisch, 
monoklin oder triklin, und 2 sind amorph. 

Betrachtet man dieses Kriterium nur innerhalb der Substanzen, die im Tem- 
peraturgebiet von 0° bis —10° keimfahig sind, so ergibt sich folgende Verteilung 
auf die 8 vorliegenden Gittertypen: 4 Substanzen sind hexagonal oder trigonal, 
7 sind kubisch und 3 sind orthorhombisch. 

3. Untersucht man die keimfahigen Substanzen beziiglich der Teilchenab- 
stande, die sich auf einer Spalt- oder Wachstumsflache auswirken kénnen, so ge- 
langt man zu einem interessanten Ergebnis: 

Berechnet man fiir keimfahige Substanzen, welche die kubische Zinkblende- 
struktur besitzen, also ein von Eis sehr abweichendes Gitter und andere Gitter- 
konstanten aufweisen, den sehr oft auftretenden Teilchenabstand a,/V'2, so gelangt 
man zu Werten, die jetzt dem Wert fiir die Tetraederkantenlange im Eisgitter 
(4,51 A) sehr nahe kommen. Ahnliche Beobachtungen kann man auch bei der 
Berechnung von besonderen Abstanden an anderen Gittern machen. 
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Man gelangt schliesslich zum Resultat, dass von den 44 aktiven Substanzen 
2/3 Abstande aufweisen, welche innerhalb einer 7 %igen Abweichung des Abstan- 
des 4,51 A im Eisgitter liegen. Unter den restlichen Substanzen sind noch 50%, 
deren Abstande unbekannt sind, das heisst, es ist durchaus méglich, dass auch 
‘diese giinstige Abstande aufweisen. So besitzen also im besten Falle noch 1/6 der 
-Substanzen von Eis abweichende Abstande. 

Betrachtet man nur die 14 aktivsten Substanzen (—4° bis —7°), so ergibt sich, 

dass fiir 12 Substanzen die Abweichungen der Abstande gegeniiber dem Abstand 

4,51 A im Eisgitter nur noch maximal 4,5 % sind. Es verbleiben noch 2 Substan- 
zen, deren Abstande stark differieren: AgNO, und V,O, (V,O, weist ebenfalls einen 
giinstigen Abstand auf, dieser kann sich aber nur bei besonderen, kaum mog- 
lichen Spaltverhdltnissen auswirken). Bemerkenswert ist, dass bei CuS und be- 

‘sonders CuSe der doppelte Wert des Abstandes Cu—S = 2X 2,19A = 4,38 A bzw. 
Cu—Se = 2x 2,275 A = 4,55A in der (0001)-Flache nur eine geringe Differenz 
zum Abstand 4,51 A im Eisgitter zeigt. 

Dehnt man die Berechnungen der Teilchenabstande, die sich auf Spalt- oder 
“Wachstumsflachen auswirken kénnen, auf die als inaktiv gefundenen Substanzen 

aus, so gelangt man aber auch hier zum Ergebnis, dass noch viele Substanzen 
-giinstige Abstande aufweisen. Von den 61 Substanzen besitzen 31% Abstande, 
deren Abweichungen geringer, ja sogar viel geringer als 7% sind. 

Vergleicht man noch die Impfsubstanzen InSb, CdTe, AgI (isoelektrische 
Reihe) und HgTe, die alle denselben Gittertypus, dieselbe Symmetrie und sehr 
ahnliche geometrische Abmessungen aufweisen, so miisste die Reihenfolge der 
Keimfahigkeit sein: 


HgTe, besser als CdTe = InSb, besser als Agl. 


Die experimentellen Keimfahigkeitsuntersuchungen zeigen aber, dass die 
Reihenfolge eine ganz andere ist: 


CdTe = AglI, besser als HgTe, viel besser als InSb. 


Uberlegt man sich den Fall von orientiertem Aufwachsen von Eis auf einer 
Impfsubstanz, so muss noch folgende Erscheinung beriicksichtigt werden: 

Vergleicht man zum Beispiel die Teilchenabweichungen, die sich ergeben 
wiirden beim Aufwachsen einer Eisgittergeraden auf einer Gittergeraden auf der 
(0001)-Flache der sehr gut keimfahigen Impfsubstanzen CuS oder AglI, so gelangt 
man zum Resultat, dass die Abweichungen zu den nachsten Nachbarn anfanglich 
wohl giinstig sind (besonders bei AgI), sich aber dann verschlechtern, um erst 
nach relativ grosser Periode wieder besser zu werden (besonders bei AgI ausge- 

- pragt). Bei der bis —18° nicht keimfahigen Substanz ZnS zeigt sich auf der (0001)- 
Flache anfanglich eine relativ grosse Abweichung, die aber sehr rasch auf einen 
geringen Wert vermindert wird. 

Es wird aus diesen Uberlegungen klar, dass nicht nur die Abweichungen des 
Teilchenabstandes 4,51 A im Eisgitter von einem bestimmten Elementarabstand 
auf dem Impfkristall betrachtet werden diirfen, sondern dass das gegenseitige 
Verhalten der Teilchenlagen wahrend des ganzen Aufwachsvorganges studiert 
werden muss. 

Es ist klar, dass solche Vergleiche nur approximativ méglich sind. Fir nahere 
Vergleiche miissten die betreffenden Gitterabstande viel genauer bekannt sein. 

Zusammenfassend muss festgestellt werden, dass zwischen den Teilchenabstan- 
den auf Spalt- bzw. Wachstumsflachen von Impfsubstanzen und der Eigenschaft 
dieser Substanzen, den Phasentibergang von fliissigem, unterkiihltem Wasser in Eis 
zu bewirken, kein unmittelbarer primarer Zusammenhang gefunden werden kann. 
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b) Stérungsversuche der Eiskeimfahigheit 


In einer besonderen Gruppe von Experimenten wurde versucht, die Aktivitat 
von eiskeimfahigen Substanzen zu stdéren. Es ergaben sich dabei folgende 
Resultate : 7 | 

1. Wie an Hand der Réntgendiagramme festgestellt werden konnte, weist 
unter den hergestellten AgI-Préparaten eines die kubische y-Modifikation und 
eines die hexagonale f-Modifikation auf. Das kubische Praparat (Fig. 2) wurde 


Figur 2 
Debey-Scherrer-Diagramm von AgI aus konz. AgNOg, kubische Struktur. 


durch Eingiessen einer KI-Lésung in eine konzentrierte Silbernitratlésung er- 
halten; die hexagonale 6-Modifikation (Fig. 3) durch langsames Auskristallisieren- 


Figur 3 


Debey-Scherrer-Diagramm von AgI aus HI, hexagonale Struktur. 


lassen einer Lésung von AgI in HI, welche mit Alkohol iiberschichtet war. Ob-_ 
wohl der Reinheitsgrad beider Substanzen in kristallographischer Hinsicht nicht 
volistandig 100% war, hatte doch eine starke Verminderung in der Zahl der aus- 
gelésten Eisteilchen auftreten miissen, falls eine Modifikation nicht eiskeimfahig 
ware. Es stellte sich aber klar heraus, dass beide Praparate im normalen Wirkungs- 
bereich von AgI (—4° bis —6°) die Vereisung hervorriefen. Dasselbe Resultat ergab 
sich bei den Praparaten, wo Gemische beider Modifikationen vorlagen, und im 
Falle von AgI-Rauch. 

2. Sehr uniibersichtlich wurden die Verhaltnisse bei der Untersuchung des 
Einflusses-von Jod auf die Wirksamkeit der Impfsubstanzen. 

Jod selbst musste nach mehrfacher Untersuchung schliesslich als nichtkeim- 
fahig bis —18° angesehen werden. Es war bei diesen Versuchen sehr grosse Vor- 
sicht in der Beobachtung anzuwenden, um nicht leicht zu Fehlschliissen zu ge- 
langen. Das direkt in der Kiihlkammer sublimierte Jod ging sofort in die feste 
Form tiber, so dass ein «Regen» von feinsten Jodkristallen einsetzte. Diese ver- 
mochten den vorhandenen Nebel nicht zu impfen. Obwohl entstandeneEiskristalle 
nur sehr schwer von den Jodkristallen unterschieden worden waren, hatte doch 
der Nebel infolge der Vereisung sehr rasch verschwinden miissen. Dies trat aber 
nicht ein, der Nebel blieb sogar oft so lange erhalten, bis die Jodpartikel ausgefal- 
len waren. 

War Cul bei der Herstellung nicht vollstandig reduziert worden, oder wurde 
reines, weisses Cul mit I, zusammen pulverisiert, so gelang es, die Eiskeimfahig- 
keit dieser Substanz innerhalb des Temperaturintervalls 0° bis —12° zu zerstOren. 
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Rontgendiagramme zeigten keine neuen Interferenzlinien. Das dem Praparat 
beigemengte Jod musste somit adsorbiert worden sein und lag nicht kristallin Vor. 
‘Wurde dagegen AgI mit Jod zusammen pulverisiert, so konnte keine St6rung 
der Impfeigenschaften festgestellt werden. 

Um das Verhalten von AgI und Cul noch intensiver zu studieren, waren aus 
‘diesen zwei Komponenten drei verschiedene Mischkristalle hergestellt worden 
mit 7,3; 9,1 und 91 Gewichtsprozenten AglI. Die Mischkristalle waren aus HI- 
Lésungen auskristallisiert und réntgenographisch verifiziert worden. Es ergab 
sich mit zunehmendem Cul-Gehalt ein Absinken der Wirkungstemperatur von 
—4,0° bis —7,5°. Dieses Verhalten muss erneut auf den Einfluss von Jod zurtick- 
gefiihrt werden. 


c) Priifung des Weylschen Eiskeimfdhigheitsmechanismus 


Es ist schliesslich versucht worden, die Giiltigkeit des von WEYL!) postulierten 
Impfmechanismus abzuklaren. Kristalline Kérper weisen je nach Gitterverhalt- 
nissen und Polarisierbarkeiten der Gitterteilchen spezielle Oberflachenstrukturen 
auf. Nach den Weylschen Vorstellungen sind solche Oberflachen bisweilen 
imstande, die in der Oberflache eines Wassertropfens vorhandene elektrische 

-Doppelschicht (das negative Ende der Wasserdipole weist nach aussen, das posi- 
tive Ende nach dem Tropfeninneren), von welcher angenommen werden muss, 
dass sie den Vereisungsprozess stark hemmen wird, bedeutend zu schwachen. Auf 
dieser Eigenschaft soll die Eiskeimwirkung eines Stoffes beruhen. 

Bei der Priifung der untersuchten Substanzen beziiglich der Polarisierbarkeit 
der am Gitteraufbau beteiligten Ionen ergaben sich folgende Verhaltnisse: 

Als keimfahig erwiesen sich vorzugsweise Impfsubstanzen, welche Anionen 
mit grosser Polarisierbarkeit im Gitter enthalten: O', eS <er be 5d = 
stanzen mit F~ oder Br~ wirken iiberhaupt nicht, solche mit Cl~ nur teilweise. 

Wurde in einer Verbindung das Kation konstant gehalten und das Anion 
durch ein anderes ersetzt, so ergab sich, dass mit der Zunahme der Polarisierbar- 
keit des Anions die Keimfahigkeit der betreffenden Substanz sich verbesserte. 


Zum Beispiel: 


Cdk, CuCl AgCl | HgBr, PbCl, NaCl | th ZnS | 
Cdl, CuBr Ag,O Hel, Phig-4 Nal iy KI LDC IA) 
CdS CuO { Agl HgTe 
CdSe Cu,O | Ag.S 
CdTe | , Cul 
: Cus 

CuSe 


¥ 
Im iibrigen verweisen wit auf den II. Teil dieser Arbeit. 


d) Allfallige Auswivkung auf die W olkenimpfuersuche 


Das Auffinden von neuen Substanzen, die sich im gleichen Temperaturinter- 
vall als keimfahig erwiesen, in welchem auch AgI aktiv ist (—4° bis —6°), offnet 
neue Perspektiven fiir die kiinstliche Auslésung von Niederschlagen in der freien 
Atmosphare und fiir die Bekampfung der Hagelbildung, in welchen Versuchen 
grésstenteils AgI als Impfstoff verwendet worden ist. AgI besitzt aber den Nach- 
teil, dass es teuer und relativ unbestandig ist (photolytische Zersetzung unter 


1) W. A. WEVL, J. Coll. Sci. 6, 389 (1951). 
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dem Einfluss von UV.-Licht unterhalb 4400 A). Von den neu gefundenen Sub- 
stanzen zeichnet sich neben CdTe und Cu,O besonders CuS aus. Letztere Impf- 
substanz hat die gleiche Keimfahigkeit wie AgI, ist aber aus sehr billigen Aus-— 
gangsmaterialien in einfacher Weise darstellbar und zudem sehr bestandig. 


In diesem ersten Bericht haben wir uns auf die Schilderung der angewandten 
Untersuchungsmethoden und erhaltenen Ergebnisse beschrankt und sind auf 
Erklarungsversuche im grossen und ganzen nicht naher eingegangen. Dies soll in 
einer zweiten Notiz geschehen, in welcher versucht wird, durch Erweiterung der 
Weylschen Theorie eine Erklarung der Beobachtungen zu geben. 

Es ist unsere angenehme Pflicht, der Werkzeugmaschinenfabrik Oerlikon, 
Biihrle & Co., Ziirich-Oerlikon, fiir die zur Verfiigung gestellten finanziellen Mittel, 
Herrn PD. Dr. W.Epprecut fiir die gelegentlichen Ratschlage, Herrn Prof. H.WE- 
BER, Institut fiir Fernmeldetechnik, fiir die Bereitstellung geeigneter Raumlich- 
keiten, herzlich zu danken. Ferner méchten wir auch der Migros-Genossenschaft 
Ziirich fiir die Uberlassung einer Kiihltruhe «Hill» und der Firma Trion fiir die 
Uberlassung eines Elektrofilters unseren Dank aussprechen. : 


Summary 


An apparatus has been developed which enables the temperature conditions 
to be studied in a freezer before, during, and after a seeding procedure with an 
inaccuracy of less than 1°C. With these means it was possible to fix with a fair 
degree of accuracy the temperature at the moment of phase transition of the 
liquid undercooled water droplets to ice and to study the dependence of this 
temperature on the properties of the seeding substance. An illuminating system 
suitable for the purpose, an improved method for the fine dispersion of the 
seeding substances and for the formation of clouds, and finally two carefully 
tested temperature measuring systems formed the basis for this result. The icing 
nucleability of altogether 104 different substances in dependency on their crystal- 
lographic properties and the special arrangement of their grid elements in the 
surface have been examined. No distinguished crystallographic influence has 
been found. 


(Eingegangen: 26. Mai 1955.) 


Note on the Torsional Rigidity of Cylinders of Circular 
Sector Cross-Sections 


By Per-Lin SHENG, Bloomington, Indiana, USA?) 


Numerical values for the torsional rigidities of certain circular sectors have 
been given by St. VeENANT some 100 years ago [5], [4], [4A]?); in general these 
values are inaccurate (cf. [3], [2]). Simple explicit results for four cases have been 
given by STEVENSON [3]. The present author published a short note [1] on this 
subject before learning of StEVENSON’s prior paper. Three of STEVENSON’s cases 


1) Indiana University, Graduate Institute for Mathematics and Mechanics. 
2) Numbers in brackets refer to References, page 419. 
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are also re-obtained recently by Herzic[7]. The present paper gives 16 cases 

obtained along the same elementary lines as in [1] by using either WHEELON’S 

summation method [6] or a classical elementary summation method (cf. [2], 
where the eight cases in which the sector angle is a multiple of 2/4 are worked 
‘out in detail and other cases are indicated). Still other cases can be obtained 
similarly, although it is probable that more than one L-function will be needed. 

As is well-known, the torsional rigidity coefficient K detined by T= wa a” i 

‘is given by 


2 
K= 2 f[waa, (1) 
A 


' where T is the transmitted torque, # the shearing modulus, « the angle of twist 
per unit length, a the radius of the sector, A the cross-sectional area, and ¥Y the 
St.Venant stress function for torsion which vanishes at the boundary. After 
substituting the well-known expression for W into (1) and performing the inte- 
gration we get 


; 512 k4 = 1 
ae = -suek® HG 
K qz (tan2ka 2k x) = wae ates ame a 


where 2k a is the sector angle (hence 0< <1). Splitting the fraction after 
_ the sign of summation into partial fractions, not of the first order in ” as usual 
(this leads to STEVENSON’S formulation), but all of the second order, 1. e., 


1 1 
(n— 4k) n® (n+ 4k)? 256 k3 (3) 
il 5 . ¢42 % v) | 
Osu | net m2 (n+ 4h)? )’ 
we get 
K= i (hat tan 2 k zm) 
(4) 
2k = 1 - 1 — 1 
> path > ee) 2 EN oe ees 
ae rien ee (n sik ) i ee ke! (n 1g 4 R) Ter Oe “ae! k)* 


The split (3) is not unique; other similar decompositions can also be used, giving 
equivalent results. The first two series in (4) can be summed either by a classical 
method or by WHEELONn’s method, while the third one can be expressed in terms 
of L-functions (cf. [2]). By WHEELON’s method [6], (4) can also be transformed to 
several integral expressions. As an example, we have 


K= (2ka+tan2km)— 37> ° Gah G— 4h) 


1 4 k (5 — 8k) 
‘ | 


k > 4 ( 4k) 4k > (4k-+1)e (5) 
f i) eee gs w* (1 — x xe 
er Bf aera — [ax ae aon dy 6 

0 0 0 
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The results, obtained by using (5) except for the last term, for which the last 


_ term (in the brackets) in (4) is substituted, are presented in the table, where 


eae 1 1 
oe eee gh ey NPCS a «5 
ye age 1 ee 
= 52 = 72 howe t= 02976) oh. (6) 
i 1 1 
w=1 ie Bf OUT BL, ais. « 


All these are L-functions of the form L(2, x); in particular, s is the well-known 
CATALAN’s constant. The two cases for k = 1/4, 3/4 are obtained by a limiting 
process, since ¥Y% has no meaning in such cases. 
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Zusammenfassung 


Fiir 16 besondere Kreisausschnittwinkel wird die Drillungssteifigkeit eines 
Stabes mit Kreissektor-Querschnitt durch L-Funktionen ausgedriickt. 


(Received: March 25, 1955.) 
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Friihjahrstagung der Schweizerischen Physikalischen Gesellschaft 
vom 7. bis 8. Mai 1955 in Vevey 


Berichte tiber angewandte Physik und Mathematik 


Zur Messung des Rauschens von Transistoren, von E. BALDINGER und 
F. LEUENBERGER, Basel?). 

Transistoren erzeugen einen Stérspannungsuntergrund, der eine Grenze dar- 
stellt beziiglich der sinnvollen Verstarkung kleiner Signale. Vom Standpunkt der ~ 
Anwendungen aus interessiert uns die Kenntnis folgender Gréssen: Rauschspan- 
nung des Emitters bzw. des Kollektors und die Korrelation dieser Rauschspan- 
nungen. Die in Figur 1 skizzierte Messeinrichtung dient zur Messung dieser 
Grossen. 

Die beiden Signale werden verstarkt, und, nachdem sie auf gleichen Pegel 
gebracht worden sind, auf eine Additionsstufe (B) gegeben. Das Ausgangssignal 
derselben wird iiber ein Filter (C) einem linearen Detektor (D) zugefiihrt. Mittels 
der Phasenumkehrstufe (A) kénnen die Signale zur Bestimmung des Korrelations- 
koeffizienten sowohl addiert wie auch subtrahiert werden. Zwei Teile der Mess- 
schaltung, die erste Verstarkerstufe und der lineare Detektor, sollen hier kurz 
skizziert werden. 

Bei der Verstaérkung kleiner Signale im Tonfrequenzgebiet hat man fiir jeden 
Verwendungszweck die geeignete Réhre wie auch deren optimale Betriebsbedin- 
gungen zu finden?). Die Réhre E80F (Philips) weist einen geniigend kleinen 
aquivalenten Rauschwiderstand auf und ist ausserdem sehr wenig mikrophonisch. 
Das gemessene Rauschspektrum’) des besten Exemplars aus einer Anzahl von 25 
liess sich darstellen durch Ry = 1670 + 750/f. Der 1/f-Term entspricht dem Fun- 
keleffekt. Bei 1000 Hz betragt der Rauschwiderstand 2,35 RQ, ein Wert, der nur 
um den Faktor 1,35 hdher liegt als der asymptotische Wert tiir hohe Frequenzen 
(Schroteffekt). 

Ublicherweise werden bei Rauschmessungen quadratische Detektoren ver- 
wendet, da diese unmittelbar den quadratischen Mittelwert U.;; messen. Der bei 
unseren Messungen verwendete lineare Detektor hat gegeniiber einem quadra- 
tischen Réhrendetektor den Vorteil einer erheblich besseren Nullpunktkonstanz. 
Ist die Wahrscheinlichkeit W(U) dU, dass der Spannungswert im Intervall 
U,U + aU liegt, durch 


W(U) dU : 


ee pO eran in 

V2 2 Unsy 
gegeben, dann gilt folgende Beziehung zwischen dem Effektivwert und dem Mit- 
telwert U des linearen Detektors (schraffierte Flache in Figur 3) 

v= Goss 
= V2 It 

1) Physikalisches Institut der Universitat. 
) Vgl. zum Beispiel A. VAN DER Z1EL, Noise (Prentice Hall, New York 1954), S. 230. 
3) Gemessen zwischen 200 Hz und 8 kHz. 


to 
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Figur 1 


U2 =4KTRaf 


100 200 300 400 500kn. | 
Figur 2 
Kontrolle der Apparatur mit Hilfe des Nyquist-Rauschens. 


UCD 


Zeitlicher Verlauf der Rauschspannung. 
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Figur 4 


Prinzipschema des linearen Detektors. 


Anzahl Anzahl 
g 0C71 4 0c70 
9 10 ih 12 WJ 14 15 db a 10 | 12 13 14 15 db 
Rauschzahl Rauschzaht 
Figur 5 


Verteilung der Rauschzahlen. Frequenz: 1000 Hz. Effektive Bandbreite: 45 Hz, Emitterstrom: 
2mA. Kollektorspannung: —2 V. Quellenwiderstand R, 3500/2. T: 295° K. 


Als Test fiir die Zuverlassigkeit der gesamten Apparatur dient die experimentelle 
Uberpriifung (Figur 2) der Nyquist-Formel U2 = 4k T R Af. Hatte man neben 
dem Widerstandsrauschen etwa den Schroteffekt eines Gitterstromes als weitere 
Stérquelle, so wire U® nicht mehr eine lineare Funktion von R. 

Figur 5 zeigt Rauschzahlmessungen an p — n— p-Transistoren der Typen 
OC 70 und OC 71 (Philips). Es wurden je 50 Exemplare untersucht?), 

Fiir die finanzielle Unterstiitzung der vorliegenden Untersuchungen méchten 

wir der Stiftung Hasler-Werke, Bern, unseren besten Dank aussprechen. 


2) Vgl. entsprechende Untersuchungen von W. GuccENBUHL und M. J.O. Strutt, Arch. 
Elektr.-Ubertragung 9, 103-108 (1955) 
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Knickung verwundener Stabe, von E. Hut, Eschlikon. 

Ein gerader elastischer Stab mit zwei verschiedenen Haupttragheitsmomenten 
des Querschnittes, dessen Hauptachsen im spannungslosen Zustand um einen 
festen Winkel je Langeneinheit des Stabes gegeneinander verdreht sind, wird an 
- beiden Enden durch eine in der Langsrichtung wirkende Kraft konstanter Rich- 
tung und konstanten Betrages auf Druck beansprucht. Bei verschiedenen Fallen 
der Lagerung des Stabes an beiden Enden wird die Knicklast berechnet und die 
elastische Linie beim Beginn der Knickung bestimmt. 

Die Arbeit wird als Dissertation (ETH) erscheinen. 

Die Ergebnisse sind bisher durch Experimente nicht iiberpriift worden. 


Entwurf eines elektrischen Analogiegerates fiir den schweizerischen Ver- 
suchsreaktor, von Maurice Martin, Ziirich. 

Es soll ein elektrisches Analogiegerat zur Untersuchung des Betriebsverhal- 
tens des Versuchsreaktors entworfen und gebaut werden. Das Gerat wird zur 
Zeit bei Landis & Gyr in Zug aufgebaut und soll, wenn Messergebnisse vorliegen, 
in einer spateren Ver6dffentlichung genauer beschrieben werden. Hier soll zunachst 
nur die prinzipielle Arbeitsweise erlautert werden. 

Der mit schwerem Wasser als Moderator und Kiithlmedium angefiillte Reaktor- 
tank aus Aluminium befindet sich innerhalb eines als Neutronenspiegel dienenden 
Graphitreflektors und einer 2 m dicken Betonabschirmung. 200 Uranstabe von 
27 mm Durchmesser tauchen in den Tank ein und werden von oben nach unten 
vom Moderatorwasser umspiilt und gekiihlt. Die Kettenreaktion mit U??> wird 
mit Hilfe von Regulierstaben gesteuert, indem durch die Absorption von Neu- 
tronen in Kadmium der Neutronenmultiplikationsfaktor k verandert wird. Im 
Gleichgewichtszustand ist die mittlere Dichte thermischer Neutronen » konstant 
und k = 1. Abweichungen yom Gleichgewichtszustand werden durch die soge- 
nannte Reaktivitat 6k = k — 1 charakterisiert. Das dynamische Verhalten des 
Reaktors muss also durch eine Differentialgleichung beschrieben werden, in der 
die Reaktivitat als Ursache und die Neutronendichte als Wirkung eingehen. Man 
kann dabei von der makroskopischen und der mikroskopischen Neutronendichte- 
verteilung absehen und mit einer mittleren, zur Warmeleistung proportionalen 
Dichte  rechnen, wenn nur die Verteilung leistungsunabhangig angenommen 
werden darf. 

Bezeichnen wir mit ¢, die mittlere Lebensdauer thermischer Neutronen — sie 
betragt etwa 1 Millisekunde — so verschwinden im Mittel n/t) pro Sekunde. Jede 
Kernspaltung erzeugt im Mittel wieder 2,5 neue Neutronen, von denen aber im 
Gleichgewicht nur eines wieder von einem U?#5-Kern eingefangen wird. Ein 
kleiner Bruchteil von etwa 0,8% der neu erzeugten Neutronen erscheint nicht 
sofort bei der Spaltung, sondern mit einiger Verzégerung. Es handelt sich um 
Neutronen, welche beim Zerfall von Spaltprodukten mit Halbwertzeiten zwischen 
0,1 und 80s entstehen. Dieser kleine Bruchteil verzdgerter Neutronen gegentiber 
den prompten ist fiir die Steuerung der Kettenreaktion sehr bedeutend und er- 
leichtert die Aufgabe der Regulierung wesentlich. Ist c; die Konzentration des 
Spaltproduktes mit der Zerfallskonstanten /,, so werden pro Sekunde A; ¢; ver- 
zogerte Neutronen von der Gruppe 7 gebildet. Praktisch kommen sechs Gruppen 
in Frage; dann kommen noch zwei Gruppen verzégerter Neutronen hinzu, welche 
durch Photospaltung des Deuteriums entstehen. Mit einer Haufigkeit 6, fiir die 
Bildung des Spaltproduktes i konnen wir fiir die Konzentration c; schreiben: 


dc; ; : 
M8 Fen he Ss ea 
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Fiir die Anderung der Neutronendichte ” k6nnen wir schreiben: 


a 5 (tees py nape pl pane eta eee 
dt boy matte 


Hierin ist s die Quelle von Fremdneutronen, welche durch spontane Spaltung und 
aus der kosmischen Strahlung entsteht. 

Fir den Entwurf einer elektrischen Analogieschaltung zur Nachbildung der 
reaktorkinetischen Gleichungen liegt es auf der Hand, die Aufladung eines Kon- 
densators zu beniitzen und als Mass fiir die Neutronendichte die Kondensator- 
spannung zu wahlen. Der Ladestrom muss der rechten Seite der jeweils nachge- 
bildeten Gleichung entsprechen. In der Differentialgleichung fiir die Neutronen- 
dichte ” ist die Produktion prompter Neutronen der Dichte m selber proportional, 
das heisst, in der Simulatorschaltung muss der entsprechende Beitrag zum Lade- 
strom der Kondensatorspannung proportional gemacht werden. Um dies zu er- 
zielen, wird die Ladekapazitat so zwischen Eingang und Ausgang eines linearen 
Gleichstromverstarkers geschaltet, dass die Kondensatorspannung gleich der 
halben Ausgangsspannung wird. Am Ausgang hat man dann symmetrisch zum 
Erdpotential eine positive und eine negative Spannung, die beide zur Neutronen- 
dichte ” proportional sind und zur Speisung eines schnell schreibenden Registrier- 
instrumentes herangezogen werden kénnen. Der Reaktorgleichung entspricht 
nunmehr die elektrische Ladegleichung: 


du ai : 
Ter aero: 


Es muss nun eine Schaltung gewahlt werden, mit Hilfe derer wir die aus einem 
stabilisierten Gleichspannungsnetz entnommenen Stréme leistungslos so steuern 
konnen, dass sie den Gréssen proportional sind, die in den reaktorkinetischen 
Gleichungen vorkommen. Hierfiir eignet sich eine Kathodenfolgerschaltung, bei 
der eine Triode so zwischen + 300 V geschaltet wird, dass eine an richtiger Stelle 
des Kathodenwiderstandes abgegriffene Spannung mit recht guter Genauigkeit 
der Gitterspannung proportional ist. Die Proportionalitatskonstante ist nahezu 1. 
In einem von Ber und Srrauss [Rev. Sci. Instr. 27, 760 (1950)] beschriebenen 
Pile-Simulator wird die Doppeltriode 6SU7GTY verwendet, welche eine mit 
hochisolierendem Sockel ausgeriistete, ausgesuchte 6SL7 ist. 

Die gewtinschte Reaktivitat (prompte Neutronen) wird an einem zwischen den 
Ausgangsklemmen des Verstarkers liegenden Helipot-Potentiometer eingestellt. 
Die verz6gerten Neutronengruppen werden durch eine entsprechende Anzahl 
gleichartiger Schaltelemente nachgebildet, bei welchen im Gitterkreis des Katho- 
denfolgers ein R — C-Glied mit richtiger Zeitkonstante der Gruppe mit der Zer- 
fallskonstanten A; entspricht. Die Dimensionierung der Widerstande und der 
Kondensatoren ergibt sich durch einen Vergleich der reaktorkinetischen und 
elektrischen Gleichungen. 

Die Nachbildung des Einflusses der Temperaturaénderungen in den Uranstaben 
und im Moderator auf die Reaktivitat und ebenso derjenige der Vergiftung des 
Urans durch das die Neutronen sehr stark absorbierende Spaltprodukt Xenon 
(Xe135) geschieht durch ahnliche Schaltungen wie fiir die Nachbildung der promp- 
ten und der verzégerten Neutronen. Allerdings miissen Produkte zwischen zwei 
Gréssen bzw. Spannungen gebildet werden k6nnen, was mit Hilfe von elektro- 


mechanischen Multiplikationssystemen kleiner Reaktionszeit bewerkstelligt 
werden kann. 


E. 
: 
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Mit Hilfe dieses Pile-Simulators soll die fiir die Regulierung wichtige Uber- 
_tragungsfunktion ermittelt werden. Sodann kann ein Modell fiir den Servomecha- 
-nismus der Regulierung damit untersucht werden. Schliesslich erméglicht die 

Untersuchung des zeitlichen Verhaltens der Neutronendichte bei verschiedenen 
- Reaktivitaitsanderungen die Aufstellung eines Programms fiir das Anlassen und 
Abstellen des Reaktors, was jeweils mit geniigender Vorsicht zu geschehen hat. 


Jahrestagung GAMM 1956 


Die Gesellschaft fiir Angewandte Mathematik und Mechanik (GAMM) ver- 
anstaltet ihre nachste wissenschaftliche Jahrestagung in der Zeit vom 22. bis 
26. Mai 1956 (Pfingstwoche) in Stuttgart. Die drtliche Tagungsleitung liegt in 
Handen von Professor Dr. GUNTHER ScHULZ, Mathematisches Institut A der 
Technischen Hochschule Stuttgart. H. GORTLER 


Erratum 


-Asymptotische Entwicklungen fir Fresnelsche Integrale und ver- 
wandte Funktionen und ihre Anwendungsm6glichkeiten bei der Berech- 
nung spezieller Raketenbahnen. Von JosEF ZBORNIK, Chur (ZAMP 5, Fasc. 4, 
345 [1954]). 


Wie P. Henrici (ZAMP 6, Fasc. 2, «Kleine Bemerkungen. . .») hervorhob, 
weist der Aufsatz Asymptotische Entwicklungen fir Fresnelsche Integrale..., 
ZAMP 5, Fasc. 4, 345ff. (1954), zwei Inkorrektheiten auf. Zu ihrer Berichtigung 
sind folgende Korrekturen vorzunehmen: 


Es ist zu streichen: S. 347, 7. Zeile von unten: «... das Restglied...» 


Es ist zu setzen: 


S, 347, 6. Zeile von unten: 4” ,,(%) anstatt ,,(%)1). 
7 co co 
5. Zeile von unten: lim ad anstatt lim a 
x—>0CO x %—>0O x 
S. 348, 3. Zeile von oben: lim +” Yi anstatt lim JZ,, 
x—>0O %—>0O 
lim #” L, anstatt lim Ly. 
*%—>00 %—>0O 
4. Zeile von oben: #” P,,(*) anstatt f, (7). 
S. 349, 10. Zeile von unten: «Da in den betrachteten...» anstatt «Dain...» 


1) Vgl. zum Beispiel: Lenser, Reihenentwicklungen in der mathematischen Physik, 2. Aufl. 
(Walter de Gruyter & Co., Berlin 1947), S. 6 
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Technische Dynamik. 2. Auflage. Von C. B. Biez—ENo und R. GRAMMEL 
(Springer-Verlag, Berlin 1953). Band I: XII + 700 S., 413 Abb., 3 Tafeln; 
DM 66.-. Band II: VIII + 452 S., 315 Abb., 3 Tafeln; DM 44.-. 

Die moderne Technik stellt an den konstruktiv tatigen Ingenieur immer 
héhere Anforderungen. Die «Sicherheitsspanne» zwischen Betriebs- und Bruch- 
beanspruchung schrumpft zusammen, so dass naherungsweise Berechnung nicht 
mehr gentigt und die héhere Festigkeitslehre herangezogen werden muss, um die 
auftretenden Spannungen mit geniigender Genauigkeit im voraus zu berechnen. 
Es ist also zweifellos ein Bediirfnis nach einem Handbuch iiber schwierigere Teil- 
gebiete der Mechanik vorhanden, und damit erklart sich der Erfolg, welcher der 
im Jahre 1939 erschienenen Technischen Dynamik beschieden ist. 

Die zweite Auflage dieses Werkes liegt nun vor. Es handelt sich hier nicht 
um ein Lehrbuch fiir den Studenten, sondern um eine Auswahl statischer und 
kinetischer Probleme, die sich durch die Methoden der theoretischen Mechanik 
untersuchen lassen. 

Im ersten Band werden zundchst die Grundlagen der Elastomechanik behan- 
delt, wobei mehr als 100 Seiten den verschiedenen Lésungsmethoden gewidmet 
sind. Als wertvolle Erginzung gegeniiber der ersten Auflage ist die Darstellung 
der Southwellschen Relaxationsmethode hervorzuheben sowie die allerdings 
etwas stiefmiitterlich behandelte Beschreibung der Dehnungsme®8streifen-(Strain- 
Gauges-) Methode. 

Der Hauptwert des Buches liegt zweifellos im zweiten Abschnitt, der verschie- 
dene der Technik entnommene Einzelprobleme behandelt. Dieser Teil ist eine 
wahre Fundgrube fiir den auf dem Gebiet der héheren Festigkeitslehre tatigen 
Ingenieur. Nicht nur sind die mathematischen Ableitungen von einer muster- 
giiltigen Klarheit, sondern die vielen numerischen Beispiele sowie die zahlreichen 
Kurventafeln erleichtern dem mathematisch wenig geschulten Ingenieur die 
Durchfiihrung der Berechnungen. 

Stab und Welle, Feder und Ring, Platte und Schale bilden die Hauptabschnitte 
dieses Bandes, wobei nicht nur Gleichgewichtszustande betrachtet werden; auch 
Ausweichprobleme, Knickung und Umkippung werden sehr ausfiihrlich, zum 
Teil nach neuen Gesichtspunkten behandelt. 

Der zweite Band ist den Fragen der Belastung durch Massenkrafte gewidmet, 
insbesondere den rotierenden Scheiben. Die Berechnung der Turbinenschaufeln 
bildet den Gegenstand eines Kapitels, wobei die Biegungsschwingungen beson- 
ders ausfiihrlich beriicksichtigt werden. Leider fehlt jeder Hinweis auf Torsions- 
schwingungen sowie auf komplexe Schwingungsformen, obwohl gerade letztere 
dem Praktiker viel Kopfzerbrechen verursachen, 

Der Abschnitt iiber kritische Drehzahlen ist eine mustergtiltige Darstellung 
dieses klassischen Kapitels der technischen Dynamik. Der Referent vermisst hier 
einen Abschnitt iiber Theorie und Praxis des Auswuchtens sowie einen Hinweis 
auf selbsterregte Schwingungen rotierender Wellen. Die zweite Halfte des zweiten 
Bandes behandelt ausfiihrlich Probleme der Kolbenmaschinen, wie Massenaus- 
gleich, Leistungsausgleich und Drehschwingungen. 

Ein gutes Inhalts- und Sachverzeichnis erleichtern den Gebrauch des vom 
Springer-Verlag wie iiblich ausgezeichnet gedruckten Buches, das dem theoretisch 
geschulten Ingenieur bestens empfohlen wird. P. de Halley 
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: Baustatik. Volume 1, 2e édition. Par Fritz Sttssi1 (Verlag Birkhauser, Bale 
(1953). 370 p., 336 fig.; 40 fr.s 05. 
Rappelons que cet ouvrage est la reproduction, sous une forme élargie, du 
cours de «Statique des constructions I» donné par le professeur STUssI aux 
éléves de la Section du Génie civil de l’Ecole polytechnique fédérale. L’auteur y 
expose avant tout les théories et les méthodes les plus récentes permettant de 
-calculer les tensions et les déformations des systémes statiquement déterminés. Il 
'y donne aussi une introduction 4 la théorie de la stabilité des constructions. La 
derniére partie est consacrée a l’équilibre des fils. 

La seconde édition de ce livre, tout en étant une mise au point soignée de 
louvrage, contient également un certain nombre d’améliorations et de complé- 
ments utiles. Signalons, en particulier, les modifications apportées au chapitre 

‘sur les arcs & trois articulations (III, 3), 4 celui sur les moments d’inertie (VEEL) 
et ala IXe partie concernant |’étude de la stabilité. 

Cet ouvrage se distingue non seulement par une excellente présentation, mais 
aussi par ses qualités d’exactitude et de clarté, qui le recommandent vivement 
aux étudiants ainsi qu’aux ingénieurs de la pratique. Henry Favre 


Proceedings of the Eastern Joint Computer Conference, gehalten vom 
8. bis 10, Dezember in Washington, D.C. Herausgegeben durch The Institute of 
Radio Engineers, Inc., New York. 125 S., 64 Fig.; $3.00. 

Die Broschiire gibt 25 Vortrage wieder, welche tiber das Sammelthema «In- 
formation Processing Systems—Reliability and Requirements» gehalten wurden. 
-Dementsprechend werden hauptsachlich Betriebserfahrungen mit bestehenden 
elektronischen Rechenanlagen sowie mit Einzelteilen gegeben, und es sind Ver- 
fahren diskutiert, die die heute noch unbefriedigende Betriebssicherheit verbes- 
sern sollen. Ferner werden einige Schaltelemente sowie neue Anwendungsgebiete 
fiir solche Maschinen beschrieben. — Fiir den Praktiker ist dieses Protokoll mit 
seinen ausfiihrlichen, zahlenmassigen Angaben von grossem Wert. A. P. Speiser 


Crystal Growth and Dislocations. Von A. R. Verma (Butterworth Scientific 
Publications, London 1953). XII + 182 S., 76 Pigt5 30s. 
Auf Grund der morphologischen Betrachtung von Kristallen war es schon 
lange bekannt, dass die Hauptwachstumsflachen in der Regel niedrig indizierbar 
‘sind. W. Kosset und I. N. Stranski haben vor etwa fiinfundzwanzig Jahren 
zeigen kénnen, dass diese Erscheinung energetisch leicht verstandlich ist. Wenn 
man namlich die Anlagerung von neuem Material an eine bereits bestehende 
Kristallflache Schritt fiir Schritt verfolgt, so zeigt es sich, dass der grdésste Ener- 
giegewinn dann erfolgt, wenn sich an irgendwelche zufallig vorhandene ober- 
flachliche Hécker (insbesondere monoatomare Schichtstiicke) seitlich Material 
anlagert, die Oberflachenerhebungen sich somit ausbreiten, bis sie véllig ver- 
schwinden. Demgegeniiber ist es energetisch sehr unwahrscheinlich, dass sich auf 
idealen Wachstumsflachen neue Keime solcher Schichten bilden, die seitlich 
wiederum wachsen kénnten. Man konnte deshalb bis vor kurzem nicht verstehen, 
wieso Kristalle mit glatten Grenzflachen trotzdem weiterwachsen kénnen. Durch 
einen Gedanken von F. C. FRANK hat nun die Wachstumsforschung einen starken 
Auftrieb erhalten, indem er 1949 die Hypothese aufstellte, dass immer dann, 
wenn eine Schraubenversetzung (bestimmte Gitterst6rung). an die Oberflache 
eines Kristalls reicht, dort eine beim Wachstum nie verschwindende Stufe ent- 
steht, an welcher sich nach dem Kossel-Stranski-Mechanismus standig neues 
Material anlagern kann. FRANK zeigte, dass um den Ausstichpunkt einer solchen 
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me 


Versetzung eine monoatomare Wachstumsschicht sich spiralenférmig nach aussen — 


bewegt, die zu einer dusserst flachen Wachstumspyramide fiihrt. Seit 1950 sind 


nun mittels elektronenmikroskopischer Untersuchungen und mit dem Phasen- 
kontrastmikroskop, dem Interferometer und anderen optischen Hilfsmitteln an 
zahlreichen Kristallarten tatsachlich derartige Wachstumsspiralen nachgewiesen — 
-worden, so dass heute kein Zweifel mehr dariiber bestehen kann, dass der Frank- ‘ 


sche Spiralwachstumsmechanismus tatsdchlich ablaufen kann. 


A. R. VERMA hat insbesondere an Siliziumkarbid zahlreiche Varianten des — 


| 
| 


Spiralwachstums eingehend untersucht und berichtet nun in der vorliegenden 


Monographie an Hand sehr guter Abbildungen iiber diese Phanomene, ihre Er- 
klarung und die angewandten Untersuchungsmethoden. Er bezieht dabei alle bis 
1953 bekanntgewordenen Falle des Spiralwachstums mit in seine Arbeit ein und 
gibt dadurch einen ausserordentlich eindriicklichen Uberblick iiber dieses neue 
Forschungsgebiet, das zwischen Realbauforschung und Kristallmorphologie ver- 


mittelt. Das Biichlein ist daher jedermann sehr zu empfehlen, der sich mit diesen — 


Gebieten beschaftigt. Der sehr leicht versténdliche Text macht es jedoch auch 


_ jedem Nichtfachmann méglich, sich iiber die moderne Wachstumstheorie zu _ 
_ orientieren. W. Epprecht — 


Journal of the Association for Computing Machinery. Herausgegeben 
durch die Association for Computing Machinery (Prasident S. B. WILLIAMS). 
4 Hefte jahrlich. 

Die Zeitschrift ist als Fachorgan auf dem Gebiet der elektronischen Rechen- 
maschinen gedacht. In einzelnen Artikeln soll hauptsachlich tiber die Fortschritte 
in der Technik der programmgesteuerten Maschinen berichtet und die einzelnen 
Neukonstruktionen beschrieben werden. Auch die mathematische Organisation 
und die Anwendung derartiger Gerate gehért in den Rahmen der Zeitschrift. 


Inhalt von Nr.1: S. B. Wixtrams, The Association for Computing Machinery. 


— J. W. Backus, The IBM 701 Speedcoding System. — R. T. WisEMan, Life 
Insurance Premium Billing and Combined Operations by Electronic Equipment. — 
F. E. Hamitton and E.C. Kusre, The IBM Magnetic Drum Calculator Type 
650.— H. Jacoss, Jr., Equipment Reliability as Applied to Analogue Computers. — 
C.M. Epwarps, Survey of Analogue Multiplication Schemes. - RICHMOND PERLEY, 
Automatic Strain-Gage and Thermocouple Recording on Punched Cards. — News 
and Reports. Supplement: ONR Digital Computer Newsletter. E. Stiefel 


Planen und Auswerten von Versuchen. Von A. LINDER (Verlag Birkhau- 
Ser, Basel 1953), 182 S.; Fr, 21,—. 

Das Buch ist in erster Linie fiir Biologen und auch fiir Ingenieure bestimmt, 
die Versuche statistisch auswerten miissen. Es enthalt die Darstellung der ein- 
facheren Methoden der statistischen Planung, und zwar durch genaue Beschrei- 
bung durchgefiihrter Experimente und ihrer Auswertung. Das Buch verlangt 
praktisch keine Vorkenntnisse in der mathematischen Statistik. Am Schlusse 
sind einige mathematische Anmerkungen beigefiigt. 

Es handelt sich um das erste in deutscher Sprache geschriebene Buch iiber 
statistische Planung. Mit Riicksicht auf die einfache, zweckmAdssige und klare 
Darstellung und die vielen Beispiele aus der Praxis zweifeln wir nicht daran, dass 
es vor allem bei den Biologen, Medizinern, aber auch bei einem weitern Kreis 
von Wissenschaftern grossen Anklang finden wird. Tatsachlich stellt es eine wert- 
volle Einfiihrung und Anleitung auf dem Gebiet der statistischen Planung dar, 
die auch fiir experimentell arbeitende Physiker, Ingenieure und andere von 
Interesse sein diirfte. W. Saxer 
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Band II 
Anwendungen der Laplace -Transformation 
1. Abteilung 


436 Seiten mit 48 Figuren 
Ganzleinen Fr./DM 56.15, broschiert Fr./DM 52.- 


Aus dem Inhalt: 
Asymptotische Entwicklungen — Konvergente Entwicklungen — 
Gewohnliche Differentialgleichungen. 


Wahrend der erste Band die theoretischen Grundlagen der Laplace-Trans- 
formation zum Gegenstand hat, behandelt der vorliegende zweite und der 
nachfolgende dritte Band die Anwendungen, wobeies sich natiirlich nicht nur 
um sogenannte «angewandte Mathematik», sondern um die verschiedensten 
Gebiete der reinen und angewandten Mathematik handelt, in welche die 
Laplace-Transformation als Hilfsmittel eingreift. 
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Blaschke, Geometrie der Waben. Von Professor Dr. Wilhelm Blaschke, Hamburg. 108 Sei- 
ten mit 42 Figuren. Broschiert Fr./DM 15.25. Sammlung «Elemente der Mathema- 
tik vom héheren Standpunkt aus», Band 4. : 


Doetsch, Handbuch der Laplace-Transformation. In drei Banden. Von Professor Dr. Gu-— 
stav Doetsch, Freiburg. Band II: Die Anwendungen der Laplace-Transformation. 
1. Abteilung. 436 Seiten mit 48 Figuren. Ganzleinen Fr./DM 56.15. Broschiert 
Fr./DM 52.-. ; 


_ Haack, Elementare Differentialgeometrie. Yon Professor Dr. Wolfgang Haack, Berlin. 


239 Seiten mit 25 Figuren. Ganzleinen Fr./DM 22.—. Broschiert Fr./DM 18.70. 


Schlafli, Gesammelte mathematische Abhandlungen. In drei Banden. Herausgegeben im 
Auftrag der Schweizerischen Mathematischen Gesellschaft vom Komitee Steiner- 
Sehlafli der Schweizerischen Naturforschenden Gesellschaft, erganzt und erlautert 
durch Nachworte der Herausgeber. Band III: 402 Seiten mit 28 Figuren. Ganz- 
leinen Fr./DM 59.30. 


Speiser, Die Theorie der Gruppen von endlicher Ordnung. Mit Anwendungen auf alge- 
_ braische Zahlen und Gleichungen sowie auf die Krystallographie. Von Professor 
Dr. Andreas Speiser, Basel. Vierte berichtigte und erweiterte Auflage. XI und 

271 Seiten mit 43 Figuren, einer Farbtafel und einem Anhang. Im Druck. - 


Weyl, Selecta Herman Weyl. Herausgegeben zu seinem siebzigsten Geburtstag von der 
ETH. Ziirich und vom Institute for Advanced Study in Princeton. Rund 600 Seiten 
mit Figuren. Ganzleinen Fr./DM 48.90. 


Wissenschajt und Kultur 


Speiser, Die geistige-Arbeit. Von Professor Dr. Andreas Speiser, Basel. 216 Seiten. Fr./DM 
19.25. Ganzleinen. Sammlung «Wissenschaft und Kultur», Band 9. 


Van der Waerden, Erwachende Wissenschaft. Von Prof. Dr. Bartel van der Waerden, 
Ziirich. Ubersetzung aus dem Hollandischen von Helga Habicht. Mit Zusitzen vom 
Autor versehen. Im Druck. Sammlung «Wissenschaft und Kultur», Band 8. 


Weyl, Symmetrie. Von Professor Dr. Hermann Weyl, Ziirich und Princeton. Ubersetzung 
aus dem Englischen von Professor Dr. Lulu Bechtolsheim. 157 Seiten mit 72 Ab- 
bildungen. Im Druck. Sammlung «Wissenschaft und Kultur», Band 11. 


Weitzenbéck, Der vierdimensionale Raum. Von Professor Dr. Roland W. Weitzenbick, 


Zelhem. 224 Seiten mit 52 Figuren. Fr./DM 19.55, Im Druck. Sammlung «Wissen- 
schaft und Kultur», Band 10. 
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